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PREFACE 

The  computer  program  described  in  this  report  was  developed  as 
a  part  of  a  project  which  was  a  continuation  of  work  performed 
by  the  Kelly  Scientific  Corporation  on  Forecasting  of  Heavy 
Loading  Patterns  on  Highway  Bridges  (1) .   The  computer  dynamic 
stress  analysis  program  developed  in  that  work,  BRGSTRS,  is 
described  in  this  report. 

The  cooperation  of  Dr.  Conrad  Heinz  of  the  University  of 
Maryland  in  furnishing  strain  gage  recordings  and  correlated 
truck  data  on  bridges  in  Maryland  is  gratefully  acknowledged. 
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INTRODUCTION 

The  work  performed  and  the  results  obtained  under  Contract 
FH-11-7904  to  the  Federal  Highway  Administration,  Structures 
and  Applied  Mechanics  Division,  for  a  study  on  "Prediction  of 
Loadings  on  Highway  Bridges — Phase  II",  is  described  in  four 
separate  reports,  FHWA-RD-73-42 ,  43,  44  and  45. 

A  primary  objective  of  this  study  was  to  extend  the  traffic 
simulator  computer  program  developed  by  the  Kelly  Scientific 
Corporation  to  a  useable  engineering  tool  capable  of  generating 

bridge  loads.   It  was  further  the  purpose  of  the  study  to  develop 
a  finite  element  stress  analysis  program,  including  dynamic 
effects  due  to  the  live  load,  which  directly  interfaced  with 
the  load  generator,  produce  several  stress  histograms  for 
various  bridge  and  traffic  configurations,  and  to  develop  and 
implement  an  alternative  method  based  upon  analytic  methods 
rather  than  traffic  simulation. 

Report  FHWA-RD-73-42  describes  the  work  performed  in  the 
study.   This  includes  the  background  which  led  to  the  work 
performed,  a  description  of  the  work  performed  and  problems 
encountered  in  revising  BRIGLD1 .   The  results  of  sensitivity 
testing  of  BRIGLD1,  development  of  the  stress  program,  genera- 
ted histograms,  description  of  the  analytic  methods  investigation, 
and  the  conclusions  and  recommendations  are  also  included. 

Report  FHWA-RD-73-4  3   is  the  Users  Manual  which  provides 
utilization  instructions,  data  preparation  instructions,  output 
and  variables  definitions  and  a  description  of  the  BRIGLDl 
computer  program.   The  description  includes  a  narrative 
description  section,  flow  charts  and  program  listings  on  the 
main  line  program  and  each  subroutine.   The  use  of  this 
program  for  stress  range  prediction  is  illustrated  in  Figure  1. 
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Figure  1.   Traffic  Simulation  Based  System 


Report  FHWA-RD-73-44   is  the  Users  Manual  for  the  dynamic 
stress  analysis  computer  program,  BRGSTRS.   This  report 
contains  the  same  descriptive  type  matter  as  indicated  above 
for  BRIGLD1. 

Report  FHWA-RD-73-45   contains  the  Users  Manuals  for  two 
computer  programs  which  operate  as  a  system  with  BRGSTRS. 
The  first  is  the  Synthetic  Load  Generator,  SYNGEN,  which 
generates  single  axle  loads  for  the  dynamic  stress  analysis 
program,  which  in  turn  generates  a  stress  signature  curve,  trace, 
for  each  defined  axle  load.   The  second  is  the  histogram 
computer  program,  HISGEN,  which  generates  long-term  stress 
range  histograms  from  the  synthetic  single  axle  stress  trace 
data  generated  by  the  dynamic  stress  analysis  program,  as  shown 
in  Figure  2.   This  is  accomplished  by  first  forming  composite 
truck  stress  traces  for  a  given  truck  population  from  the 
single  axle  data.   Then,  forming  composite  truck  platoon  stress 
traces  for  a  given  platoon  population.   Long-term  effects  are 
estimated  from  traffic  density  estimates  and  the  estimated 
incidence  of  each  platoon  configuration.   The  information 
contained  in  this  report,  for  both  SYNGEN  and  HISGEN,  is  of  the 
same  form  as  described  above  for  BRIGLDl. 
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Figure  2.   Axle  Signature  Based  System 
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OPERATING  INSTRUCTIONS 

Operating  instructions  for  the  Bridge  Dynamic  Stress  Analysis 
Program  are  divided  into  four  sections,  (1)  Data  Input,  which 
includes  descriptions  of  all  data,  default  parameters,  and  data 
input  coding  forms,  (2)  Data  Output  with  examples  of  printed 
output  and  a  description  of  the  Bridge  Load  Data  Block,  (3)  Deck 
Set  Up  description,  and  (4)  a  Sample  Case  run. 

Prior  to  using  BRGSTRS,  if  a  user  utilizes  BRIGLD1  to  generate 
traffic  data  for  BRGSTRS,  he  should  verify  the  following: 

1.  That  he  is  getting  the  correct  distribution  of  trucks 
by  type  for  the  expected  traffic. 

2.  That  he  is  getting  the  correct  density  of  trucks,  for 
the  simulated  period  of  time,  for  the  expected  traffic. 

3.  That  he  is  getting  truck  weights  that  correlate  with 
the  expected  real  truck  weights. 

Subsequent  to  verifying  the  above,  he  can  then  use  BRGSTRS  on 
the  generated  load  data  stored  on  magnetic  tape. 


Data  Input 

Data  Input  is  divided  into  program  header,  single  element  data, 
and  tabular  data.   Tabular  data  consists  of  Longitudinal  Beam 
Parameters,  Lateral  Beam  Parameters,  and  Output  Specifications. 

The  program  header  consists  of: 

Card   Col    Item    Description 

1    1-2    ITAPE1    Bridge  Load  Data  input  device  number 
1    3-4    ITAPE2    Point  Stress  Trace  output  device  number 
1    5-72   HEAD     6  8  character  run  header  for  printout  and 

output  on  Point  Stress  Trace  tape 

The  following  figures  illustrate  some  of  the  data  elements  to  be 
input. 
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Figure  3.   Bridge  Deck 
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Figure    4.      Diaphragm   Section 
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Figure  5.   Steel  Girder 
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Figure    7.      Simple    Span 
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Figure    8.      Continuous    Beam 


Single  Element  Data  -  Single  element  data  items  are  input 
under  namelist  INPUTA.   Only  those  elements  with  no  default 
values  need  be  entered  and  the  order  of  entry  within  the 
namelist  is  arbitrary.   Data  Form  1  may  be  used  for  the 
NAMELIST  input  (see  Figure  9).   The  distances  CRUD  and  CRLD 
need  only  be  entered  for  concrete  lateral  beams  and  the 
distances  C2LRU  and  C2LRL  need  not  be  entered  if  effective 
longitudinal  beam  parameters  are  entered. 

The  first  card  $DATA  and  last  card  $END  must  be  provided  even 
if  there  is  no  data  to  be  entered.   A  comma  must  follow  each 
data  element  except  the  last.   The  following  elements  are 
input  under  namelist  DATA: 


ITEM 

DEFAULT 

XLB 

NL 

2 

YBCL 

EC 

ES 

NUC 

NUS 

DT 

Calcula 

THRSH 


UNITS    MAX. 


Ft. 


Ft. 


Lb/in' 
Lb/in' 


300 


Psi 


DTL 


DESCRIPTION 

Length  of  the  bridge. 
Number  of  lanes. 
Distance  to  center 
line  of  the  bridge. 
Young's  Modulus  for  concrete, 
Young's  Modulus  for  steel. 
Poisson's  Ratio  for  concrete 
Poisson's  Ratio  for  steel. 
Integration  interval  for 
dynamic  calculations.  This 
value  is  calculated  if  it 
is  not  input. 

Stress  threshold  value  for 
steady  state  and  end  of  a 
stress  event. 
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ITEM 

DEFAULT 

UNITS 

WEIGHT 

Lbs  . 

DR 

In. 

CRUD 


CRLD 


NB 

NSUP 
C2LRU 


C2LRL 


DRT 


NLRU 


NLRL 


In. 


In. 


In 


In 


In 


MAX.   DESCRIPTION 

Weight  of  the  bridge. 
Diameter  of  the  reinforcing 
steel  rod  used  in  the  deck. 
Distance  from  upper  deck 
surface  to  center  of  upper 
layer  of  reinforcing  steel 
rods. 

Distance  from  upper  deck 
surface  to  center  of  lower 
layer  of  reinforcing  steel 
rods . 

10  Number  of  main  longitudinal 
beams . 

11  Number  of  main  beam  supports. 
Distance  from  the  top  of  the 
deck  to  the  neutral  axis  of 
the  upper  layer  of  reinforcing 
steel  rods ,  concrete  diaphram. 
Distance  from  the  top  of  the 
deck  to  the  neutral  axis  of 
lower  layer  of  reinforcing  steel 
rods,  concrete  diaphram. 
Diameter  of  reinforcing  steel 
rods  in  transverse  direction. 
Number  of  reinforcing  steel 
rods  in  the  upper  layer  of 

the  lateral  cross  section  of 
the  deck  (DELX  Wide) . 
Number  of  reinforcing  steel 
rods  in  the  lower  layer  of 
the  lateral  cross  section 
of  the  deck  (DELX  Wide) . 
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ITEM     DEFAULT 


CLRU 


UNITS 
In. 


CLRL 


In. 


ILSC 


Tl 


In. 


IL 


SIGTHS 

300. 

psi 

DSIG 

3000. 

psi 

NDBUG 

0 

MAX    DESCRIPTION 

Distance  from  the  deck  surface 
to  the  neutral  axis  of  the  upper 
layer  of  reinforcing  steel 
rods  in  the  lateral  cross 
section  of  the  deck. 
Distance  from  the  deck  surface 
to  the  neutral  axis  of  the  lower 
layer  of  reinforcing  steel  rods 
in  the  lateral  cross  section 
of  the  deck. 
Steel  or  concrete 
diaphram  indicator. 

0  =  steel,  Not  0  =  concrete 

Mean  thickness  of  the  bridge 

deck. 

Read  control:   if  =  0,  read 
diaphram  parameters,  if 

not  =  0,  read  effective 
diaphram  parameters. 

Stress  peak  decay  threshold. 
Histogram  class  interval  size. 
Debug  printout  selection  (10 
element  array) .   If  =  0  no 
printout.   If  not  =  0  print 
intermediate  values  for  the 
following  programs: 


Element 

Program 

1 

MAIN/GRDNET 

2 

MODULI 

3 

STIFF/ASEMBL 

4 

DMOMNT 

5 

PNTLOD 

6 

GRDLOD 

7 

LODDIS 

8 

BROUT 

9 

BRGDYN 

10 

INCRMT 

12 


CARD 
3 


COL. 
1-8 

9-16 
17-24 


ITEM 
XCS(l) 

XCS(2) 
XCS(3) 


DESCRIPTION 

Distance  to  the  coverplate, 

see  Figure  8;  ft.,  F8.3. 

Distance  to  the  end  of  the 

coverplate,  see  Figure  8;  ft., 

F8.3. 

Distance  to  the  end  of  the  beam, 

see  Figure  8;  ft.,  F8.3. 


Beams  may  be  defined  in  three  different  manners,  (1)  grossly  by 
specifying  the  total  parameters  (cards  4A,  Figure  10) ,  or  (2) 
to  be  calculated  in  detail  as  a  steel  beam  (cards  4B,  Figure  10) , 
or  (3)  to  be  calculated  in  detail  as  a  concrete  beam  (cards  4C 
and  4D,  Figure  10).   If  coverplates  are  included,  one  card  each 
of  4A,  4B  or  4D  must  be  defined  for  each  cross  section  variation. 

Effective  Longitudinal  beam  parameters  are: 


CARD 
4A 

4A 

4A 

4A 


COL. 

1-10 

11-20 
21-28 
29-38 


ITEM 
INRT 

EBAR 


DESCRIPTION 

.  4 
Effective  Moment  of  inertia  in 

E10.3. 

Effective  Young's  Modulus  lb/in' 

E10.3. 

Effective  neutral  axis  distance 

from  deck  surface,  in.,  F8.3. 

Steel  or  Concrete  Modulus 

depending  on  composition  of 

2 

girder  lb/in  ,  E10.3. 


Longitudinal  Steel  Beam  parameters  are 


CARD 
4B 


COL. 
1-8 


ITEM 
AI 


DESCRIPTION 


Girder  cross-sectional  area, 


in  ,  F8.3.' 


13 


o 

z 

~ 

o 

-O 

». 

- 

<- 

m 

•* 

^ 

- 

- 

"""    I 

n 

i 

1 

!     1     ! 

i         !         '         i 

w 

vjn 

<*> 

•O 
ro 

r\> 

ro 

\*> 

v; 

w 

s 

rj 
J* 

1: 
ft 
ft 

•O 

>c 

o 

r:    - 

£ 

cr 

03" 

"2 

f: 

Sd~IH 

h3  jU  jZ  !Z 
S   rJ  ic!  i  Si 

jW  In   ihrl  !M   !F  r 

„ 

% 

ii  [W 

* 

IT  \Ui  |0 
III     III 

II 

II       JO     |ll 

u  h 

i 

>— 1 

1-1 

1> 

Q.LY' 

jr 

-— - 

i'^rS 
14    PJ 

'U 

" 

ir 

EG        [      i 

j  ii 

> 

H 

n 

II 



— 

n 

•Si       I      ! 

1— 1 

i 

!       !       ^       !                    !       ! 

— 

.  __j  _ 

I      i      ;            i 

_.}.._ 

-,0Q 

— 

— 

| 

!.      1    .._!_    1 

.-.-4 

— 

— 

— 

! 

1        !        ! 

— 

] 

.  ...L    L. 

1        1 

1 

i 

i 

"HI 

>:  o 

1      j 

p 

..  Oi 

i  " 

^  h 

i 

1 

i 

■ 

~ 

"  T" 

1 

a 

1 



— 

..J 

1 

! 

(o 

1 

K 

.... 

Si 

i 

— 

— 

<* 

..  .J.    . 

.   __j. 

-°. 

! 

| 

— 

&j 

| 

2 

►3 

t 

' 

1 

8 

i 

•-' 

c 

r    j 

:zl  t 

H 

i    t 

__.!_■  |._ 

— 

— 

!i! 

fl> 

"  i ■ 

'  J 

VD 

— 

— 

. 



..... 

■   r~ 
—  i — 

j 

j.  .. 

— 

\J% 

• 

— 

._... 

.... 

i 

i 

1 

J*£ 

— 

i 

! 

w 

— 

i 

'    j 

:.•; 

O 
ft 

— 

— 

— 

— 

w 

! 

— 

ft 

OJ 

— 

— - 

i 

..._. 

— 



— 



i: 

— 

— 

— 

To 

0 

— 

....... 



— 

__!_.... 

i^ 

— 

. — 



— 

— 

— 

! 

"1 

3; 

3 

— 

cS 

1 

ft 

H 

^_ 



_  .. 

........ 

........ 

I 

1 

— 

1 

ii 

— 

_..._ 

.... 

-|- |- 

ft 

. — 

— 1— 

1 
1 " 

i 

'<?> 

..... 

— 

... 

I 

M' 

.. 

. — 



— 

— 

_... 

...... 

'/- 

1 

i 

Z.1..Z 

I 

I 
" r~~ 

— 

<j\ 

:: 

— - 

— 

._._ 

— 

"T 

.._. 

E 

« 

_.._.. 

■x'! 

— 

— 



-— t - 

i 
._1 ._ 

r 

z'."iz: 

... 

O 

• — 

— - 

_... 

.>- 

-■ 

— 

.....  | 
ztz; 

■_-- 

_  j 

._  j _, 



i 

-"  i 

z".iz: 

i 
"  "I"-'" 

£ 

— 

._.._ 

_.. 

& 

-— 

.... 

ro 

-— 

— 

1 

- 

..... 

! 

| I I 1 ; 

.... . 

^ 

... 

- 

— - 

..._. 

— 

- 

: 

i 

"Hz: 

i     ..  |  -..-..- 

....  j...... 

dp. 
._ 

8 

14 


a 

k 

rO 

fc 

■— 

— 

— 

-o 

o~- 



s 

— 

c 

PJ 

— 

£ 

o 

— 

— 

s 

O 

s 



3 

(Xi 

-o 

s 

o- 

— 



u-\ 

a 

~ 

R 

t/\ 

K 

i 

R 

5 

3 

^ 

1 

8 

\/-. 

5 

-T 

fVJ 

T 

Sf 

S^ 

r^ 

S 

1 

_|_ 

£ 

_J 1 

r 

1 

i*» 

! 

CM 

i     |     i 

K 

i               j 

R 

I 

!     i 

r*i 

1 

CI 

i 

S 

i 

i    « i 

S 

I 

J 

SQ 

i 

" 

L 

T 

i 

P3 

I 

i 

"  L 

CM 

I 

I 

i 

CNI 

i 

■ 

i       !     - 

a 

1 

!        i 

2 

1 

^_ 

2 

1       1       i 

1 

^ 

1     !     i 

1       1     - 

•o 

1     '     ' 

1 

2 

1     ' 

1       1     - 

~ 

__     ■  |     i 

[ 

2 

1     i     i 

. 

a 

i 

~ 

1       !       I       ]       1       I 

- 

1       1       |               1 

o> 

00 

!  "  I 

*■*" 

I             ! 

i 

col       i     ii 

•o 

it 

II  1       ll|       II     j 

i 

Xi  ill  u 

u> 

oi  -Ji  n\ 

H1  tli  3 

A 

""7 

pq|  PhI  co i  ii  i 

ii 

r,l|  > — ■■  5! 

s 

ti 

~  K 

-5  ^i 

J         ~[ 

,,] 

»— 1 1 

C0|   1-1 

1 

n 

ft 

p— i  ,    E-r     i— i 

CO 

-Qi~Z  <& 

1       i       !       I 

— I 

1      1 

u 

0) 

id 

0) 

U 

<u 
C 
d 

U 


o 

Oh 

+-> 

o 
fll 
a, 
d 
co 

CU 
co 


-O 

w-v 

i—i 

r^ 

^-' 

>o 

"O 

*C> 

S 

•o 

3 

*2 

^ 

o 

• 

oo 

a 

ir\ 

^ 

1 

— 

5C 

ss 

o- 

g 

s 

—j. 

rj 

•"» 

S3 

• 

S\ 

co 

» 

3 

<■•"» 

fVj 



'*■* 

p; 

o- 

zz 

vc 

Jp 

• 

u"\ 

^ 



r, 

,-v, 

OJ 

R 

— • 

■ 

i2 

,-_ 

•O 

^ 

zi 

a 

~ 

— - 

" 

— 

*** 

r- 

•o 

-J-. 



— 

— 

T3 

d 

c 
o 

u 


s 

o 


Q 

0) 
h 

d 

60 


15 


CQ 


CARD 

COL. 

ITEM 

1A 

1-10 

YCL(l) 

1A 

11-20 

YCL(2) 

IB 

XSUP(l) 

Tabular  data  input  follows  the  NAMELIST  data  input  as  follows: 

Lane  Center  Line  Positions 

DESCRIPTION 

Center  line  position  of  Right  Lane 

input,  ft.,  Format  F10.4. 

Center  line  position  of  Left  Lane 

input,  ft.,  Format  F10.4. 

Position  of  first  Bridge  Support 

(normally  =  0.0) ,  Format  11F7.2. 

Longitudinal  Beam  Data 

The  following  set  of  data  must  be  repeated  for  each  effective 
longitudinal  beam  up  to  10  beams.   (See  Figure  10.) 

CARD      COL.      ITEM        DESCRIPTION 

2       1-3       NRU        Number  of  upper  reinforcing  steel 

rods  in  the  concrete  deck  -  13. 

2       4-6       NRL        Number  of  lower  reinforcing  steel 

rods  in  the  concrete  deck  -  13. 

2       7-16       IUR        Moment  of  inertia  of  a  steel  rod  in 

.  4 
the  upper  part  of  the  beam  in  ,  E10.3. 

2      17-26       ILR        Moment  of  inertia  of  a  steel  rod  in 

4 
the  lower  part  of  the  beam  in  ,  E10.3. 

2       27-34      YB         Distance  to  center  of  beam, 

ft.,  F8.3. 

Effective  beam  width,  ft.,  F8.3. 
Number  of  variations  in  the  girder 
cross  section,  Max.  =  3  (w/wo  coverplates! 
Steel  or  concrete  girder  index 
Beam  data  index,  0  =  Read  all  beam 
parameters,  not  0  =  Read  effective 
beam  parameters . 
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4B        9-16      IGO 
4B        17-24      CI 


4B 


4B 


4B 


25-32 


CT 


33-40     CCP 


41-48     TCP 


Moment  of  inertia  of  the  nth  cross- 

4 
section  of  the  Jth  girder,  in   ,  F8.3 

Distance  to  the  neutral  axis  of  the 

Nth  cross  section  of  the  Jth  girder, 

in.,  F8.3 

Total  distance  from  deck  surface  to 

the  bottom  of  the  Nth  cross  section, 

girder  bottom,  in.,  F8.3 

Distance  to  the  neutral  axis  of  the 

bottom  cover  plate  of  the  Nth  cross 

section,  in.,  F8.3 

Thickness  of  plate  in  the  Nth  cross 

section  if  cover  plate  exists,  in., 

F8.3 


Longitudinal  Concrete  Beam  Parameters  are: 


CARD 


4C 


4C 


4C 


4C 


4C 


COL. 


1-8 


9-16 


17-24 


25-27 


28-30 


ITEM 


BG 


CG 


TG 


NGRU 


NGRL 


DESCRIPTION 

Width  of  the  Mth  section  of  the  Nth 
cross  section  of  the  girder,  in., 
F8.3 

Distance  to  the  neutral  axis  of  the 
Mth  section  of  the  Nth  cross  section 
of  the  girder,  in.,  F8.3 
Thickness  of  the  Mth  section  of  the 
Nth  cross  section  of  the  girder,  in. , 
F8.3 

Number  of  the  reinforcing  steel  rods 
in  the  upper  layer  of  the  Mth  section 
of  the  Nth  cross  section  of  the  gird- 
er, 13 

Number  of  reinforcing  steel  rods  in 
the  lower  layer  of  the  Mth  section 
of  the  Nth  cross  section  of  the  gird- 
er, 13 
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4C 


4C 


4C 


4D 


31-38 


CRU 


39-46 


CRL 


47-54 


IRGO 


1-8 


DRG 


Distance  to  the  neutral  axis  of 
the  upper  layer  of  reinforcing  steel 
rods  in  the  Mth  section  of  the  Nth 
cross  section  of  the  girder,  in., 
F8.3 

Distance  to  the  neutral  axis  of 
the  lower  layer  of  reinforcing  steel 
rods  in  the  Mth  section  of  the  Nth 
cross  section  of  the  girder,  in.,  • 
F8.3 

Total  Moment  of  inertia  of  the  rein- 
forcing steel  rods  of  the  Mth 

section  of  the  Nth  cross  section  of 

.  4 
girder,  in  ,  F8.3 

Diameter  of  the  reinforcing  steel 

rods  of  the  Nth  cross  section  of 

the  girder,  in,  F8.3. 


Diaphragm  Parameters  may  be  input  (cards  5B)  or  Effective 
Diaphragm  Parameters  may  be  used    (See  Figures  11  and  12.). 

Effective  diaphragm  Parameters  are; 


CARD 


COL. 


ITEM 


DESCRIPTION 


5A 


5A 


5A 


1-8 


9-18 


19-28 


DELX 


IT 


ET 


Grid  incremental  distance  per  span 
L,  ft. ,  F8.3 

Effective  Moment  of  inertia  per 

.  4 
DELX  width  in  each  span  L,  in  , 

E10.3 

Effective  Young's  Modulus  per 

2 
DELX  width  in  each  span,  lb/in  , 

E10.3 


Diaphragm  Parameters  are: 
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CARD 


COL. 


ITEM 


DESCRIPTION 


5B 


5B 


5B 


5B 


5B 


5B 


5B 


1-8 
9-16 

17-26 
27-29 

30-32 

33-40 
4  —  48 


DELX 


ALB 


ILB 


N2LRU 


N2LRL 


DELXP 


T2 


Grid  incremental  distance  per  span 
L,  See  Figure  4,  ft.,  F8.3. 

Area  of  steel  diaphragm 

2 

in  the  lateral  direction  in  , 

F8.3. 

Moment  of  inertia  of  steel  real 

4 
diaphram,  in  ,  E10 . 3 

Number  of  reinforcing  steel  rods  in 

the  upper  layer  of  a  concrete 

diaphram,  13 

Number  of  reinforcing  steel  rods  in 

the  lower  layer  of  a  concrete 

diaphram,  13 

Width  of  the  concrete 

diaphram,  in.,  F8.3 

Thickness  of  the  concrete  diaphragm, 

in. ,  F8.3 


Output  options  are  specified  on  cards  6  and  7  (see  Figure  13, 
Data  Form  4) .  One  card  7  must  appear  for  each  beam  whether  or 
not  points  are  specified.  Output  options  are: 


CARD 


COL 


ITEM 


I0PT1 


IOPT2 


IOPT3 


IOPT4 


DESCRIPTION 

Stress  graph  for  specified  points 
1  =  printout,  0  =  no  printout 
Stress  printout  for  specified  points 
1  =  printout,  0  =  no  printout 
Stress  histogram  for  specified  points 
1  =  printout,  0  =  no  printout 
Full  stress  printout  (not  recom- 
mended for  long  runs  because  of  the 
amount  of  print  generated)  1  =  print- 
out, 0  =  no  printout 


en 
£ 
O 

a 
O 
+-> 
d 

a 

O 


^ 
"i 


. — i 
£ 

0) 

PQ 

s 

0) 

pq 

CO 

£ 

CD 

pq 

£ 

0$ 

CD 

in 

£ 

CD 

pq 

s 

rd 

CD 

a 

cd 

pq 

00 

£ 

CD 

PQ 

£ 

CD 

pq 

o 

£ 

CD 

pq 

•  r 

CNj 

— 



CO  RJ 

O  2 

Oh  2 

— - 

— -- 

en 

a 
o 

Oh 

o 
u 

£ 

— 

5 

— 

— 

— 



— 

— ■ 

Ph  S 

... 

— 

— - 

— -•- 

....-•-. 

•_ 

_.._. 

1  \     ^> 

.2 

Oh    -pi 

— 



— ~- 



— 

— 

— 

-!     1 

£ 

o 


Q 


ro 
H 

CD 


vO 


24 


7 

2-9 

7 

10-17 

7 

18-25 

CARD     COL      ITEM     DESCRIPTION 

7     1       NPTS    Number  of  points  to  be  sampled  this  beam, 

up  to  3  points.  If  no  points  are  desired, 
enter  zero.  II 
PTS     First  sample  point  F8.3 
Second  sample  point  F8.3 
Third  sample  point   F8.3 

Sample  points  will  be  converted  to  the  nearest  grid  point  used 
for  calculation  within  the  program.  The  converted  sample  points 
will  be  printed  out  at  the  beginning  of  the  run. 

A  user  should  always  specify  sample  points  at  midspan  on  those 
beams  directly  under  the  two  right  hand  lanes.   Further,  when 
specifying  the  center  lines  of  the  two  allowed  lanes,  on  a 
bridge  of  more  than  two  lanes,  these  center  lines  should  be 
the  actual  center  lines  of  the  two  right  most  lanes,  i.e., 
where  the  heaviest  traffic  density  will  occur.   If  sample  points 
are  specified  only  at  locations  and  beams  where  light  loading 
will  occur  its  possible  that  no  output  will  be  generated  because 
of  the  stress  threshold  value  for  steady  state  and  for  the  stress 
range  threshold  value  for  determining  a  point  of  maxima.   A 
user  can  essentially  out  smart  himself  and  obtain  no  output. 
There  will  be  nothing  wrong  with  the  program,  only  the  user's 
input. 
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Bridge  Load  Data  Input 

Bridge  load  data  is  input  to  the  program  in  data  blocks  correspond- 
ing to  a  generation  integration  interval  DTL  and  Event  number 
IEVNTX.   The  Bridge  Load  Data  block  consists  of  the  following  ele- 
ments : 


Item 

Dimension 

SUMHRX 

1 

DTL 

1 

IEVNTX 

1 

NOAXLX 

1 

NTRUKX 

1 

LAST 

1 

LTYPE 

20 

WGTX 

20 

SPDX 

20 

LLANE 

20 

TIMEX 

20 

XPOSX 

50 

LNUMX 

50 

WEITX 

50 

DXPOSX 

50 

ACCLRX 


50 


Description 

Simulation  time 

Integration  interval 

Event  number 

Number  of  axle  loads  on-going  on 

the  bridge 

Number  of  trucks  this  event 

Logical  =  True  if  this  is  the 

last  block  of  this  event 

Truck  type 

Weight  of  truck 

Speed  of  truck 

Lane  of  truck 

Time  truck  went  on  the  bridge 

Axle  position 

Axle  lane 

Axle  weight 

Distance  axle  will  travel  during 

integration  interval 

Acceleration  of  axle 


The  truck  identification  data,  LTYPE,  WGTX,  SPDX.  LLANE,  and 
TIMEX,  was  captured  as  a  truck  went  on  the  bridge  and  does 
not  change  during  an  event  except  that  new  trucks  are 
added  to  the  list  as  they  occured  during  the  event.   The  axle 
load  data,  XPOSX,  LNUMX,  WEITX,  DXPOSX,  AND  ACCLRX  are  current 
values  for  each  integration  step. 
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Data  Output 

Output  from  the  Dynamic  Stress  Analysis  Program  consists  of: 

•  Printout  of  Run  Input  Parameters 

•  Optional  printout: 

1.  Stress  graph  for  specified  sample  points. 

2.  Stress  printout  for  specified  sample  points. 

3.  Stress  histogram  for  specified  sample  points. 

4.  Full  stress  printout. 

•Point  Stress  Trace  output  tape  -  if  device  number 
is  specified 

•  Debug  output 

Printout  of  Run  Input  Parameters  -  Run  Parameter  printout  depicted 
belonging  under  sample  case  consists  of  the  NAMELIST  DATA 
printout,  bridge  deck  specifications,  longitudinal  beam  data, 
lateral  beam  data,  output  options,  and  sample  points  on  each 
beam  as  input  and  as  converted  to  the  nearest  computational  grid 
point. 

The  stress  graphical  output  consists  of  a  plot  of  stress  vs. 
time  for  up  to  three  specified  points  on  each  beam.   All  points 
on  the  same  beam  will  appear  on  the  same  graph.   One  set  of 
plots  will  be  generated  for  each  event.   Three  distinct  symbols 
will  be  used  for  sample  points,  i.e.,  point  1  =  "*",  point  2  = 
"+",  point  3  =  ".". 

The  stress  printout  for  specified  sample  points  will  be  a  tabu- 
lar printout  of  the  minimum,  maximum  and  range  stress  data  for 
each  cample  point  by  event. 
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The  stress  histogram  will  consist  of  a  histogram  of  maximum  stress 
points  for  each  sample  point  broken  into  10%  class  intervals. 

The  full  stress  printout  will  consist  of  a  printout  for  each  time 
frame  as  follows : 

Time  (t) 

Span  #  (1) 

Beam  #  (1) 

Position:  xn  x~  x_ x 

12   3  n 

StreSS:  °B1  aB2  QB3 QBN 

Dynamic  Deflection:    <5   ^   •$ ° 

J  12   3  n 

Dynamic  Load  Factor:   Dlfn  Dlfn  Dlf_ Dlf 

2  12     3  n 

Beam  #  (2) 

* 

Beam  #  (J) 

Span  #  (2) 
Beam  #  (1) 

Span  (N) 

The  point  stress  output  tape  will  consist  of  the  following: 

Header: 

Nev,  W,  V,  NWI/  At,  BL,  Nfi/  Np(l),  .  ..,  Np (NB) 

Data  Block: 

0(1,1),  ...,  a  (NP(1) ,1),  a  (1,2),  ...,  a  (NP(2),2) ,  ... 

a  (1,NB) ,  ...,  a(NP(NB),  NB) 
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Each  data  block  is  repeated  for  each  time  point  in  the  stress 
calculations  for  event  Nev:   The  next  event  header  follows 
immediately  after  the  last  data  block  of  the  previous  event. 
Each  header  and  data  block  is  a  separate  record. 

Data  elements  are  defined  as  follows: 

W  Weight  of  axle  (real) 

V  Speed  of  axle  (real) 

At  Time  increment  (real) 

BL  Bridge  Length  (real) 

Nev  Event  Number  (integer) 

NNL  Lane  Number  (integer) 

N  Beam  Number  (integer) 

a 

Np(n)     Sample  points/Beam,  l^=n4=3     (integer) 

a,   ,   -,     Stress  value  l^nl^NP,  l4n24NB  (real) 
(nl ,nZ) 

Debug  output  will  be  generated  only  on  user  specification  and 
will  consist  of  all  pertiment  parameters  used  or  generated  by 
the  specified  subroutine.   Debug  output  of  a  subroutine  is 
specified  by  setting  the  corresponding  input  parameter  NDBUG(n) 
equal  to  1.   Program  subroutines  correspond  to  the  following 
NDBUG  elements: 


NDBUG  (1) 

GRDNET 

(2) 

MODULI 

(3) 

STIFF 

(4) 

DMOMNT 

(5) 

PNTLOD 

(6) 

GRDLOD 

(7) 

LODDIS 

(8) 

ASEMBL 

(9) 

BRGDYN 

(10) 

INCRMT 
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Default  setting  of  the  NDBUG  array  is  equal  to  zero  specifying 
no  debug  printout . 

Sample  output  may  be  found  in  the  section  Sample  Case  below. 

It  should  be  noted  that  the  output  of  the  input  data  includes : 

1.  An  initial  printout  of  the  default  data  contained  in 
the  program,  and 

2.  A  subsequent  printout  of  the  new  data  input  by  a  user. 
This  data  overlays  the  default  data.   Consequently, 
there  is  sometimes  a  residue  of  default  data  printed 
beyond  the  user  specified  input  data.   This  occurs  when 
the  default  data  contained  in  any  particular  table  ex- 
ceeds that  which  is  specified  by  a  user.   This  excess 
data  in  no  way  affects  the  operations  of  the  computer 
program. 
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Deck  Set  Up 

The  program  deck  should  be  set  up  as  follows  z 

Users  job  card 

JCL  specifying  program  source  and/or  program  source  deck 
or  object  deck  or  combination  of  both 

Data  output  specification  -  JCL  consisting  of  a  card 
specifying  printer  output  on  device  10  if  other  than 
device  03  ,  a  card  specifying  input  (bridge  loading),  on 
device  ITAPE  1  and  a  card  specifying  output  on  device 
number  ITAPE2 . 

Data  input  specification  -  JCL  consisting  of  a  card  spe- 
cifying data  deck  input  IN. 

Data  Deck  set  up  as  described  in  section  on  data  input 
above 

Figure  14,   Deck  Set  Up,  further  illustrates  the  above.  An  ex- 
ample is  found  below  under  the  section  Sample  Case. 

As  a  result  of  limited  testing,  BRGSTRS  appears  to  require, 
for  simple  span  bridges,  from  1/6  to  approximately  1/30  of 
the  real  time,  simulated  in  the  truck  load  data  generation 
by  BRIGLD1,  for  computer  time. 
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Data  Deck 


JCL  -  Data  Deck  Input  Device 
IN  Specification 


JCL  -  Output  Tape  or  Disk  on  Device  #  ITAPE2 


JCL  -  Input  Tape  or  Disk  on  Device  #  ITAPE1 


JCL  -  Printout  Specified  on  Device  10 
(Optional  if  default  device  =  10) 


Program  Source  or  Object  Deck 
(Optional  -  May  be  on  Disk  or  Tape) 


JCL  -  Program  Source  identification 


Job  Card 


Figure  14.   Deck  Set  Up 
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Sample  Case 

A  sample  case  is  provided  here  as  an  example  for  the  user. 
It  is  not  recommended  that  the  source  be  recompiled  for  each 
run.   An  object  program  should  be  stored  on  disk  if  frequent 
use  is  to  be  made  of  the  program.   An  "off-line"  listing  is 
provided  of  the  cards,  JCL ,  and  data,  used  to  run  the  program 
(see  Figure  15) .   An  "on-line"  program  printout  is  provided 
as  Figure  16. 
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//R4111PF5    JOB    (0118, FHRS)f 'PAT    FOTHERG I LL • , Tl ME=20, CLASS=F,REGION=250K 

/*SETUP  R4111PF5,REM0TE8 

/^MESSAGE       INPUT,  007987  *** 

/^MESSAGE       OUTPUT,  007989  *** 

//CU  EXEC  F0RTGCLG,PARM.F0RT=(S0URCE,MAP),TIME=20 
//FORT.SYSIN  DO  * 


//LKED.SYSIN  DD  * 


//GO.SYSUDUMP  DD  DUMMY ,D  ISP= ( NEW  , DELETE ) 
//GO.FT06F001  DD  SYSOUT=A , DCB=( RECFM=UA, BLKSI Z E= 133 ) 
//G0.FT07F001  DD  DSN=FOT HGL ,UN IT=2400-4 , V0L=SER=007987 , 
//       DISP= (OLD, KEEP), LABEL  =  ( 1 , SL  )  , 
//      DCB=(RECFM=VBS,LRECL=1692,8LKSIZE=677  2) 


//GO.FT08F001  DD  DSN=FOT HGL ,UN IT=2400-4 ,V0L=SER=007989 , 
//       DISP=(NEW,KEEP) , LABEL  =  ( 1  ,SL ) , 
//      DCB=(RECFM=VBS,LRECL=1692,BLKSIZE=6  772) 
//GO.FT05F001  DD  * 


Figure  15.   JCL 

34 


1    o 

o 

i      3 

o  « 

at 

o  e 

•    o 

o 

o 

n 

o 

W) 

ii 

o 

3 

i- 

• 

Z 

oc 

r- 

:  »m 

o 

:    o> 

fr 

i     <*• 

ii 

<f 

r+ 

l     o> 

K    • 

m  o» 

►  O 

o>  o- 

o 

O  <M 

if.   • 

o© 

o> 

'o  n 

cm  •'i 

,  •  3  © 

©Z 

1 

» 

o 

III 

SO 

II 

II  o 

3 

_l  o 

Z 

at 

l/> 

►  o 

-i 

-1 

o 

CM 

;    o  o 

► 

o 

» 

'     o 

•  ir\ 

o  CM 

o  -0 

o 

o   • 

o   » 

o  o 

oo 

o 

o 

O   II 

o 

oct 

o 

co  o 

o 

a 

o  _ 
CO  o 


a. 
o 

GO  O 


o  d  a  » 

•  oc  _l  o 

OJU 

O  <M     » 

1 

o  o  o 

•  Ob 

O  in  CM 

O  *-■ 

o          o 

o         o 

O  II        o 

col-      o 

<*  X         Oil 

o      o  o 

►*        O  3 

||   UJ            •  to  O 

o  3:      ->a 

ui   »  »       z 

►      a.        ► 

3    II 

GO  3 

Z  at 

O     ►  _l 

CM  O  CM  O  O 

a  o       »o 

0*  O         CM  O 

•*  O          »H  O      • 

CM  O              •  a 

; 

-<  •         o 

•o  o          o 

•  m             o 

CO                     CO 

1 

o 

II 

ii  x  ii       n 

1 

_)  <^>  00        o 

1 

u  a  z  ii  « 

col    ►  _l  to 

o 

>-  1-        a:  q    < 

z 

CM  <-<        Z 

1 

ai 

o        •■ 

a. 

1    o  0* 

< 

111  o 

i 

i- 

c-  o       o 
&■  O       o 

i- 

a*  o      o 

3 

<7-  o       o 

a 

o-o      o 

H 

o    .         • 

3 

O  .0        o 

k 

a 

II  •*           o  c 
_i  •       -  co 

Z  O    II     II 
II    _i  ct  ii 

> 

in 

1-  CX  _i  to 
o  o  Z  I 

►     »     »|- 

o             o 

, 

o              — 

o 

o                uo 

z 

o                  •• 

o  m  o  o  o  c 

> 

UJ 

O  0-  o  o 

a. 

•  a>  o  o 

< 

Off1  OO 

t- 

<  n0  o  o  o 
1—       o  o  u"» 

I- 

3        a  o  CM 

3 

a.  ii  cm    .  o 

o 

Cl 

ZU       IH       I 

z 

Z 

-JO       o 

UJ 

<J  X 

u 

co 

>-  ■© 

o 


_l  <M 

■-.  O 

I 


ex 

3  cm 

-,  o 

I 


CD 

1 

J© 

>- 

o>     ! 

UJ      • 

a. 

** 

a 

a  © 

o 

o> 

o 

o  o 

• 

oo 

• 

CM       | 

• 

o 

at 

•4        1 
1 

1 

1 

o 

-i 
at 

-1  o 
CM 

Z 

_l  CM 

1-  o 

m 

o 

1- 

o 

o  o 

1     i 

oo 

3 

o 

Ul 

o 

a.  o 

• 

0* 

t 

_i 

•r 

•t 

J- 

CM 

CO 

o     | 

CO 

Z 
CO 

_J 
<  o 

• 

«  o 

M 

o 

o 

o  o 

oo 

•        3  CM 

O        -  © 

CM  I 

Ul 


o  o 
o  o 
-■  o 


oo 

CO 

Q. 

o  o 

-J 

o 

•-.  o 

>-  o 

• 

• 

irv 

© 

CO 

CO 

•o 

1- 

OO 

a.  cm 

o 

z 

z  ^ 

•o 

o 

o 

»■* 

CO 

h- 

t 

3 

<r 

l/l 

oc  o<      o 

o 

z          © 

0. 

o 

t- 

T  O 

cc 

X           <■  -0 

—  a 


3 

_. 

CM 

X  O 

CM 

or  o>      o 

< 

in 

_l  o 

1- 

z          © 

CO 

UJ  O 

a. 

o 

• 

o    • 

O 

-0 

o 

-jo 

(M 

rH 

X            --O 

<  CM  Z 

p* 

UJ         ~ 

1- 

co      o 

0. 

o 

o 

X 

m 

O  CM 

o 

o  r- 

•  z 

m  •- 

-•  o 


O  O         -' 

o  o 

O  O        I- 

•   •      z 

vr\  cm       '-, 

o 


o  o 
o  o 
o  o 


I 
I 

oo 
oo 

i 

I 
I 

i 
oo 

oo 

Oo 

oo 
oo 
oo 

•to 

CMO 

~>o 


-p 
o 
c 

■H 

0: 

c 

•H 

1^1 
I 

C 
O 


$-1 

•H 


35 


c 

Q 

a 
-*1 


t4-t 


»      :   •     ■  — 


a 
r 

•H 

o 
u 


4-1 

o 

-p 

P; 

a; 
G; 
■Hi 

t  ■ 
ci 

ni 


i 

3  I 
-H 


36 


STRESS-EVENT  2  LUAI)    EVENT    493 


BEAM    NO.  2  POINT 

STRESSES:    MINIMUM  TIME 

1.3  3  0.050 

1184.04  0.600 

-2841.16  0.850 

BEAM    IMG.  2  POINT 

STRESSES:    MINIMUM  TIME 

-1290.49  0.250 

-2452.19  0.500 

-29919.02  0.900 

BEAM- -NO-..      ~2-  _  --RX3I-N-T-  -  .. 

STRESSES:    MINIMUM         TIME 

-512.86      0.400 

-3349.30       0.700 


15.0 00 

MAX  I MOM 

TIME 

RANGE 

213  2.35 

0.400 

? 1  31.02 

2230.21 

0.650 

1046. ] 7 

885  .78 

1.05  0 

3 72 6 ,94 

35.000 

MAX  I  MUM 

TIME 

RANGE 

5893.46 

0.400 

718  3.96 

3  0  310.67 

0.7  00 

3  2  762.86 

33  2  7.03 

1.05  0 

33246.05 

4  5.0  00 

MAX  I  MOM 

T  I  M  E 

RANGE 

577  .31 

0.5  00 

]  09 (J.  17 

4129.27 

0.9  00 

747  8.5  7 

figure   16.      On-Line   Printout    (Continued) 
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PROGRAM  DESCRIPTION 

Introduction 

The  Dynamic  Stress  Analysis  Program  is  described  here  in 
terms  of  its  individual  subroutines  and  data  interfaces. 
The  section  on  each  subroutine  will  contain  a  description 
of  the  routine,  the  program  flow  chart,  and  a  program  list- 
ing. Preceeding  the  subroutine  descriptions  is  a  list  of 
program  variables.  This  list  contains  all  variables  that  are 
used  by  more  than  one  subroutine  or  are  used  for  program 
data  input.  Units  or  data  type  associated  with  the  variable 
are  provided  where  applicable,  table  dimensions,  and  a  list 
of  subroutines  using  each  variable  are  included. 

List  of  Variables 

The  program  List  of  Variables  contains  a  list  of  all  pro- 
gram names  and  an  identifying  letter,  a  list  of  common  in- 
dices and  their  size  and  meanings,  and  a  list  of  all  data 
elements  which  provide  interfaces  between  program  subrou- 
tines or  those  data  elements  which  are  required  for  data 
input  or  output.  Consequently,  program  data  interfaces  will 
not  be  restated  in  the  Program  Routines  section. 
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List  of  Variables 

C 
C 

c    PROGRAM  SYMBOLS     MEANING 
C 

c  A     ASSMBL    CALCULATE  INTERNAL  MOMENTS,  BENDING  MOMENTS, 
£  DEFLECTIONS,  STIFFNESS  COEFFICIENTS  AND 

C  ASSEMBLE  STIFFNESS  MATRICES 

C 

c         B     BRGDYN    CALCULATE  BRIDGE  DYNAMICS 
C 

c         D     DMOMNT    CALCULATE  UNIT  LOAD  DISCONTINUITY  MOMENTS 

C   ~~ 


B 

BRGDYN 

D 

DMOMNT 

E 

ERROR 

G 

GRDNET 

I 

..  INPTBD 

J 

INIT 

L 

LODDIS 

M 

MODULI 

N 

INCRMT 

P 

PNTLOD 

R 

GRDLOD 

S 

STIFF 

T 

STIFFL 

U 

STIFFA 

ERROR  MESSAGE  PRINTOUT  ROUTINE 


INPUT  BRIDGE  DATA 


INCRIMENT  TIME  AND  LOAD  POSITIONS 


C 
C 

c         G     GRDNET    ESTABLISH  GRID  NETWORK 
C 

C 

C 

c         J     INIT      INITIALIZE  DATA  BLOCKS 
C 

C         L     LODDIS    CALCULATE  ACTUAL  DISCONTINUITY  MOMENTS 
C 

c         M     MODULI    CALCULATE  MODULI  AND  STRUCTURAL  PARAMETERS 
C 

C 

c         P     PNTLOD    DEVELOP  POINT  LOADS 
C 

c         R     GRDLOD    CALCULATE  GRID  LOADS  /  SPAN 
C 

c         S     STIFF     CALCULATE  LONGITUDINAL  STIFFNESS  COEFFICIENTS 
C 

c         T     STIFFL    CALCULATE  LATERAL  STIFFNESS  COEFFICIENTS 

C 

C         U     STIFFA    ASSEMBLE  THE  STIFFNESS  MATRIX 
C 

C 

C         Q     MAXEXT    EXTRACT  STRESS  PEAKS 

C 

C 

C 

C 

C 

c 

c    INDICES 


X     BREXEC    MAIN  ROUTINE 


c         I    100    GRID  PIONT  INDEX 
C  II         TEMPORARY  INDEX 


u         J  10  BEAM  INDEX 

C         K  20  GRID  PIONTS  (  ONE  SPAN  ) 

C         L  10  SPAN  INDEX 

C         LL  11  SUPPORT  INDEX 

C         M  3  CONCRETE  GIRDER  SECTION  INDEX 

N  3  CROSS  SECTION  INDEX 
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List  of  Variables  (Continued) 

C 

C    LABEL  DIM   TYPE  DESCRIPTION  PROGRAM 

C 

C    AI     N,J        GIRDER  CROSS-SECTIONAL  AREA  IM 

C 

C    ALB    L  AREA  OF  STEEL   REAL  OR  EQUIVALENT  BEAM  IN  THE         IM 

C  LATERAL  DIRECTION 

C 

C 

C 

C 

C 

C 

C    B      J  EFFECTIVE  BEAM  WIDTH  (FEET)  IM 

C 

C    BCP    N,J        WIDTH  OF  COVER  PLATE  IN  THE  N  TH  CROSS  SECTION  OF     IM 

C  THE  J  TH  BEAM 

C 

C    BG     M,N,J      WIDTH  OF  THE  M  TH  SECTION  OF  THE  N  TH  CROSS  SECTION   IM 

C  OF  THE  J  TH  CONCRETE  GIRDER 

C 

C 

C 

C 

C    C      N,J        EFFECTIVE  NEUTRAL  AXIS  DISTANCE  FROM  DECK  SURFACE     IM 

C 

C    CCP    N,J        DISTANCE  TO  THE  NEUTRAL  AXIS  OF  THE  BOTTOM  COVER      IM 

C  PLATE  OF  THE  N  TH  CROSS  SECTION  OF  THE 

C  J  TH  BEAM 

C 

C    CG     M,N,J      DISTANCE  TO  THE  NEUTRAL  AXIS  OF  THE  M  TH  SECTION      I 

C  OF  THE  N  TH  CROSS  SECTION  OF  THE  J  TH 

C  CONCRETE  GIRDER 

C 

C    CI     N,J        DISTANCE  TO  THE  NEUTRAL  AXIS  OF  THE  N  TH  CROSS        IM 

C  SECTION  OF  THE  J  TH  GIRDER 

C 
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List  of  Variables  (Continued) 

C  LABEL  DIM   TYPE  DESCRIPTION                                            PROGRAM 

C 

C  CLRL             DISTANCE  FROM  THE  DECK  SURFACE  TO  THE  NEUTRAL         IM 

C  AXIS  OF  THE  LOWER  LAYER  OF  REINFORCING 

C  STEEL  RODS   IN  THE  LATERAL  CROSS  SECTION 

C  OF  DECK 

C 

C  CLRU             DISTANCE  FROM  THE  DECK  SURFACE  TO  THE  NEUTRAL         IM 

C  AXIS  OF  THE  UPPER  LAYER  OF  REINFORCING 

C  STEEL  RODS   IN  THE  LATERAL  CROSS  SECTION 

C  OF  DECK 

C 

C  CLT              TOTAL  DISTANCE  TO  BOTTOM  OF  LATERAL  BEAM  =T1+T2       IM 

C 

C  CRL    M,N,J      DISTANCE  TO  THE  NEUTRAL  AXIS  OF  THE  LOWER  LAYER  OF    IM 

C  REINFORCING  STEEL  RODS  M  TH  SECTION  OF  THE 

C  N  TH  CROSS  SECTION  OF  THE  J  TH  GIRDER 

C 

C  CRLD             DISTANCE  FROM  UPPER  DECK  SURFACE,  TOP  TO  CENTER  OF    IM 

C  LOWER  LAYER  OF  REINFORCING  STEEL  RODS  OF 

C  THE  DECK 

C 

C  CRU    M,N,J      DISTANCE  TO  THE  NEUTRAL  AXIS  OF  THE  UPPER  LAYER  OF    IM 

C  REINFORCING  STEEL  RODS  M  TH  SECTION  OF  THE 

C  N  TH  CROSS  SECTION  OF  THE  J  TH  GIRDER 

C 

C  CRUD             DISTANCE  FROM  UPPER  DECK  SURFACE,  TOP  TO  CENTER  OF    IM 

C  UPPER  LAYER  OF  REINFORCING  STEEL  RODS  OF 

C  THE  DECK 

C 

C  CT     N,J        TOTAL  DISTANCE  FROM  DECK  SURFACE,  TOP  TO  BOTTOM  OF    IM 

C  THE   N  TH  CROSS  SECTION,  GIRDER  BOTTOM, 

C  OR  BOTTOM  COVER  BOTTOM,  OF  J  TH  BEAM 

C 

C 

C  C2LRL             DISTANCE  FROM  THE  TOP  OF  THE  DECK  TO  THE  NEUTRAL      IM 

C  AXIS  OF  THE  LOWER  LAYER  OF  REINFORCING 

C  STEEL  RODS,  CONCRETE  DIAPHRAM 

C 

C  C2LRU             DISTANCE  FROM  THE  TOP  OF  THE  DECK  TO  THE  NEUTRAL      IM 

C  AXIS  OF  THE  UPPER  LAYER  OF  REINFORCING 

C  STEEL  RODS,  CONCRETE  DIAPHRAM 

C 

C 
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List  of  Variables  (Continued) 

C    LABEL  DIM   TYPE  DESCRIPTION  PROGRAM 

C 

C    DELW   L  GRID  SQUARE  WEIGHT  PER  SPAN  (LB)  GB 

C 

C    DELX   L  GRID  INCREMENTAL  DISTANCE  PER  SPAN  (FT)  GMRA 

C  IN  LONGITUDINAL  DIRECTION 

C 

C    DELXP  L  WIDTH  OF  THE  CONCRETE  DIAPHRAM  IM 

C 

C    DLF    I  DYNAMIC  LOAD  FACTOR  XB 

C 

C    DFLND  I  DYNAMIC  VELOCITIES  XB 

C 

C         DFLNDDI  DYNAMIC    ACCELERATIONS  XB 

C 

C    DFLN   I  DYNAMIC  DEFLECTIONS  XB 

C 

C    DR  DIAMETER  OF  REINFORCING  STEEL  ROD  USED  IN  DECK         IM 

C 

C    DRG    N,J        DIAMETER  OF  THE  REINFORCING  STEEL  ROD   IN  THE   N  TH   IM 

C  CROSS  SECTION  OF  THE  J  TH   CONCRETE  GIRDER 

C 

C    DRT  DIAMETER  OF  REINFORCING  STEEL  RODS  IN  TRANSVERSE      I 

C  DIRECTION 

C 

C    DT  INTEGRATION  INTERVAL  FOR  DYNAMIC  CALCULATIONS  IX 

C 

C    DTL  INTEGRATION  INTERVAL  FOR  LOAD  GENERATOR  X 

C 

C    DXPOSX  II        DISTANCE  OF  TRAVEL  DURING  SIMULATION  INTERVAL  IX 

C 

C 

C 

C 

C    E      N,J        STEEL/CONCRETE  MODULUS,  DEPENDING  ON  COMPOSITION  OF   EIM 

C  J  TH  GIRDER 

C 

C    EBAR   N,J        EFFECTIVE  YOUNGS  MODULUS  IM 

C 

C    EC  YOUNG'S  MODULUS  FOR  CONCRETE  IM 

C 

C    EEI    LtJ        AVERAGED  PRODUCT  OF  THE  EFFECTIVE  BEAM  MODULUS        DA 

C  OF  ELASTICITY  AND  MOMENT  OF  INERTIA 

C 

C    EI     N,J        PRODUCT  OF  THE  J  TH  EFFECTIVE  BEAM  MOMENT  OF  MSDA 

C  INERTIA  AT  CROSS  SECTION  N 

C 

C    EIBL   N,J        EI  FOR  EFFECTIVE  BEAM  -  LONGITUDINAL  M 

C 

C    EICG   N,J        EC  *  ICG  M 

C 

C    EICP   N,J        ES  *  ICP  M 

C 

C    EID    N,J        EC  *  ID  M 

C 

C    EIG    N,J        ES  *  IG  OR  EICG  +  EIRG  M 

c 

C    EIL    L  EI  FOR  LATERAL  SECTION  L  TH  SPAN      43  MS 


List  of  Variables  (Continued) 


EI  LB  L 

ELLBC  L 

EILBR  L 

EILD  L 

ELLDC  L 

EILR  L 

EIR  N,J 

E I  RG  N  »  J 
EPSLON 

ES 

ET  L 

F  I, J 

IB  J 

IBL  N,J 

ICG  N,J 

ICP  N,J 


ES  *  ILB  OR  EILBR  +  EILBC 

EC  *  ILBC 

ES  *  ILBR 

EILDC  +  EILR 

EC  *  ILDC  . 

ES  *  ILR 

ES  *  IR 

ES  *  IRG 

PERCENT  ERROR  ALLOWED  IN  THE  NUMERICAL  SOLUTION 
OF  THE  NON-HOHOGENOUS  DIFFERENTIAL 
EQUATIONS 

YOUNGS  MODULUS  FOR  STEEL 

EFFECTIVE  YOUNGS  MODULUS  PER  DX  WIDTH  IN  EACH  SPAN 


ID 
IE 
IEVENT 


N,J 


IF  =0  READ  ALL  BEAM  PARAMETERS 

IF  NOT  =0  READ  EFFECTIVE  BEAM  PARAMETERS 

CALCULATED  MOMENT  OF  INERTIA  -  LONGITUDINAL  -  OF 
EFFECTIVE  BEAM  WIDTH  =  INRT(N,J) 

MOMENT  OF  INERTIA 

MOMENT  OF  INERTIA  -  COVER  PLATE 

MOMENT  OF  INERTIA  -  DECK 

EVENT  NUMBER  -  BRIDGE 

EVENT  NUMBER  -  LOAD  GENERATOR 


M 
M 
M 
M 
M 
M 
M 
M 
I 

IM 
IM 

RLA 


M 

M 

M 

M 

XO 

WZ 
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List  of  Variables  (Continued) 

C    LABEL  DIM   TYPE  DESCRIPTION  PROGRAM 

C 

C 

C    IG     N,J     R   MOMENT  OF  INERTIA  M 

C    IGO    N,J     R   MOMENT  OF  INERTIA  OF  N  TH  CROSS  SECTION  OF  J  TH       IM 

C  GIRDER 

C 

C    II     N,J     R   TOTAL  MOMENT  OF  INERTIA  M 

C 

C    IL  IF  =0  READ  LATERAL  BEAM  PARAMETERS  I 

C  IF  NOT  =0  READ  EFFECTIVE  BEAM  PARAMETERS 

C 

C    ILB    L       R   MOMENT  OF  INERTIA  OF  STEEL  REAL  OR  EQUIVALENT  IM 

C  BEAM 

C 

C  ILBC       L  R       MOMENT    OF     INERTIA  M 

C 

C    ILBR   L       R   MOMENT  OF  INERTIA  M 

C 

C    ILDC   L  .      R   MOMENT  OF  INERTIA  M 

C 

C    ILR    J       R   TOTAL  MOMENT  OF  INERTIA  OF  A  STEEL  ROD  IN  THE  IM 

C  LOWER  PART  OF  THE  DIAPHRAM 

C 

C    ILSC  STEEL  OR  CONCRETE  DIAPHRAM  INDICATOR  IM 

C  0  =  STEEL,  NOT  0  =  CONCRETE 

C 

C  IN  INPUT    DEVICE    NUMBER  I 

C 

C    INRT   N,J     R   EFFECTIVE  MOMENT  OF  INERTIA  FOR  THE  J  TH  EFFECTIVE    IM 

C  BEAM 

C 

C         10  OUTPUT    DEVICE    NUMBER  IMGPR 

C 

C    IR     N,J     R   MOMENT  OF  INERTIA  M 

C 

C    IRG    N,J     R   MOMENT  OF  INERTIA  M 

C 

C    IRGO   M,N,J   R   TOTAL  MOMENT  OF  INERTIA  OF  THE  STEEL  REINFORCING      IM 

C  RODS  OF  THE  M  TH  SECTION  OF  THE  N  TH  CROSS 

C  SECTION  OF  THE  J  TH  GIRDER 

C 

C    IS     L  NUMBER  OF  GRIDS  LONGITUDINALLY  ON  A  GIVEN  SPAN        GDESRLA 

C 

C    ISC    J  STEEL  OR  CONCRETE  GIRDER  INDEX  IM 

C  NOT  0  =  STEEL,   0  =  CONCRETE 

C 

C    ISS  STEADY  STATE  INDICATOR  0=  STEADY  STATE, 

C  NOT  0  =  TRANSIENT  STATE 
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List  of  Variables  (Continued) 

C    LABEL  DIM   TYPE  DESCRIPTION  PROGRAM 

C 

C    IT     L       R   EFFECTIVE  DIAPHRAM  MOMENT  OF  INERTIA  PER  DX  WIDTH     IM 

C  IN  EACH  L  SPAN 

C 

C    ITIME  OUTPUT  TIME  POINT  INDEX  XOO 

C 

C    IUR    J       R   MOMENT  OF  INERTIA  OF  A  STEEL  ROD  IN  THE  UPPER         IM 

C  PART  OF  A  LATERAL  BEAM 

C 

C    KTIME  OUTPUT  TIME  POINT  INDEX  XOO 

C 

C 

C 

C    MB     I  BENDING  MOMENTS  AR 

C 

C    MD     LL,K,J  R   UNIT  DISCONTINUITY  MOMENTS  AT  SUPPORTS  LL,  FOR        DL 

C  EACH  BEAM  J 

C 

C    MDT    LLtJ    R   DISCONTINUITY  MOMENTS  AT  SUPPORTS  LL,  FOR  EACH        LA 

C  BEAM  J 

C 

C    MI     I  INTERIOR  MOMENTS  A 

C 

C 

C    NA  NUMBER  OF  AXLES  FROM  DATA  BLOCK  PR 

C 

C    NAP  NUMBER  OF  LONGITUDINAL  POINT  LOADS  PGR 

C  (USUALLY  EQUAL  TO  NA) 

C 

C    NB  NUMBER  OF  MAIN  LONGITUDINAL  BEAMS  IMASRLG 

C 

C    NCS    J  NUMBER  OF  VARIATIONS  IN  THE  J  TH  MAIN  GIRDER  CROSS    IMS 

C  SECTION 

C 

C    NDBUG  DEBUG  PRINTOUT  OPTIONS 

C 

C    NGRL   M,N,J      NUMBER  OF  REINFORCING  STEEL  RODS  IN  THE  LOWER  LAYER   IM 

C  OF  THE  M  TH  SECTION  OF  THE  N  TH  CROSS 

C  SECTION  OF  THE  J  TH  CONCRETE  GIRDER 

C 

C    NGRU   M,N,J      NUMBER  OF  REINFORCING  STEEL  RODS  IN  THE  UPPER  LAYER   IM 

C  OF  THE  M  TH  SECTION  OF  THE  N  TH  CROSS 

C  SECTION  OF  THE  J  TH  CONCRETE  GIRDER 

C 

C    NI  DTLOAD  /  DT  XY 

C 

C    NL  NUMBER  OF  LANES  I 

C 

C    NLRL  NUMBER  OF  REINFORCING  STEEL  RODS  IN  THE  LOWER  LAYER   IM 

C  OF  THE  LATERAL  CROSS  SECTION  OF  DECK 

C 

C    NLRU  NUMBER  OF  REINFORCING  STEEL  RODS  IN  THE  UPPER  LAYER   IM 

C  OF  THE  LATERAL  CROSS  SECTION  OF  DECK 
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List  of  Variables  (Continued) 


c 

LABEL 

DIM   TY 

c 

c 

NPTS 

J 

c 

c 

NRL 

J 

c 

c 

c 

NRU 

J 

c 

c 

NSI 

c 

c 

NSUP 

c 

c 

NUC 

c 

c 

NUS 

c 

c 

NXPTS 

c 

c 

N2LRL 

L 

c 

c 

c 

N2i.RU 

L 

c 

c 

c 

PTS 

JJ,J 

c 

c 

SIG 

I 

c 

c 

SIGB 

I, J 

c 

c 

SIGP 

JJ,J,KK 

c 

c 

SK 

II, II, J 

c 

c 

SKP 

J, J 

c 

c 

STEADY 

c 

c 

c 

T 

c 

c 

TCP 

N,J 

c 

c 

c 

TG 

M,N,J 

c 

c 

c 

THRSH 

c 

c 

Tl 

c 

c 

T2 

L 

PE  DESCRIPTION!  PROGRAM 

NUMBER  OF  SAMPLE  POINTS  PER  BEAM  IGBO 

NUMBER  OF  LOWER  REINFORCING  STEEL  RODS  IN  EACH        IM 
BEAM  SECTION 

NUMBER  OF  UPPER  REINFORCING  STEEL  RODS  IN  EACH        IM 

BEAM  SECTION 
NSUP  -  1  IXYGMDASRL 

NUMBER  OF  MAIN  BEAM  SUPPORTS  2  <=  NSUP  <=  11        IXYL 
R   POISSON'S  RATIO  FOR  CONCRETE  I 

R   POISSON'S  RATIO  FOR  STEEL  I 

NUMBER  OF  X  GRID  POINTS  GBOD 

NUMBER  OF  REINFORCING  STEEL  RODS  IN  THE  LOWER  LAYER   IM 
OF  THE  REAL  OR  EQUIVALENT  CONCRETE  DIAPHRAGM 

NUMBER  OF  REINFORCING  STEEL  RODS  IN  THE  UPPER  LAYER   IM 
OF  THE  REAL  OR  EQUIVALENT  CONCRETE  DIAPHRAGM 


SAMPLE  POINTE  FOR  STRESS  OUTPUT 

DYNAMIC  BENDING  STRESSES 

STATIC  BENDING  STRESSES 

DYNAMIC  BENDING  STRESSES  AT  SAMPLE  POINTS 


STEADY  STATE  INDICATOR 

TIME 

THICKNESS  OF  EACH  PLATE  IN  THE  N  TH  CROSS  SECTION 
OF  THE   J  TH  BEAM 

THICKNESS  OF  THE  M  TH  SECTION  OF  THE  N  TH  CROSS 
SECTION  OF  THE   J  TH  CONCRETE  GIRDER 

THRESHOLD  VALUE  TO  DETERMINE  STEADY  STATE 

MEAN  THICKNESS  OF  BRIDGE  DECK  (FEET) 

THICKNESS  OF  A  CONCRETE  DIAPHRAGM 
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IGBO 

XB 

AB 

XO 

SA 

SA 

XB 

IM 
IM 

XIB 

IM 
IM 


List  of  Variables  (Continued) 

C    LABEL  DIM   TYPE  DESCRIPTION  PROGRAM 

C 

C    WP     IP, J       WEIGHT  PER  LOAD  POINT   1.=IT.=4  PR 

C 

C    WEIGHT  WEIGHT  OF  THE  BRIDGE  (LBS)  IG 

C 

C 

C    X      I  GRID  CENTER  POSITIONS   (FT)  GDSRA 

C 

C    XCS    N,J        LONGITUDINAL  DISTANCE  TO  THE  N  TH  DISCONTINUITY  IN    ISD 

C  THE  J  TH  GIRDER  (FEET) 

C 

C    XIL    L  MOMENT  OF  INERTIA  -  LATERAL  SECTION  PER  L  TH  SPAN     M 

C 

C    XILB   J,L        AVERAGED  MOMENT  OF  INERTIA  OF  EFFECTIVE  BEAM  A 

C 

C    XILR   L  MOMENT  OF  INERTIA  OF  THE  RODS  IN  A  LATERAL  SECTION^ 

C  OF  THE  L  TH  SPAN 

C 

C    XLB  LENGTH  OF  BRIDGE  (FEET)  IGR 

C 

C    XP     II         POSITION  OF  AXLE  LOADS  ZPR 

C 

C    XSUP   L  DISTANCE  TO  THE  CENTER  OF  THE  L  TH  SUPPORT  (FEET)     IXGSDA 

C 

C 

C    YB     J  DISTANCE  TO  CENTER  OF  MAIN  BEAM,  GIRDER  IR 

C 

C    YBCL  DISTANCE  TO  CENTER  LINE  OF  BRIDGE  (FEET)  ISR 

C 

C    YCL    J  DISTANCE  TO  CENTER  LINE  OF  THE  J  TH  LANE(FEET)         IP 

C 

C 

C    YP     JT         TRACK  POSITION  PR 

C 

C    ZBL    N,J        CROSS  SECTIONAL  MODULUS  -  EFFECTIVE  LONGITUDINAL      M 

C  BEAM  -  LOWER  HALF 

C 

C    ZBU    N,J        CROSS  SECTIONAL  MODULUS  -  EFFECTIVE  LONGITUDINAL      MD 

C  BEAM  -  UPPER  HALF 

C 

C    ZCL    L,J        AVERAGED  CROSS  SECTIONAL  MODULUS  -  LOWER  DA 

C 

C    ZCU    L,J        AVERAGED  CROSS  SECTIONAL  MODULUS  -  UPPER  DA 

C 

C    ZL     N,J  M 

C 

C 

C 

C 

C 

C 
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Program  Routines 

The  following  describes  the  Dynamic  Stress  Analysis  Pro- 
gram in  terms  of  its  individual  subroutines.  The  program 
functions,  equations  solved,  program  interfaces  are  de- 
scribed, followed  by  the  program  flow  chart  and  program 
listing.  Figure  17,  Program  Interfaces,  depicts  the 
hirearchy  of  subroutine  calls. 
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Initialize 


BREXEC 


Loop 


INIT 


INPTBD 


GRDNET 


MODULI 


STIFF 


JConti 


nuous 


Beam 


DMOMNT 


STIFFL 


STIFFA 


MATINV 


PNTLOD 

LODDIS 

GRDLOD 

Cont 

.inuous 

Beam  Only 

ASEMBL 

BRGDYN 

MATMUL 

►BROUT 

MAXEXT 

INCRMT 

STIFFL 


STIFFA 


MATINV 


MATMUL 


MATINV 


n-i 


MATMUL 


Figure   17.   Program  Interfaces 
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BREXEC 

The  main  program,  BREXEC  reads  in  program  input  output  device 
selection  options  and  output  tape  header.   It  calls  the  sub- 
routine INIT  to  clear  data  buffers  and  reads  all  input  data 
through  subroutine  INPTBD.   It  performs  all  program  initializa- 
tion through  calls  to  GRDNET,  MODULI,  STIFF,  and  DMOMNT.   The 
structural  event  number  is  initialized. 

The  main  program  loop  is  then  entered  and  proceeds  as  follows: 

1.  Read  bridge  load  data  block. 

2.  If  it  is  the  first  pass,  a  number  of  intervals  and 
an  integration  interval  is  calculated. 

3.  Point  loads  are  calculated  P^TLOD. 

4.  Loads  at  grid  points  are  calculated  GRDLOD. 

5.  If  a  continuous  beam  calculate  actual  load  dis- 
continuity moments-LODDIS ,  interior  moments  and  deflections, 
bending  moments  and  stresses  and  assemble  stiffness 
coefficient  -  ASEMBL  . 

6.  Calculate  bridge  dynamics  aril  dynamic  stress  and  load 
factor  -  BRGDYN  . 

7.  If  a  steady  state  has  been  ceachad  the  event  stresses 
are  output  -  BROUT  and  the  event  number  is  updated. 

8.  Load  positions  are  incremented  INCRMT  and  the  loop 
proceeds  from  3  above  or  a  new  data  block  is  read  -  1  above. 


51 


The  program  terminates  when  the  bridge  load  data  is  exhausted, 
The  following  subroutines  are  called  by  this  program: 


INIT 

INPTBD 

GRDNET 

MODULI 

STIFF 

DMOMNT 

PNTLOD 

GRDLOD 

BRGDYN 

BROUT 

INCRMT 
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f     BREXEC    ) 


INPUT 

PROGRAM 

OPTIONS 


WRITE  OUT 

PROGRAM 

OPTIONS 


INIT 


INITIALIZE 
DATA  BLOCKS 


TMPTRn 


INPUT  RUN 
PARAMETERS 


"7 


ITAPE1 
ITAPE2 
HEADER 


A 


"7 


HEADER 
ITAPE1 
ITAPE2 


GRHNET 


w 


Figure   18.    BREXEC  Program  Flow  Chart 
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EVENT  NUM 

NUMBER 

IE=1 


Figure    18.   BREXEC  Program  Flow  Chart  (Continued) 
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READ  BRIDGE 
LOAD  DATA 
BLOCK 


Yes 


Figure   18.   BREXEC  Program  Flow  Chart  (Continued) 
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Yes  ^X 

OUTPUT  ALL  STRESSES 
DEFLECTIONS  AND 
DYNAMIC  LOAD 
FACTORS 

\ 

0,6     , 

XX 

DLF  ^ 

/ 

Figure   18.   BREXEC  Program  Flow  Chart  (Continued) 
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Yes    >^LAST 
JLOAD? 


UPDATE  EVENT 

NUMBER 

IE=IE+1 


OUTPUT  STRESS 
HISTOGRAM 


B 


;  HISTOGRAMS 


(^   EXIT    J 


Figure   18.    BREXEC  Program  Flow  Chart  (Continued) 
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ORTRAN  IV  G  LEVEL   21  MAIN  DATE  =  73172  20/34/5" 


C 
__C.   BR EXEC 
C 

OJ)_OL REAL  NUC,NUS,IUR,ILR,IGO,IRGO,INRT,  ILB,  IT 

0002  REAL  MD,  MOT 

0003 I  nGTf.AI   I  ASTP 

0004  LOGICAL  STEADY 

0005 _JNTEGER  SIGC 

C 
000ft COMMON/ BLK1/CLRL,CLRU,  CLT,  CRLD,  CRUD,  C2LRL,  C2LRU,  DRt  DRT, 

1    EC,  ES,  1LSC,  NLRL,  NLRU,  Tl ,  IL,  NUC,  NUS, 
2    NRU(IO),  NRL(10)t  B(10),  I_SC  (  10 )  ,_  A  LB  ( 10  )  ,    N2  L  RU (  1 OJ  ,  _ 

3    N2LRLQ0),  T2(10),  IT(10),  IUR(IO),  ILR (  10 ) , I LB[ 10 ) , 
4_   D_ELXP(10),  ETI10),_ 

5    AI<3,10>,  CIO, 10),  CT(3,10),  CCP(3,10),  TCP(3,10),  BCP(3,10), 
6    DRG(3,10),  C(3,10),  IG0(3,10),  INRT(3,10),  EBAR(3,10), 

7    BG(3,3,10),  CG(3,3,10),  TG(3,3,10),  NGRU( 3, 3, 10)  ,  NGRL (3 ,3  ,10 ) , 
a    r.Rtl(3f3f10)T  r.RI  (3f  3.10),  IRG0(3.3,10) ....._  _ _ 

9    ,DUM(1568) 


0007  COMMON  /BLK2/  DT ,       THRSH,   IE,  IEVENT,  IN,  10,  ISS,  ITER, 

. 1 NAPj.  NB,  NO,  NI,  NL,  NSI,  NSUP,  STEADY,  T,  WEIGHT,  XLB,  YBCL, 

2    NXPTS,  NDBUG(IO), 

_3 DELW(IO),  DELX(  10),  EE.IJ10j  10  )_j_  EIL  (  10)  ,  IS<10),  NCS(IO), 


4    XSUP(IO),  YB(10),  YP(4),  YCL(IO),  ZCL(10,10),  ZCU(10,10), 

6  __E(3,10),  EI(3,10),  XCS(3,10),  ZBL(3,10),  ZBU<3,10), 

7  f"(80,10),  MD(  11,20,  10),  MDT(  11,10),  WP(50,4),  X(80),XP(50) 

C  _ 

0008  COMMON  /BLK3/  SUMHRX,  DTL,  IEVNTX,  NOAXLX,  NTRUKX,  LAST, 

1      LTYPE(20),  WGTX(20),SPDX(20),  LLANE(20),  TIMEX(20), 


2      XP0SX(50),  LNUMX150),  WEITX(50),  DXPOSX( 50 ) , ACCL RX ( 50 ) 


0009  COMMON  /BLK5/  S IG ( 80 ) , DFLN ( 80 ) , DFL ( 80 , 10 ) , I T IM E, KT IME , 

1  SIGP(120,3,10),  SIGC(10,3,10),  SIGTHS,  DSIG 

2  ,  I0PT1,  I0PT2,  I0PT3,  I0PT4,  NPTS(IO),  PTS(3,10) 
1_ -^_LPJ"S(_3_,101._ _ 

C 

0010  COMMON    /BLK6/     IBL(3,10) 
C 

0011  COMMON    /BLK7/WGT,VEL,DELT,BLEN,NEV,NNL,NBEAM,NP( 10) ,DMM(3) ,ITAPE2 
1,SIG0UT(3,10) 


0012  DIMENSION    HISHED(20) 

0013  DIMENSION    HEAD(17) 

0014  DIMENSION    STRESI(IO) 
Q0.15. DIMENSION    L0ADBL(376) 

0016  EQUIVALENCE     ( HI  SHED ( 1 ) ,WGT ) 

QQU £flUIVAL£NCElLDADBLLUi5UMHKX)    .  _       . 

0018  NAMELIST/HISNAM/    WGT, VE L , DELT, BLEN , NE V, NNL , NB EAM , MP 

QQ19 __    DATA     IFTST    /O/ 

C 
._..._.  ..     CALL     IN  IT 

10    =    6 

W^  =  _5_ - 
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-0023- 


0024 

0025 
-00-2.6.- 


0032 
0033 

UQ3A. 

0035 
0036 


QQ3JL 


0038 
G02L9_ 

-0040 


C 

X 

C 

£_.. 

C 


..-INITIALIZATION 

JRPUT    PROGRAM  OPTIONS 

REAR     (  IN, 1001  )    -JJAP-ElOXAPEZ^hlEAD 


1001    FORMAT    ( 


212,     17A4) 


WRITE     (10,1002)     HEAD,     ITAPE1,     ITAPE2 
inn?    FORMAT     aHl.i.JJ7A4,//i'     INPUT    TAPE    NO.',     14, 
1  •  OUTPUT    TAPE    NO. ' , 14) 

C 


0027 . KTIME  _=.  1 

0028  ITIME    =    1 

nn?9 _M£it  j=__0 _ 

0030                                        THRSH    =    300. 
0031 SIGTHS  -=-  300.0 


DSIG    =    3000.0 

ifveni=  a 


£AJJ_  INPTBD 
CALL  GRDNET 
CALL    MODULI 


__J_F_LNSUPJ,EQ.2J 

CALL  DMOMNT 
_.GD_JQ  60 


GO  TO  50 


50  C.AI  I  ST-LEF_LU- 


0041 


0042 


0043 
0044 
0045 
0046 


0047 
0048 
0049 
00-50 
0051 


0052 
0051 
0054 
005  5 
0056 


C 


6Q  IE 


LASTP 


.FALSE, 


PROGRAM  MAIN  LINE 

100  CONTINUE 
130  CONTINUE 

IF(  ITAPE1.EQ.0) 
IE-  I  HARE  _L*£  1 .  5.  J 


GO  TO  140 
GO  JQ_  135 


READ  ( ITAPE1,1120,END=165)  LOADBL 
.1120  FORMAT  (20A4) 

GO  TO  137 
135  READ!  I  TAPE1 , END= 165 )   LOADBL 
137  CONTINUE 


— -G- 


IF  (NDBUGd  J.EQ.O)   GO  TO  150 
__  -K  =  NOAXLX 

WRITE(6,1145)  SUMHRX 
DO  146  1=1, K 
146  WRITE(6,1146)  I,LNUMX(I),  WEI TX ( I  )  , XPOSX(  I  ) , DXPOS X(  I )  ,  ACCLRX ( I  ) 
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0057  1145  FORMAT! •OBRIDGE  LOAD  TIME  = ' , Fl 2. 4, • SEC • // 10X, ' LANE ' ,4X , • WE IGHT 

1j  'POSITION' ,2X, 'D I  STANCE ',2X, 'ACCELERATION   ORDER'/) 

0058  1146  F0RMAT<7X,I2,3X,I2,F10.0,F10.2,F10.2,F14.2) 
O055 WJUIE- .16,1050)  IEVNTX,  NTRUKX  . 

0060  1050  FORMAT  COEVENT  NUMBER ', 14 ,8X ,' NO.  OF  TRUCKS ', 14// , 

1 IX,  'TYPE WEIGHT SPFFP LANE LIME_ENTER1NG  BRIDGED 

0061  DO  151  1=1, NTRUKX 

J1D62- 131  WRITE  (6,1051)  LTYPE (  I  ),WGTX(I  ),SPDX(I  ),LLANE(  I),TIMEX(I) 

0063  1051  FORMAT( IX , 14, F 10.0,F8. 0, I8,F8.0) 

M)±A ISO  CONTINUE 

0065  IF  ( IEVENT.NE.IEVNTX  )   STEADY  =  .TRUE. 

0066 IEVENT  =  IEVNTX 

0067  140  CONTINUE 

0068 I F_U £T. ST  .  GT  .  0 ) .  GO .  TO  _  200 

0069  IFTST  =  IFTST  +  1 

0070 IF  (DT.NE.O)    GO  JO  90 

C 

0071 - PI  ST  =  XSUP(7)  -  XSUP(l) . _ 

0072  IF (NSUP.EQ.2)   GO  TO  80 

0073 D0_75  I=3,JiSUP 

0074  DIS  =  XSUP(I)  -  XSUP(I-l) 

0075 ; I_F_  (DIS.LT.DIST)   DIST  =  DIS 

0076  75  CONTINUE 

0077 80  DT  =  DIST/1200.0 

0078  90  CONTINUE 


0079  XNI  =  DTL/DT  +  .5 

£080 _  N I  =  XN I 

0081  DT    =    DTL/NI 

0082  IF(NDBUG(l).NE.O)    _  W  R  LI  EliQj  1 200  L   D  T ,  D  T  L  t  N I 

NI',2F10.3,I6) 

GO    TO    200 


0083 

1200 
C 

FORMAT< '0  DT   DTI 

0084 
0085 
0086 
0087 

IF  ( ITAPE2.E0.0) 
DELT  =  DT 
BLEN  =  XLB 
NBEAM  =  NB 

0088 
0089 

199 

DO  199  1=1,10 
NP( I)  =  NPTS( I ) 

0090  200    CONTINUE                                                                                                                                                           C 

C  ( 

0_Q9_I LF.LLT APE 2_._E_Q>  0 )       G0_  JO  .  20 1 

0092  IF     ( IEVNTX. EO.NEV)       GO    TO    202 

0093  NEV    =     IEVNTX 

0094  WGT    =    WEITX( 1) 

0095 VEL  =  SPDX(l) 

0096  NNL  =  LNUMX( 1) 


C    CALCULATE  NUMBER  OF  TIME  INCREMENTS 
.  __ X 

0097  NJ  =  (XLB-1.0) /(DT*VEL ) 
C ._ 

0098  WRITE     (ITAPE2)     HISHED 
HCL9S , WRIIE-  LLO ..tLLSNAM J. 
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0100  GO  TO  201 

010 1 202  „N J_ ..=  NJ  -  NI 

0102  IF(NI.EQ.O)   GO  TO  100 

0.1.01.  ZO.l  CONTINUE  .... 

C 
0104 M  =  1 ._ .  

C  i 

01.05 -.       210-CONT INUE 

0106  CALL    PNTLOD  ( 

QJLQ1 IE__l-NkP-.EQ..Q,  AND.SJEADY  J       GO    TO    270 

0108  CALL    GROLOD  i 


0109  IF(NSUP.GT.2)       CALL    LODDIS 


0110  1X4    =    0 

mil DO    230    L  =  1,NSI     _  ' 

0112  1X3    =    1X4    +    1 

Oil  3 1X4    =    TX4    +_IS_LL1 .._..       

0114                                       IF     (NSUP.EQ.2)    GO    TO  250                                                                                                                < 

0115 CALL    ASEMBL(L)      ' 

0116                            250    CONTINUE  < 

JJlli : XA  LL    B  RGDYN  (  L  )  _  ' 

C  < 

0118  IF  (  1 0  P  T  4  .  E  oToT"  ~  GO  To  ~26  8  ~ 
C__     __ 

C    OUTPUT  ALL  STRESSES,  DEFLECTIONS,  AND  DYNAMIC  LOAD  FACTORS 

C    CONVERT  DEFLECTIONS  TO  INCHES 
c_ , _..  .  . 

0119  DO    261     I=1,NXP~TS 

0120     261  DFLNU)  =  DFLN(I)  *  12.0 

C 

0121  _  TIME2  =  DT  *  (ITIME  +  KTIME  -  1) 

0122  WRITE  (10,1261)   IE,  IEVENT,  TIME2 

012  3 1261  FHRMAX-.llML.t_i.S  IR.E.SS-  EVEN  T_»  ,  1 4_,  * .    LOAD  EVENT  ',14,'     TIME  =', 

1       F8.3! 
C 

0124  1X2  =  0 
C 

0125  DO  267  J=1,NB 

_0126 T  XI  -  TX7  +  1 

0127  1X2  =  1X2  +  IS(L) 

0128  WRITE  (10,1262)   L,J 

0129  1262  FORMAT  ( 1  HO, 12X , ' SP AN  NO.  ',13,'     BEAM  NO. • , 13/ , IX » • POS IT  ION '  ) 
.  .  _        C 

0130  WRITE  (10,1263)   (X(I),  1=1X3,1X4) 

01  3J . 1  ?ft3..EHRMAL-LLQF-LZ...3-l _  .    ..„_. 

0132  WRITE  (10,1264) 

0133 1264  FORMAT  ( 1 X , /, 1 X , 'DYN AM  I C  BENDING  STRESSES') 

0134  WRITE  (10,1263)  (  SIG (  I  ) , I =1 X  1 , I X2 ) 
C 

0135  WRITE  (10,1266) 

-013.6 1  2  66  FORMAI— <  J..QD Y.N.AM  IC. DEFLECT.  IONS'). 
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0137 
0L38  . 

Q13SL„ 

0140 
0141 


0142 

0143 
0_14_4_ 

fM45 


0146 
014X-.- 
0148 
_Q_L49__ 
0150 


0151 


0152 
0_153_ 
0154 
JLL5_5_ 


0156^ 

015  7 

0158_ 

0159 

0160 

0161 


WRITE  (  10,1267)  (  DF  HI, L), 1  =  1X1, 1X2) 

1267  FORMAT  (10F12.8) 

WRITE  (10,1268) 

1268  FORMAT(IHO) 

267    r.niMTTNIIF 


230    CONTINUE 

268    CONTINUE 
IE    LJTAPE2.GT.0)       WRI TE ( I TAPE2 )     S I  GOUT 

TF     (KTIME.LT.120) GjQ  _IQ_2_69  


CALL    BROUT 

II IME=1T1ME+120 

KTIME    =    0 
_269    KTIME    =    KTIME    +    1 

IF     (ITAPE2.NE.0)       GO    TO    270 


IF  (STEADY. AND. KTIME. GT. 3)   GO  TO  280 

1 70  CONTINUE 

IF     (NDBUGd  J.EQ.O)       GO    TO    271 
WRITE     (  10,1271  )    M,    NI,    T,     IE 
1271    FORMAT    (  ' 01 NCREMENT    M',14,'  _  NI • , 1 4, • 

1  '  STRESS    EVENT*, 16) 

_271    CONTINUE 

CALL     INCRMT 

M    =    M    +    1 

IF     < ITAPE2.NE.0.AND.M.GT.NJ)       GO    TO    280 

IF     (M.GT.NI  )    GO    TO    290 


TIME'  ,F12.3, 


GO    TO    210 


0162 

0163 

0164 

QJA5_ 

0166 

0167 

0168 

0169 

0170 


0171 
0172 


£LL23_ 


0174 
017S. 
0176 


280    CONTINUE 

CALL    BROUT 

IE    =     IE    +    1 
KTIME    =     l__  _ ._    . 

ITIME    =     1 

IF     (  ITAPE2.NE.0.  )       GO    TO     100 

GO    TO    270 
290    CONTINUE 

GO    TO    100 


165    CONTINUE 

IF     ( I0PT3.EQ.0)        GO    TO    300 

OUTPUT  STRESS  HISTOGRAMS 

D_D_L6A  1=X,  LQ .._     

164    STRESI ( I )     =    DSIG    *    I 

WRITE     (10,1165)     (  STRESI  (I  ), 1  =  1, 10) 
1165    FORMAT     CI  STRESS    HISTOGRAMS',//,'     10%     INTERVALS ',  1.0  (  ' 

1     ,//,'     BEAM    POINT') 


TO'  ,F8.1 ) 


0177 


D£L3_67  ...J  =  1,NB. 
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0178  NPT    =    NPTS(J) 

0179 IF(NPT.EQ.O)       GO    TO    367 

C 

0L30 DD.-L6l2.X  =  1,NPJ 

0131  167    WRITE     (10,1167)     J,PTS(K,J>,     ( SIGC ( I , K, J ) ,  I  =  1 ,  10 ) 

0182 1167    FORMAT    ( 1  X  ,  14  ,    F  LCU3 ,_  1QJ  5X  ,  Ii> )  1 . 

0183  367  CONTINUE 

0184  168  CONTINUE 

JH8.5 300  -STQP-.QJ 

C 

01  86 EM! 
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INIT 

The  subroutine  INIT  initializes  common  blocks  BLK1  and  BLK2 
to  zero. 
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C     INIT     ) 


Zero 
Common 
BLK  1 


Zero 
Common 
BLK  2 


f    RETURN    J 


Figure    19.   INIT  Program  Flow  Chart 
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ORTRAN    IV    G    LEVE4 21 - -     -   IN  IT 

0001  —  --SUBROUTINE    LNIT 

C 
G-      SET   DATA  SLOCKS   TO    ZERO 


0002- 


COMMON    /BLK4^-  -10M  U  640a) 


0003 
00O4 
0005 
0006— 


0007 

oooa- 


C0MM0N  /BLK2/  IDM2(4005> 
COMMON  /BLK4/  -IDM4<  19520) 

COMMON  /BLK5/  IDM5C4938) 
COMMON-  /BLK7/  IDM7(51) 

DO  10  1=1,1208 
10  IQMK-I4— a-0 


0009 DO  20  1=1*4005 

0010  20  IDM2U)  =  0 

0011  DO  40  1=1,19520 
00t2 40  IDM4(I)-=-0 


0013 

0014 

0015 
0016- 


0017- 


0018 


0^40  1=1,4938 

30  IDM5U)  =  0 


DO  70  1=1,51 
J0-IDM7(I)  -=_ 0 - 


RETURN 
END 


DATE  =  73171 


21/01/13 
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INPTBD 

The  INPTBD  reads  all  input  bridge  data  as: 

Single  element  data 
Longitudinal  beam  data 
Lateral  beam  data 
Output  specifications 

All  elements  are  printed  out  after  input.   All  data  input  in 
inches  is  converted  to  feet  for  computation. 

A  detailed  list  of  input  data  is  found  in  the  section  Data 
Input  above. 

No  subroutines  are  called  by  this  routine. 
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( 

INPTBD    J 

Initial 

ize  Default 

Parameters 

NDBUG 

=  0 

IN 

=  5 

10 

=  6 

LAST 

=  TRUE 

NL 

=  2 

ND 

=  1 

READ  SINGLE 
ITEM  INPUT 


ZERO  IT 
ARRAY 


READ  IN 
LANE  CENTER 
LINES  YCL 


C 


V 


WRITE  OUT 

LANE 

CENTER  LINES 


NAMELIST 
INPUTA 


YCL 


A 


YCL 


Figure   20. 


INPTBD  Program  Flow  Chart 
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A 


J  =  1 


-* 


READ  BEAM 
PARAMETERS 


V 


k. 


PRINT  OUT 

BEAM 

PARAMETERS 


±. 


READ  DISTANCES 

TO 
DISCONTINUITIES 
THIS  BEAM 


il 


PRINT  OUT 
DISTANCES 

TO  DISCON- 
TINUITIES 


_\k 


NRU,  NRL,  IUR, 
YB,  B,  NCS,  ISC, 
IB 


NRU,  NRL, 
IUR,  ILR,  YB, 
YB,  B,  NCS, 
ISC,  IB 


XCS 


Figure   20.    INPTBD  Program  Flow  Chart  (Continued) 
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INITIALIZE 
NEUTRAL  AXES 
ARRAYS  C,  CI 


REAL  ALL 

BEAM 

PARAMETERS 


WRITE  ALL 

BEAM 

PARAMETERS 


±L 


Yes 


Concrete 


No 


\ 


AI,IGO,CI, 

CT,CCP,TCP, 

BCP 


AI,IGO,CI,CT, 
CCP,TCP,BCP 


Figure   20.    INPTBD  Program  Flow  Chart  (Continued) 


READ  EFFECTIVE 
STEEL  BEAM 
PARAMETERS 


J!L 


PRINT  OUT 
EFFECTIVE  BEAM 
PARAMETERS 


INRT,  EBAR 
C.  E 


INRT,  EBAR, 
C,  E 


READ  CONCRETE 
BEAM  SECTION 
PARAMETERS 


Z_ 


\ 


WRITE  OUT 
CONCRETE  BEAM 
SECTION 
PARAMETERS 


± 


READ  DIAMETER 
OF  REINFORCING} 
STEEL  RODS 


PRINT  DIAMETER 
OF  REINFORCING 
STEEL  RODS 


S    BG,CG,TG, 

NGRU , NGRL , CRU , 
CRL,IRGO 


BG,CG,TG, 
NGRU, NGRL, 
CRU,CRL, 
IRGO 


DRG 


/ 


DRG 


Figure   20. 


INPTBD  Program  Flow  Chart  (Continued) 
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READ 

SUPPORT 

POSITIONS 


Z_ 


V 


XSUP 


-V 


PRINT 

SUPPORT 

POSITIONS 


^ 


XSUP 


READ  LATERAL 

BEAM 

PARAMETERS 


L 


V 


-*_ 


WRITE  LATERAL 

BEAM 

PARAMETERS 


'■   DELX,ILB, 
ALB,N2LRU,N2LRL, 
DELXP,  T2 


\ 


~7 


DELX,ILB, 
;ZVLB,N2LRU, 
N2IRL, DELXP 
T2 


Figure   20.    INPTBD  Program  Flow  Chart  (Continued) 
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READ  EFFECTIVE 
LATERAL  BEAM 
PARAMETERS 


/_ 


s 


\ 


DELX,  IT, 
ET 


AL. 


WRITE  OUT 
EFFECTIVE 
LATERAL  BEAM 
PARAMETERS 


V 


■i 


\ 


/ 


DELX,  IT, 
ET 


>f- 


INPUT  OUTPUT 
SPECIFICATIONS  AND 
SELECTED  POINTS 


k 


IOPT 

NPTS 
PTS 


1-4 


V- 


WRITE  OUT  OUTPUT 
SPECIFICATIONS  AND 
SELECTED  POINTS 


IOPT 
NPTS 
PTS 


1-4 


RETURN 


Figure    20.   INPTBD  Program  Flow  Chart  (Continued) 
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INPTBD 


DATE  =  73170 


00/31/2 


0001 
00  a2- 
0003 
£004- 


SUBROUTINE  INPTBD 

REAL  NUCNUS,  IUR,ILR,  IGO , IRGO, I NRT, I  LB, I T 

REAL  MD,  MDT 

LOGICAL-L-AST 


000  5 


0006 


0007 


-0008  -- 


1 
-2- 

3 

4- 

5 

-6- 


COMM  ON  /  BLK 1  /G-LftWCL-RLK-CLT-r-CR L-Or-GR-UD  ,    C2LRLt    C2LRU*-DR,    DRT  ♦ 
EC,    ES,     ILSC,    NLRL,    NLRU,     Tl,     IL,    NUC,    NUS, 
NRUUOt,    NRL(IO),    B(10),     ISC(IO),    ALB(IO),    N2LRU(10), 
N2LRL(10),    T2(10),     IT(10),     IUR(IO),     I LR ( 10) , ILB( 10)  , 


7 
-8- 

9 


£ELXP(10),     ET-(l-0), 

AI(3,10),    01(3, 10),    CT(3,10),    CCP(3,10),    TCP(3,10),    BCP(3,10), 
-ORG  (3,1  a)-,-  C  13 ,-UD  ,-  ISQ4-3+-1Q.1-,  -INRT(3,ia) ,     EBAR(3,10), 

BG(3,3,10),    CG(3,3,10),    TG(3,3,10),    NGRU( 3, 3, 10) ,    NGRL (3 ,3 , 10 ) , 

-CRURIS, 3, 10),-  CRL(  3,3,10),    IRG0(3,3,10) 

,DUM( 1568) 


IN,     10,     ISS,     ITER, 

T,    WEIGHT,    XLB,    YBCL  , 


COMMON    /BLK2/    DT,  THRSH,       IE,     IEVENT, 

— 1 NAP,  -Nft,— NO^-EH-,    NL-,    NSI  r  NSUPr    STEADY, 

2    NXPTS,  NDBUG(IO), 

—3- DELW-(IO)^  DELX(IO),  EEI(10,10),  EIL(IQ),  IS(10),  NCS(IO), 

4    XSUP(IO),  YB(10),  YP(4),  YCL(IO),  ZCL(10,10),  ZCL(10,10), 

6- -EI 3,1-0),  EI  13, 10),  XCS(3,10),  ZBL(3,10),  ZBU(3,10), 

7    F(80,10),  MD( 11,20,10),  MDT(11,10),  WP(50,4),  X(80),XP(50) 


( 


COMMON  /BLK5/  S IG ( 80) , DFLN ( 80  ) , DFL( 80, 10) ,  IT  IME,KTIME  , 

_0 SI-GP1120.,3,10)  ,  SIGC(  10,3,  10),  SIGTHS,  DSIG 

2      ,  I0PT1,  I0PT2,  I0PT3,  I0PT4,  NPTS(IO),  PTS(3,10) 
^--   --,  IPTS(3»10) 

— NAMEL  IST74^IPJJ-T-A/J(La,-NL^YaCL»  ... .  EC,ES,-NUC,  NUS,  DT,  THRSH,  WEIGHT, 

1  DR,     CRUD,CRLD,NB,NSUP,C2LRU,C2LRL, DRT, NLRU, NLRL, CLRU,CLRL, 

2  ILSC,  Tl,  IL,  SIGTHS,  DSIG,  NDBUG 


0009 

C 

DIMENSION  IB(10) 

-00  to 

0011 

DO  20  1=1.10 

20  NDBUG (I )  =  0 

. 

C 

0012 

IN  =  5 

O013 

-  10  =  6 

0014 

LAST=.TRUE. 

-QQ4  5- 

-WL ---  2 

0016 

C 

ND  =  1 

0017 

READ  (IN,INPUTA) 

00 1-8- 

C 

WRITE  (lOtlNPUTA) 

— GONA/E-RT-UNITS  TO  £E-ET- 

C 

00 1 9  -  - 

- 

- EC  =  EC*144. 

0020 

ES  =  ES*144. 

0021 

OR  =  DR/12. 

0022 

CRUD  =  CRUD/12. 

0O2  3- 

-C-RL-D  =  CRLD/12.  -■-  ■■ 
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INPTBD 


DATE 


73170 


00/31/2" 


0024 
0025 
0026 
002  J 
0028 
-0O££- 


C2LRU  =  C2LRU/12. 
C2LRL  =  C2LRL/12. 
DRT  =  DRT/12. 
CLRU  =  CLRU/12. 
CLRL  =  CLRL/12. 
-T-4— =L-Ii/_2„ 


0030- 
0031 


DO  100  L  =  l,  10 

100  IT(L)=0. 

READ  CENTER  LINE  POSITIONS 


0032 
0H33- 
0034 
0035 
0036 
0037 


READ  (IN, 1000)  (YCL( J),J=1,NL) 
1QQ0  FORMAT  (2F10.4) 

WRITE  (10,2000) 
2Q00  FORMAT  ( 1H1 , • YCL ( J ) • ) 

WRITE  (10,2001)  (YCL( J), J=1,NL) 
2001  _E.QR.MAT  1_LX,?F  16-.4J 


0038 


_0_Q3___ 


_C READ  SUPPORT  POSITIONS 

C 

READ( IN, 3015)  ( XSUP ( L ) , L= 1, NSUP ) 
C 

_WR_IJE_(  10.2016)  ( X  SUP (L),L=1, NSUP) 


0040 
0041 


0042 


2016  FORMAT  ( 'OSUPPORT  POSI TIONS1 , 1 1 F10. 2 ) 
3015  FORMAT  (11F7.2) 


0043 

0044 
0045 
0046 

0047 

0048 


0049 
0050 

0-051 
0052 
0053- 


C 
C 


L._  READ  BEAM  PARAMETERS 

D0_  L  J=1»NB 

READ  ( IN, 1001)  NRU( J),NRL( J ) , I UR ( J),ILR( J ) , YB( J ) , B( J ) ,NCS ( J ) , ISC ( J 
1),IB(J) 
1001  FORMAT  (2I3,2E10.3,2F8.3,3I1) 
WRITE  (10,2002) 

2002  FORMAT  ( »0   BEAM    NRU    NRL  IUR  ILR', 
1 i,, _.       YB  .          B    NCS    ISC     IBM 

WRITE  (10,2003  )  J ,NRU ( J ) ,NRL( J ) , I UR ( J ) , I LR ( J ) , YB( J)  ,B ( J ) ,NCS ( J  )  , 
1    ISC(J),  IB(J) 

2003  FORMAT  (3 16 ,2E 16  .3 , 2F 12.3, 3 16 ) 

CONVERT  UNITS  TO  FEET 


IUR(J)  =  IUR( J  )/20736.0 
ILR(J)  =  ILR( J)/20736.0 


__  NCSJ=NCS(J) 

READ  (IN, 1002)  ( XCS ( N  ,  J  )  ,N=  1 , NCS J ) 
lQ0-2_fDRMAI.  L3F8.3)   ._.  - 


00  54 
0055 
0056 
0057 


WRITE  (10,2004) 

2004  FORMAT  ( 1 X,  '  XCS (N, J ) • ) 

WRITE  (10,2005)  (XCS( N, J ) , N=l, NCS J ) 

2005  FORMAT  (3F12.3) 


0O5__-- 


D.CL2  N  =  1,NCSJ 
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0059  CI(N,J)  =  0.0 

0060  _..   CM*  J)  =0.0 

0061  IF  (IB(J).NE.O)  GO  TO  5 
X-  -  -  -  -  -  . 

C    REAO  ALL  BEAM  PARAMETERS 
-C- 


0062  IF  ( ISC( J ).E0.0)  GO  TO  3 

C    STEEL  BEAM 
JC 


0063               REAO  (IN, 1003)  AI (N , J ) , IGO ( N, J ) t CI ( N, J ) , CT( N, J ) ,CCP{ N , J ) 
1 ^-JCP(N,J.)  ,..  BCPXN,.J-L.  _  


0064          1003  FORMAT  (7F8.3) 
_006S MRI  T.E  .  (  10 ,  2006  ) 

0066  2006  FORMAT  ('0  AI  IGO  CI  CT« 
__1 ,.«. _    CCP          TCP          BCPM 

0067  WRITE  (10,2007)  AI  (IM,  J  )  ,  IG0(  N,  J  ) ,  CI  (  N,  J  )  ,CT(  N,  J  )  ,  CZPi  N  ,  J  ) 
— — 4— »  JC P  ( M  ,  J4 -»-  BCP  ( N  ,-J  )  -  —..--- 

0068  2007  FORMAT  (1X,7F12.3) 

C         CONVERT    UNITS    TO    FEET 

0069  AI(N,J)    =    AI (N, JJ/144.0 

0070 T  GO  ( N  , . I )  _=  -IGQlIN.,  J  ) V  20 7  36  .  0  .  „ _ 

0071  CI(N,J)     =    CKN,  J)/12.0 

01122 CJ  (N  ,JJ    =   CT ( N , J ) / 1 2 . 0 

0073                                      CCP(N,J)     =    CCP(N,J)/12.0 
1)074 TCP  IN,  J)    =    TCP(N,J)/12.0 

C 
0075 GO    TO _2 

C 

C_        READ    EFFECTIVE    BEAM    PARAMETERS 

C 
_0£L7_6.....  .  5  J^EAD     (IN, 1007)     INRT  (N,  J  ) ,  EBAR(  N,  J  )  ,C  (  N,  J)  ,  E  (  N,  J  ) 

0077  1007    FORMAT    ( 2E10.3 ,F8 .3 ,E 10. 3 ) 

-007-8 iLRII.£_tIO,20L4)     

0079  2014  FORMAT  ( IX, • INRT ,EBAR ,C, E « ) 

0080_  WRITE  (10,2015)  INRT(  N,  J  )  ,EBAR  (  N,  J  )  ,  C(  N,  J)  ,  E  (  N,  J  ) 

0081  2015  FORMAT  ( IX, 2E 16.3 ,F12. 3,E 16. 3 ) 

0082 INRT(N,J)  =  INRT (N , J ) /20736  .  0 

0083  EBAR(N,J)  =  EBAR (N , J )*144.0 

00-84- C  Ul  r  J4-  -  -C-(  N-,-J  )  / 12  .  0 

0085  E(N,J)  =  E(N,J)*144.0 
C 

0086  GO  TO  2 
C 

C    CONCRETE 

0087  3  CONTINUE 

0088  WRITE  ( 10,2010) 

0089  2010  FORMAT  ('0  BG  CG  TG   NGRU   NGRL«. 

1    •  CRU  CRL         IRGO') 

C 

-Q09^ DO  4  M=U3 
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0091 

0092 
00-93- 

-0O94- 


READ  (  IN, 1005  )  BG(M,N,J),CG(M,N,J),TG(M,N,J) ,  NGRU(M,N ,J)  ,NGRL(M,N, 
1J),CRU(M,N, J),CRL(M,N, J) , IRGO(M,N, J) 
1005  FORMAT  ( 3F8.3 , 2  I  3  ,3F8 . 3  ) 

.  WRITE  (10,2011)  BG(M,N, J),CG(M,N, J),TG(M,N, J ) ,NGRU( M ,N , J ) ,NGRL(M,N 

1,J ),CRU(M,N,J ),CRL(M,N,J),IRGO(M,N, J) 
-^Q-ll  FORMAT  -i4X^-SF42^at2-^6t3F-12^3) 


CONVERT  UNITS  TO  FEET 


-00-95 

0096 
0097 


0098 

0099. 

0100 


BGLM,N,J)=  BG(M,JM,  J)/12.0 

CG(M,N,J)  =  CG(M,N,  J)  /12.0 
IG-Ui»ISUJI-=- XGIM-,JvLmD/  12.0-  _.._  - 

CRU<M,N,J)  =  CRU(M,N, J)/12.0 
CRL(M,N,J)  =  CRL(M,N,J)/12.0 

IRGO(M,N,J)  =  lRG0(M,N,J)/20736.0 


0101 


0102 

0103 

0104 

0JL05_ 

0106 

0107 


4  CONTINUE 


READ  (IN, 1006)  DRG  (N,J) 
.1006  FORMAT-  .(..F8.3) 

WRITE  (10,2012) 
20.12-  FORMAT  HXt'DRG'l 

WRITE  (10,2013)  DRG(N,J) 
2013  FORMAI— L1X_,-EJ.2.-3J  .._.  


0108 


0-10SL 


C CONVERT    UNITS    TO    FEET 

C 
DRG(N,J)     =    DRG(N, JJ/12.0 

C 
2  CnNJ-INllf 


0110 


0111 
0112 
-0113- 

0L14 
0115 
0116 
0117 
■0418-- 

0119 
0120 

0121 
0122- 


0123 
0124 
-0-12-5- 


1    CONTINUE 


_. .  C  __ 

c 


IF(  IL.EO.O) 
IF  (IL.NE.O) 

2018- FORMAT    (LOS 
1  •       N2LRU 

2020    FORMAT     ( '0', 
NSI=NSUP-1 
DO    6    L=1,NSI 
IF     (IL  .NE, 


WRITE 
WRITE 

N2LRL 


(10,2018) 
(  10,2020) 

DELXP 


DELX 


DELX 


ILB 


T2'  ) 


IT 


ALB 


ET»  ) 


0)     GO    TO    7 


-READ-  UN-,  1009)  -       DELX(L), 

1  ILB(L),ALB(L),N2LRU(L),N2LRL(L),DELXP(L),T2(L) 

1009    F0RMAT(2F8.3,E10.3,  2I3,2F8.3) 

WRITE     (10,2019)  DELX(L), 
..._.       1                                              ILB(L)  ,ALB(  L)  ,N2LRU(L)  ,N2LRL(  L)  ,DELXP(  L)  ,T2(L) 
2019    FORMAT    ( IX , IF  1 2 .3 ,E 16 . 3, F 1 2.3, 2 1  6, 2F 1 2. 3 ) 
1T-1L)  -=  --O— — 


C 
C 

c 


CONVERT    UNITS    TO    FEET 

ALB(L)    =    ALBID/144.0 
ILB(L)     =    ILB(L)/20736.0 

0  E4XP  14-1  -   DEL  XP  (  L  )  / 1 2.  0 
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0126  T2(L)  =  T2(L)/12.0 

--  -C—  

0127  GO  TO  6 


0128  7    READ    (IN, 1010)  DELX(L),     IT(L),     ET(L) 

0129 101#-FeRMA-T— (—F8-.3r2€-10^3)- — -■  

0130  WRITE    (10,2021)  OELX ( L ) , I T( L ) , ET( L) 

^0ia4- -2021 -FORMAT    (IX,    F12.3,2E16.3) 

C 
C CONVERT  ONT-T-S   T 0  -FEET 

C 

0132 ITtL)    =    IT(L-1^20736.0 _  _ 

0133  ET(L)    =    ET(L)*144.0 


0134  6    CONTINUE 

C         READ    OUTPUT    SPECIFICATIONS 

0135  READ    (IN, 3001)     I0PT1,     I0PT2,     I0PT3,     I0PT4 
D13i> 30£  1    FnRMAJ    14  ID 

0137  WRITE     (10,3002)     IOPTl,     I0PT2,     I0PT3,     I0PT4 

D138 300.2- FQBM AT    L1H0,     'IOPTl       I0PT2       I0PT3       I0PT4'  ,  /,  1  X,4I6  ) 

C 
0139 IF    (  IOPTl  .FGUIl^AHa^IO£T2^£CU0<AND.IOPT3<EQ.0)       GO    TO    60 

C 

_Q1_4Q DD_5_Q  _I  =  UNB     _ 

0141  50    READ    (IN, 3004)    NPTS(I),     ( PTS( J, I ) , J=l, 3 ) 

0142 30fl4_FJ3RMAT    U 1 ,3F8  .3  ) 

C 

.01.43 WRIIEJ  TO, 3005  I 

0144  3005  FORMAT  ('ONPTS      POINT  1      POINT  2      POINT  3   ', 

_  1      •    INPUT   POINTS'//) 

C 

0145 DO  55  I=1,NB 

0146  NJ  =  NPTSU  ) 

0147 5-5_WiU  IE  - 1  la,  30061 -NP.T.&U-U  ( P  TS(  J , I } , J= 1 , NJ ) 

0148  3006  FORMAT  (IX,  14,  3F12.3) 

0149  i>0  CONTINUE 
C 

0150  RETURN 

0151  END 
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GRDNET 

The  GRDNET  routine  establishes  the  grid  point  values  X  and  the 
number  of  longitudinal  grid  points  per  span  IS  (L) . 

The  number  of  points  is  first  established 

XSUPT  .n  -  XSUP0 

IS,   =      L+1       2 

jj  ■  ' 

If  ISLNB  is  greater  than  80,  ISL  is  recalculated  to 

IS,   _   80 
L   "   NB 

and  a  new  grid  incremental  distance  is  calculated 

XSUPL+1   "   XSUPL 


AXL   ■         ISL 

J-l 


A  grid  weight  is  calculated  as 


.  „  \  XLB  -  Weight 

Li 


/xsupl+1  -  xsup2\ 

\  XLB  -  Weight    / 


££) 


Grid  positions  are  calculated  as 


X.   =  XSUPT   +  At-       +   (i-D  Axt 
i  L       2  L 


Output  points  are  converted  to  the  nearest  grid  point. 
No  subroutines  are  called  by  this  routine. 
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Igrdnet  ) 


xo  =  o. 


Io=  ° 

II  =  1 

L  =  1 

ksupo=  0 


I     -  XSUP(£)-XSUPO 
S(L)       W. ' 

AX(L)=  XSUP(C,)-XSUPQ 

f      

.  S(D) 

AT,T(T^-pSUP(L)-XSUPO 
V        XLB 


± 


+  1 


■) 


W 


x.(n.)NB 


X(II)  =  X   +  AX(L)/2 
o 

*END    1s(L)  +  I0 


I 


I   =   II+l 


B 


Figure   21.    GRDNET  Program  Flow  Chart 


80 


^ 


X(I)  =  Xjj.^  +  d-DAX^j 


Figure   21.   GRDNET  Program  Flow  Chart  (Continued) 
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No 


x' 

II  =  II+l 

Xo  =  XSUP(L) 
XSUPO  ~   XSUP(L-1) 
Io  :      IEND 

__2 

r 

^ 


Yes 


Yes 


Write  out 
grid  data 


"7 


L,IS,AX, 
AW,  X. 


Convert  Output 
points  to  the 
nearest  grid 
points 


Write  out  con- 
verted points 


~7 


NPTS,  PTS, 
IPTS 


(    RETURN    J 


Figure  21.     GRDNET  Program  Flow  Chart  (Continued) 
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:ORTRAN  IV  G  LEVEL   21 


GRDNET 


DATE  =  73170 


00/31/2 


0001 


0002 


C 
G_ 


0003 


0004 


0005 
4)006 


SUBROUTINE  GRDNET 
THIS  ROUTINE  ESTABLISHES  X  GRID  POINTS  AND  ISL 
REAL  MD,  MDT 

THRSH, 


IN,  10,  ISS,  ITER, 

T,  WEIGHT,  XLB,  YBCL , 


COMMON  /BLK2/  DT,       THRSH,   IE,  IEVENT, 
_^_-NAP.,  NB,  NO,  NI,  NL,  NSI,  NSUP,  STEADY, 
2    NXPTS,  NOBUG(IO), 

3-  ---DEL-W410)-,  DELXUO),  EEI(10,10),  EIL(IO),  IS(10),  NCS(IO), 

4    XSUP(IO),  YB(10),  YP(4),  YCL(IO),  ZCL(10,10),  ZCU(10,10), 

6 E(3,10)r  EI  (3,10),  XC-S-t-3-^10) ,,.  ZEL(3, 101^  ZBU(  3, 10  )  ,  _ 

7    F(80,10),  MD( 11,20,10),  MDT(11,10),  WP(50,4),  X(80),XP(50) 


C- 


COMMON  /BLK5/  S IG ( 80) , DFLN ( 80 ) , DFL( 80, 10) , IT IME ,KT IME , 

-1 SIGP(120,3,10W  SIGC( 10,3, 10),  SIGTHS,  DSIG 

2      ,  I0PT1,  I0PT2,  I0PT3,  I0PT4,  NPTS(IO),  PTS(3*10> 
-3 ,  IPTS<3t4-Q-) 


DO    20    1=1,100 
2.Q  XX 14    =    0.0 


0007- 


-32-2- 


0008 
-QO09 
0010 

0011 

-0012- 

0013 

0014 
-0015 
0016 
0017 


_x_ 


XO    = 

_L0_  = 

II    = 


0. 

0 

1 


0018 

0019- 

0020 


DO    1    L=1,NSI 

_DJCSUe_=— JK.SU E-iL  +  1)    -_XSI)R.(U_. 
IS(L)    =    DXSUP/DELXU) 

ITST  =   IS(L)*NB 
JFiJ.TST.LE.80)   GO  TO  14 

ITST  =  80/NB 
D-EL-X-t-LJ-  =-  0-XSUP/4-TSI 

IS(L)  =  ITST 

WRITE  (10,100) 
100  FORMAT  ( 'OINPUT 
.   ._l i    FOR  L  = 


DELX(L),L 

DELX  TOO  SMALL  - 

•,I2) 


DELX  CHANGED  TO' ,F9.3, 


0021 


0022 


14  GONT  1NUE-- 


DELW(L)  =  DXSUP/XLB*WEIGHT/( IS(L)*NB)/G 


0023 
0024 


X(  II )  =  XO  + 
I  END  =  IS(L) 


DELX(L)/2. 
+  10 


-0025- 


4-L-^ 


0026 
002  7 
0028 
0029 
0030 
-00-34- 


DO  2 
X(4  ) 
II  = 
II  = 
XO  = 
10  - 


Air*       1 

I=IL,IEND 

=  X(I-l)  +DELX(L) 

I 

II  +  1 

XSUP(L+1) 

I€NO 
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0032 

0033 

£034- 

0035 

0036 


1  CONTINUE 


IF  (NDBUG(l).EQ.O) 
.WRUg_410-flO)-  — 
FORMAT  (    • OGRDNET 
00  3  L~1,NS-1 


GRONET 


GO  TO  15 


DATE 


73170 


00/31/2: 


10 


0037 
0038- 

-0£39~ 
0040 
0041 


WRITE  (10,11)  L,IS(L),  DELX(L),  DELW(L),  XSUP(L),  XSUPU+1) 
FORMAT  (IX, • NO.  SPANS=» ,I3,2X,'N0.  GRIDS=',I3,«   DELX= • ,E12 .4, 


3 
11 

1  •       DELW  =  »,E12.4,»       SUPPORTS       S2E10.3) 

_4.5— CONTINUE 

WRITE    (10,12) 
12    FORMAT    ( Lq  GRID    CENTERS- 1  0NGLIUD-1-NA1.UC  _=  Jt!) 


0042 
0043 


L3 


WRITE    (10,13)    X 
FORMAT  XIX,  LOE  12*4  J 


CALCULATE    TOTAL    X    POINTS 


0044 
0045 


0046 
0047 


NXPTS    =    0 
nn    lfl    1  =  1  »MS1. 

NXPTS    =    NXPTS 
JL8_  XONXi  NU  E 


+    IS(L) 


0048 

0049 
0050 


CONVERT    SAMPLE    POINTS    TO    NEAREST    GRID    POINT 

IF     (I0PT1.EQ.0.AND.I0PT2.EQ.0.AND.I0PT3.EQ.0)       GO    TO    60 

DO    30    J=1,NB 

NPT  _=    NPTS(.l) - 


0051 
-QH52— 

0053 
-Q054— 

0055 
£054. — 


IF    (NPT.EQ.O)       GO   TO    30 
_D0    25    1  =  1, NPT 

DO    21    K=l, NXPTS 

IF     (PTS( I  ,J).EO.X(K) )       GO    TO    24 

IF     (PTSU,  J).E0.X(K).OR.PTS(  I,  J).GT.X(K+1)  ) 
__GO_.JO  22 -   -      .-        -       - 


GO    TO    21 


0057 
QQ58 
0059 
CI  060 
0061 
0062- 


0063 
0064 
0065 

£066 
006  7 

-0£6£- 
0069 


21  CONTINUE 

22  XT1    =    PTS(I,J)    -    X(K) 
XT2    =    X(K+1)     -    PTS( I, J) 

-  IF     (XT2.GT.XT1)       GO    TO    23 
PTS( I  ,J)    =    X(K+1) 

LP_XS(-i^J4— =-    K+i        

GO    TO    25 

23  CONTINUE 
PTSU, J)    =    X(K) 

-2-4  CONTINUE 

IPTS( I, J)  =  K 
25  CONTINUE 


30  CONTINUE 


WRITE  OUT  CONVERTED  POINTS 


0070 
££7-1- 


WRITE( 10,3005) 
3005-EQ&MA-T-  1-IONPTS- 


POINT  1 
84 


POINT  2 


POINT  3 
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1      •   CONVERTED  POINTS'//) 
C    


0072  DO  55  J=1,NB 

£0-73 NJ— =-   NPTS(  J) 

0074            55  WRITE  (I0»3006)   NPTS(J),  ( PTS( I , J ) , I = 1, NJ ) 
-0O7-5 300  6  FORMAT  ( IX-r -14- r- 3£42-.3^ 

0076            60  CONTINUE 
-00-73 WRITE  (-10^3007  ) 

0078               WRITE  (10,3008)   IPTS 
-00*9- 30&7— FORMAT  {J  04  PTS'  * - 

0080  3008  F0RMAT(2X,3I4,6X,3I4*6X,3I4,6X»3I4,6X,3I4) 


0081                                       RETURN 
_00£2 £HD 
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MODULI 

The  MODULI  routine  calculates  the  effective  longitudinal 
parameters : 

INRT   Effective  Moment  of  Inertia 

EBAR  Effective  Young's  Modulus 

C     Effective  Neutral  Axis 

E     Steel/Concrete  Modulus 
from  the  input  data  for  Steel  or  Concrete  longitudinal  data  or 
accepts  the  effective  values  as  input. 

Calculates  the  upper  and  lower  cross  sectional  moduli 


(EBAR   A    /INRT   A 


ZBL   .   =    ZBU 
n,D  n,} 


Calculate  the  effective  lateral  parameters 

EIL   Effective  Young's  modulus  times  the  effective 
moment  of  inertia 

XIL   The  effective  moment  of  inertia 

from  the  input  lateral  beam  parameter  or  the  effective  lat- 
eral beam  parameters. 

The  calculated  parameter  and  intermediate  values  are 
printed  out  if  the  debug  printout  is  specified,  i.e., 
NDBUG(2)  =  1. 

No  subroutines  are  called  by  this  routine. 

86 


(      MODULI     J 


<4> 


J    =    1 


N   =    1 


No 


<£ 


A2    =    NRUj    »(S|)' 


A,   =   TCP, 

„     .TBCPn., 

No 

/CCP(N,J) 

N.     =  o 

3 

es 

a            r 

J                                              ll    /   W                        **  /  w 

A3   "   ° 

/ 

> 

/ 

% 

* 

A4   "   BJT1 

p.       1 

< 

^ 


Figure     22. 


MODULI  Program  Flow  Chart 
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Yes 


Concrete 


A5"AIN,J 


Steel 


a  =  n( — ^)  y  (ngrum  „  _  +ngrl  m  "  ) 

5     V    2   /  ^±\  M,N,J,      M,N,jy 

+  m§i  (TGm,Nf JBGm/N,j) 


m^^^^m^^^^m.N^^Um^j; 


y   (tg  m 

nk  V  m'N' 


JBGM,N,JCGm,N,J 
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CC,  =  CIM  , 
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E    A*cc 


CP,  = 


CP„  = 


K=l 


K   K 


5 

E   A 

K=l 


K 


CTN/J-CP1 


ID  = 

EID  = 

IR  = 

EIR  = 
ICP  = 


IYBJ  Ti  +  bj  Ti  (CPi-K)2 

EC  *  ID 

A2 (CP1-CC2) 2+Ax (CP1-CC1) 2+NRUJIRUJ+NRLJIRL^ 

ES  *  IR 

I±(BCPN:_TTCP3  .T)+BCPM  ,TTCP,T  T(CCP,T  T-CPn)2 


n,j   n,j 
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IRG  = 

EIRG  = 
ICG  = 

EICG  = 

IG  = 

EIG  = 


3-    /wV^cp1- 


M=l  X 


RUM,N, J>  2nGRUM,N, J+  <CP1-CELM,N, J»  2NGI 


J      r  i  3 

E.   I2"BGM,N,JT1 


ES*IRG 
3 


+BGM,N,JTGM,N,J(CP1 


"7TGM/N/J)  J 


EC*ICG 

IRG+ICG 

EICG+EIRG 


V 


IBL 
EI 


ZBU 
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6 


■zk. 
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4     L 
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o 
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6 
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6 


y 


6 
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ALAT, 
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6 
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AL 
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EILDC=EC*XILDC 
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XILD=XILDC+XLR 


Figure   22. 


MODULI  Program  Flow  Chart  (Continued) 
94 


Yes 


Steel 


No 


Concrete 


EILB=ILBTES 

Li 


XILBR=ALAT_ (CP1L-C2LRU) 

3 

EILBR=ES*XILBR 

XILBC=ALAT4/12  (T2L)  2+ALAT4  (CP1L-ALAT4  j2 

EILBC=EC*XILBC 

EILB=EILBR+EILBC 

ILBT=XILBR+XILBC 

Ll 


V 


EILL=EILB+EILD 

XIL=ILBT+XILD 
Li  Li 


K 


L=L+1 


Yes 


No 


<£> 


WRITE  OUT 

LATERAL 

PARAMETERS 


C    RETURN    J 


A 


7 


EIL,XIL,CP1L, 

CLAT,ALAT, 

XILOC,EILDC, 

XLR,EILR,EILD 

EILB,XILBR, 

EILBR,XLBC, 

EILBC,ILB 


Figure  22.    MODULI  Program  Flow  Chart  (Continued) 


95 


ORTRAN  IV  G  LEVEL   21 


MODULI 


DATE  =  73170 


00/31/2"/ 


0001 


0002 


0003 


SUBROUTINE  MODULI 
X-    THIS  SUBROUT  RECALCULATES  MOMENTS  OF  INERTIA,  NEUTRAL  AXES, 

C      EI  AND  CROSS  SECTIONAL  MODULI 
REAL  NUC,NUS,IUR,ILR,IGO,IRGO,INRT, I  LB, I T, ID, ICP, IG , IRG, ICG , IBL , IR 

C 
REAL  MD,  MD-T -   _ _.      


0004 


0005 


REAL  IBL 

C0MM0N/BLK1/CLRL,CLRU,  CLT,  CRLD,  CRUD,  C2LRL,  C2LRU,  DR,  DRT , 

-1 EC,  ES-f-  ILSC,  NLRL,  NLRU,  Tl,  IL,  NUC,  NUS, 

2         NRU(IO),    NRL(IO),     B(10),     ISC(IO),    ALB(IO),    N2LRUU0), 
3 M?  I  Rl.  (  1  0  )  j_ I  ?  (  1  0  L»— LXLIQlU.  1  UR.L1  OJL,_J  LR  (  LO  )  ,  I  LB  U  0 )  ,_   _ 


4         DELXP(IO),     ET(10), 
_5 AI.(3«lQ)j    CI(3,10),    CT(3,10),    CCP(3,10),    TCP(3,10),    BCP(3,10), 

6         DRG(3,10),    C(3,10),     IG0(3,10),     INRT(3,10),     EBAR(3,10>, 

7__BGt3t3,10),    CG(3,3,10),    TG(3,3,10),    NGRU(  3,  3, 10)  ,    NGRL  (3  ,3  ,  10  ) , 

8         CRU(3,3,10),    CRL(3,3,10),     IRG0(3,3,10) 
-9 ,DI.IM(ltiftR) ...- 


DQ06 


jQQQl 


POOR 


_£DMMDR 7BLK2/    DT,  THRSH,       IE,     IEVENT,     IN,     10,     2SS,     ITER, 

1         NAP,    NB,    ND,    NI,    NL,    NSI,    NSUP,     STEADY,     T,    WEIGHT,     XLB ,    YBCL ♦ 

2_    MXPTS,    NDBUG(  10), 

3         DELW(IO),    DELX(IO),    EEI(10,10),     EIL(IO),     IS(10),    NCS(IO), 
A XSUPtlO).    YB(IO).    VP(4).    YCL(  10)j  .ZCL(  IP,  10)  ,     ZCU(.  10, 10  ) ,  . 

6         E(3,10),    EI(3,10),    XCS(3,10),    ZBL(3,10),    ZBU(3,10), 
.7 FL8.0_»..10),.  MD111,20,1Q).,    MDT(11,10),    WP(50,4),     X(80),XP(50) 

__GOMMON  7BLK6/_   IBL  <  3, 10) 

DIMENSION    CC(5).A(5).ALAT(6).CLAT(6),     XIU10) 


£0D9_ 

QQJJL 

0011 

001?  ._ 

0013 

0014 

0015 

0016 
0OJ.X- 
0018 
0019 


0020 


EQUIVALENCE  (CPLAT,CP1L) 

DO  30  J=l,10 
DO  29  1=1,3 
_LBLLL,-JJ_=  .0.0..  


C 


29  CONTINUE 
XIL(J)  =  0.0 

30  CONTINUE 

IF  (NDBUG(2).EQ.O)   GO  TO  40 

lOO-JJIRMAJL.t .■■-  '0M0DULI1 ) 

WRITE  (10,100) 
40  CONTINUE 

CALCULATE  LONGITUDINAL  PARAMETERS 

PI=3. 141593 


0021 
0022 

0023 

0024 


DO  1  J=1,NB 
NU=NCS( J) 

DO  2  N=1,NU 

J.F_ICIN+J1.NE*0)  -GO- TO  13 
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0025  TMP    =PI*DR*DR/4. 

0026 JU1)=NRL(  J  )*TMP 

0027  A(2)=NRU<  JUTMP 
C-      ._    - 

0028  IF     (CCP(NtJ).EQ.O)    GO    TO    3 

-0O2-9 A-(-3-UTCP(N,  J-U-&CP  UU-  J  ) 

0030  GO    TO    4 

0034 3-A-<3-U0-. 

0032  4    A(4)=B(J  )*T1 

0033  IF     ( ISC( J  ).NE.O)    GO    TO    6 

C         CONCRETE 
C __. 

0034  TERM1=0. 

0035_    I£RM2  =  Q  . 

0036  TERM3=0. 

0037 T-E-RM4=0. _ 

0038  DO    5    M  =  l,3 

£113.9 I£RM1  =  TERM1+NGRU(M,N.,  J)+NGRL(M,N,J) 

0040  TERM2=TERM2+TG(M,N, J)*BG(M,N, J) 

0041 .    _.TERM3  =  TERM3+CRL(M,N,  J)*NGRL(M,N,  J)+CRU<M,N,  J)*NGRU(  M,N,  J) 

0042  5    TERM4=TERM4+TERM2*CG(M,N, J) 

004  3 TMP=PJ_*D&G11>U  JJL/2.J*CLRG  (AL,  J  )  LZ* _  

0044  A(5)=TMP*TERM1+TERM2 

00:45. C_I  (_N_iJ  )  =  (  TMP#TERM3  +  TERM4  )  /A  (  5) 

0046  GO    TO    7 


C         STEEL 
X 


0047  6    A(5)=AI(N,J) 
-    C. 

0048  7    CONTINUE 

004,9- _       CC(1)=CRLD 

0050  CC(2)=CRUD 

0051 CC4a.UCCPi.ALf  JJ  -  - - 

0052  CC(4)=Tl/2. 

0053  CC(5)=CI(N,J) 

0054  CP1=0. 
G055     -                 — -  DEN  =  0. 

C 
£0£6 DS— 5^  K=W5- 

0057  DEN=DEN+A(K) 

0058  55    CP1=CP1+A(K)*CC(K) 
C 

0059      _  CPUCP1/DEN 

0060  CP2=CT(N,J)-CP1 

-0O61 I-D=B-U)V-12^.*Tl*Ti*Tl+-B(JJ*Tl*i  CP1-T1/2.  )**2 

0062  EID=EC*ID 

0063    1R=A12)*(CP1-CC(2) )**2+A(l)*(CPl-CC(l) )**2 

1+NRU( J )*IUR(J)+NRL(J)*ILR( J) 

0064  EIR=ES*IR 

0065  ICP=BCP(N,  J)/12.*TCP(N,J)**3+BCP(N, J)*TCP(N, J)*(CCP(N,J)-CP1)**2 

.0&66 £4  CP  =  ES  *.  1C  P 
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0067 


0O6S- 


0069 
OQJO 


IF     USC(J).EQ.O)    GO    TO    8 


-STEEL 


-1G~A(5)*(CC451  =CP-1-X±*2i  IGO(Nt  J )_ 

EIG=ES*IG 
GOTO    9        - 


-CONCRETE 


0071 


0072 
0073 


0074 
0075 


0076 


AQ77 
0078 
0079 


0080 


a    TMP=PT*r)Rr,(Nf.l)  /7.*nRG  (Nf  .1)  /?. 


IRG=0. 
ICG=£U _     _.. 

DO    10    M=l,3 

IJlG_=IRG+(  JCPl-CRU(MtN,  J)  ) **2*NGRU( M, N, J ) +( CP1-CRL  (  M ,N , J )  )**2*NGRL 

1(M,N,J)+IRG0(M,N,J) )*TMP 
10    I r. G  =  I C G ± B G ULslU^l UAZ**I±** 3+EHLM, Jj,  ll*TGlM_iN-*J)*(XPl-TG(  M_iN,  J  )72.  ) 


EIRG=ES*IRG       . 

EICG=EC*ICG 
_1G=JRG+ICG 
EIG=EICG+EIRG 


0081 
0082_ 
0083 


IBL(N,J)=ID+IR+IG+ICP 

EI(N,J)=EID+EIR+EIG+EICP 
ZBU(N,J)=       EI(N,J)/(EC*CP1) 


00  84 


IF     (  ISC(J) .NE.O)    GO    TO    11 


0085 
00_8i> 


CONCRETE 

ZBL(N,J) 
_     GO   TO    12 

-STEEL 


=    EI  (N,J)/(CP2*EC) 


0Q.87 
0088 
QQ89 
0090 
009JL 


0092 
009J3. 
0094 
0Q95_. 


11  ZBL(N,J>     =    EIINi  J)/(CP2*ES) 

12  GO    TO    33 

13  ZBU(N,J)=EBAR(N, J)/C(N, J)*INRT(N, J)/E(N, J) 
ZBL(N,J)=ZBU(N, J) 

IBL(N,J)  ..=  .  INRT  (N»J)__ 

EI(N,J)     =    EBAR(N,J) 


33    CONTINUE 

IF     <NDBUG(2) .EQ.O) 
100Q    F0RMAT(1X,10E13.4) 


GO    TO    14 


0X1 96- 


0097 
0098 
0099 
01Q0 
0101 
0102 


MATTE     (ID  ,1 OQJL) 

1001  FORMAT     (  'OA'  ) 
WRITE     (10,1000) 

WRITE    (10,1002) 

1002  FORMAT    ( • OC I ' ) 
WRITE     (10,1000) 

WJUTE    110,10031 


CI 
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0103 
0104. 
0105 
OlOiu 


0107 
.010-8- 


0109 


0110 
0111 


1003  FORMAT  ( • OCC • ) 
WRITE- ( 10,1000)   CC 
WRITE  (10,1004) 

1004  FORMA!  ( 

1    •           EIR 
2- i CP2' 


•0 


ID 


EID 


IR' 


ICP 


EICP 


CP1 


WRITE  (10,1000) 
___.  __WiUT£  (10,1010)   _. 
1010  FORMAT  (  '0         IG 
1 •--         ICG! L) 

WRITE  (10,1000)   IG, 
14  CONTINIIF 


ID,  EID,  IR,  EIR,  ICP,  EICP,  CP1,  CP2 


EIG 


EIG,  IRG,  EIRG,  ICG 


IRG 


EIRG' 


0112 
-Q1JL3-  . 

.0114 

2    CONTINUE 

1-CONTINUE.                        ......... 

C 

LF_  (_NDBUG_(2).EQ.O)       GO    TO    15 

0115 
0116 

WRITE    (10,101) 
101       FORMAT    (  '0          J         N 

IBL 

0117 

1          ,'                               ZBU 
no   2^    J=1 ,    NP 

ZBL'  ) 

0118 

0115 

0120 
ni  ?1 

DO    26    N=l,3 

WRITE    (10,102)     ( J,N,IBL(N,J), 
102    FORMAT    ( IX, 214, 4E 14.4) 
?h    r.nNTINtlF 

EI  (N 

0122 

C 

.15.  CONTINUE  ...       

EI' 


EI (N,J),ZBU(N,J),ZBL(N,J) ) 


-CALCULATE  LATERAL  PARAMETERS 


0123 


0124 
0125 
0126 
J3121 
0128 
0129 


0130 
Q13.L 


-Da-2fl...L=l»MSl 

IF  (IT(L).NE.O.)  GO 

ALAT(1)=T1*DELX(L) 

AR=DR/2.#DR/2.*PI 

ALAT(2)=NLRU*AR 

ALAT(3)=NLRL*AR 

CLATU)-=T-l/2. 


TO  21 


CLAT(2)=CLRU 
CLAT(3)=CLRL 


0132 


IF  ( ILSC  .EQ.O)  GO  TO  22 


-CGNCRE-T-E-  


013  3- 

0134 

01-35 

0136 

0137 

0138 

0139 


.QIAO-  


ALAT(4)=DELXP(L)*T2(L) 
ALAT(5)=N2LRU(L)*AR 
ALAT(6)=N2LRL(L)*AR 
CLAT(4)=.5*(T2(L)+T1) 
-CLAT(5)=-C2LRU 


CLAT(6)=C2LRL 
-GO.  TO  2  3 

STEEL 

22  ALAI14-)  =ALB-(-U- 
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0141  ALAT(5>  =0. 

0142-   ....  ALAT(6)=0.  . 

0143  CLAT(4)=.5*(T2(L)+T1) 

0444 —  C1AT(6J  =  0. 

0145  CLAT(5)=CLAT(6) 

0446^ — 2-3    AL=0. — .  --      -     - 

0147  CP1L=0. 

0448 DO    200    K= 1,6 

0149  AL=AL+ALAT(K) 

04-50 2Q0-  CP1L=CP1L±CLAT(K)#ALAT(K) 

0151  CP1L=CP1L/AL 

0JL52 X1LDC=DFI.  XXL.)  / 1 2  .  *X1**3+QELXJJ. 1  *  I.  1*  (  CP1 L- 11/  2  . .)  *  *  2 

0153  EILDC=EC*XILDC 

0154 XLR=AJ-AT(2)*(CP1L-CLRU)**2+ALAT(3)*(CP1L-CLRL)**2+N2LRU(L)*IUR(1  )  + 

1N2LRL(L)*ILR(1 ) 

0155 EILR=ES*XLR 

0156  EILD=EILDC+EILR 

ni57 xxlq=xii  nr.+xi  R __.. 

0158  IF     (ILSC       .NE.    0)G0    TO    24 


STEEL 


0159  EILB=ES*ILB(L) 

0160 GQ_KL2_5_ ..       _      ._ 

C 
C CONCRETE 

C 

0161 24_XILBR=ALAT(5)*(CP1L-C2LRU)**2+ALAT(6)*(CP1L-C2LRL>**2 

0162  EILBR=ES*XILBR 

0163 XILBC=ALAT(4)/12.*T2  (_L .)  **  2+ AL  AT  (  4 )  _  *  (  C  P 1 L-CL  A  T  (  4 )  )  **2 

0164  EILBC=EC*XILBC 

0165     EI  L  B  =  E  I L  B  R+ E  I L  BC 

0166  ILB(L)=XILBR+XILBC 
C_          ._ 

0167  25    EIL(L)=EILB+EILD 
0168 XJL_LLJ.=1LB(L)+X1LD _.   __ 

C 

_C CALCULATE    Z    -      MODULI    FOR    LATERAL    EFFECT 

0169  GO    TO    20 

017Q_      21    EIL(L)  =  IT(L)*ET(L) 

0171  XIL(L)=IT(L) 

0172 20- XQNX1NUE  __ .. 

C 

0173  IF     (NDBUG(2  ).EQ.O)       GO    TO    60 

0174  WRITE     (10,103) 

0175 1Q3    .FORMAT.   t.VQ         L  EIL  XIL') 

0176                                      WRITE     (10,302)     ( L ,E I L ( L  ) , X  I L( L ) , L= 1 , NSI ) 
0177 302     FOJlMA.T_-(  1  X  ,  1.4 ,_.  ,.2.E  12^3 ) __  _     

0178  WRITE     (10,1017) 

0179      1017    FORMAT     ('0  CP1L  CLATM 

0180  WRITE     (10,1000)       CP1L,    CLAT 
0-18L     .-.__    _     WRITE     (10,1018) 

0182  1018    FORMAT    (•OALAT1) 

0-1-83 wjuxe  na,iooa)     alax_._ . 
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0184  WRITE  (10,1019) 

Oia5 UU9  -FORMA-T  (  «0      XILDC         .EILDC  XLR  EILR', 

1    •  EILD  EILB         XILBR         EILBR', 

. -2_.  • _   XILBC         EILBC) 

0186               WRITE  (10,1000)   XILDC,  EILDC,  XLR,  EILR,  EILD,  EILB,  XILBR, 
1 E4t£R*r-X-i-L-&C-r-E-I-4.-BC-      


0187 

04*8- 

0189 


0190 
01  SI 


0192 


WRITE  (10,1030) 
^030  FORMAT  ( • 0 1  LB •  ) 

WRITE  (10,1000)   ILB 


60  CONTINUE 
RETURN 


END 
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STIFF 

The  subroutine  STIFF  assembles  the  inverse  stiffness  matrix, 
the  flexibility  matrix,  for  simple  spans. 

To  assemble  the  flexibility  matrix  the  following  definitions 
are  first  required. 

1.   6-?   is  the  deflection  for  a  1#  load  applied 
lm  rr 

at  the  i   grid  point  and  on  the  j  l  beam,  i.e.,  at  the 
intersection  of  X.  and  Y  .,   (for  the  L   span).   It  is 
evaluated  at  all  longitudinal  points  for  the  j   beam  of 


th  /a-i-x  »  L 

the  L        span   for  [  ^        I      +1]   <m      ^^      Isp.      Hence, 

p=i 


h- 


there  is  a  vector  of  these  deflection  coefficients  for  each 
longitudinal  beam  stiffness  relationships. 

'  i  .  th 

2.  6  .   is  the  deflection  of  the  r   beam  for  a  1#  load 

applied  at  the  j   beam,  and  the  i   ,  grid  point,  i.e., 

at  the  intersection  of  the  X.  and  the  Y„ . .   It  is  evaluated 

at  all  lateral  points  along  the  i    grid  line  for  l<r<N  . 

Hence,  there  is  a  vector  of  these  deflections  for  each 

lateral  point  on  each  grid  line.   These  reflect  the  lateral 
deck,  primarily,  deflection  relationships. 

3.  The  applied  loads  are  assembled  as  a  vector  of  se- 
quential loads  along  the  grid  lines  with  the  grid  lines 
sequenced  by  their  ascending  order,  e.g., 


102 


Given  the  deflections  at  each  grid  center,  for  the  above  example, 
the  deflection  vector  is 


16 


where  6   =  6  ,  .  ,  x  ,  .   •  _    „-/-;  i\i 
n     (3-1)4-1,  i.e.,  n=(3~l;i 


>sed  of  all  the  deflections  due  to  loads  imposed  at 


6,  is  compos 

points  1,2,3,4,5,9  and  13.  Similarly,  the  virtual  load   imposed 
at  each  grid  center  are  directly  related  to  the  composite 
stiffness-deflection  relationships . 


This  routine  calculates  the  uncoupled  effective   main  beam  unit 
load  deflections  for  a  simple  span, 
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(XSUPL~V  (XSUPL+1    V 


2  2 

(XSUPL+l"XSUPL}  "(Xm"XSUPL)  "XSUPL+l'Xi) 


6{XSUPL+1"XSUPL)EI 


lm 


where  a  unit  load  is  applied  at  i  and  its  effect  measured  at  m, 

•i 
This  routine  thus  calls  STIFFL  to  compute  6  .  ..  the  uncoupled 

unit  load  deflection  of  the  lateral  effective  beam.   These  are 
then  used  to  calculate  the  coupled  unit  load  deflections. 

If  a  continuous  beam,  with  more  than  two  supports,  is  to  be 
evaluated,  this  subroutine  first  estimates  the  deflections  in 
a  simple  span  manner,  as  described  above  and,  then,  adds  the 
unit  load  discontinuity  bending  moment  contribution  to  the 
deflection  for  each  grid,  i.e., 


jj        m        j   +  XSUPL+1-Xm 

lm      lm      EIT  .       L 

10 


M 


DLij 


( 


Xm  XSUPL 


XSUPL+1~XSUPL 


) 


where  <5.   is  used  as  a  replacement  of  the  simple  span  value 
lm  r  c  tr 

6-?  in  estimating  the  final  flexibility  matrix. 

Subroutine  STIFFA  is  then  called  to  assemble  the  flexibility 
matrix,  that  is,  the  inverse  of  the  stiffness  matrix,  and  to 
include  the  coupled  effects. 

The  following  subroutines  are  called: 

STIFFL 
STIFFA 
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f         STIFF    J 


ZERO  XK,  XKI,  AND 
ATEMP  MATRICES 


Form  SK  Matrix 

X   =  X   -  X 

m 

i    ♦  (X 
m  v  SUPL+1 


SK. 


m   _m    SUP(L)    r 
^♦(x^^^.-x.)  [(X 


2  —2 

SUPL+l"XSUPL)  "Xm"(XSUPL+l 


-v2] 


i,m,3 


6(XSUPL+1  XSUPL)EI 


Yes 


I  =ML+1 
o 

I  =1  +IS (L) 
1   o 


1=1 


o 

31 


m=I 


SK.  .=  SK.   .+   XSUPL+1  Xm 

im}     im]      EIT  _.  L 


Lj 


M 


-iij\: 


m   SUPL 


/   DLlA  XSUPL+1  XSUPL/ 


01 ';    !  jusupl+i"xsupl/  'iiDLiJ  V  xsupl+i  xsupl; 


m=m+l 


No 


Figure  23 


STIFF  Program  Flow  Chart 
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1=1+1 


No 


Yes 


Yes 


PRINT  OUT 

PROGRAM 

VARIABLES 


STIFFL 


STIFFA 


IX,IO,SLS, 
B,A,XSUP, 
EI  AND  S] 


f     RETURN    J 


Figure  23.   STIFF  Program  Flow  Chart  (Continued) 
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:ORTRAN  IV  G  LEVEL   21 


STIFF 


DATE 


73170 


00/31/2 


0001 


--C-- 


0002 
0O03 

0004 


SUBROUTINE  STIFF(L) 

INTEGER  P,  Q,  R 
REAL  MD,  MOT 

-£©MWGW-  /-&4K  2V  -DT-, 


1 
_2__ 

3 
-A — 
6 

-7 — 


THRSB+   IE, 

NL,  NSI,  NSUP, 


IEVENT*  IN*  10,-ISS,  ITER, 
STEADY,  T,  WEIGHT,  XLB,  YBCL , 


NAP,  NB,  ND,  NI, 

NXPTS,  NDBUG( 10), 

DELW(IO),  DELX(IO),  EEI(10,10), 

XSUPUO),  Y£(L04,  YP14),  YCL(IO) 

E(3,10),  EI(3,10),  XCS(3,10),  ZBL(3,10),  ZBU(3,10), 

F(80»10)-+-  MD1  11,20,  10),  -MDItll,aOU.  W.P(-50,41,  X(80)  ,XPt50) 


EIL(  10) ,  IS( 10) ,  NCS( 10), 
ZCL(10,10),  ZCU(10,10), 


0005 


0006 
-0007 


.COMMON  /BLK4/  X1U80, 801, XKI  (80,80),  ATEMP(  80,  80  )  ,  DTEMP(80) 
1  ,MI(20),  MB(20),  SIGB(20,10) 


DIMENSION  SK(20,20,10) ,  SKP(10,10) 
EQUIVALENCE  4--SJU  1,1,1  )+AIEMP_i  1  *1_U- 


LSKPCL,-1),ATEMP(1,51)) 


-Z_ER0_.QUJ_BLK.4 


Q.Q08_ 

0009 

0010 


DO  2  1=1,80 
DO  2  J=l,80 
XK(  Tf.l>  =  O.n 


0011 

J1CL12 

0013 


XKI  (I  ,J)  =  0.0 

AXEMPJL,.  J  )__==_  0^0 

2  CONTINUE 

CALCULATE  STIFFNESS  ELEMENTS 


0014 
0015 
0016 
0017 
0018- 


100 


ML 
LL 
DO 
ML 
-ML- 


=  0 
=  L 

100  NN 
=  ML  + 
=l  .ML-  -- 


=  1 ,  LL  _ 
IS(NN) 
XSIL1 


0019 
0020 
0021 
0022 


002  3 


JX 

IY 

JU 

XO 


=  IS(L) 
=  IS(L)/2 

=  XSUP(L+1) 
=  XSUP(L) 


DO  60  J=1,NB 


0024 

D025 

0026 

0027 


002  8 
1)029- 
0030 
0031 


IB  =  1 

I  =  1S(L) 
10  CONTINUE 
M_=  J 

MC  =  IB 
_  B  =  XI  -  X(  I+ML) 
20  A  =  X(M+ML )  -  XO 
AM  =  XIM+ML) 


-&032- 


-N. 
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STIFF 


DATE 


73170 


00/31/2 


0033 
0034 
0035 

0036 

0037 
-0038 

0039 
J)04Q.- 


0041 
004? 


30 


IF(AM.LT.XCS(N, J) )  GO  TO  50 

N  =  N  +  1 

IF(N.LE.NCS( J))  GO  TO  30 

STOP  13 


50  CONTINUE 
._.  SKiI,M,J)  =  (B*A*((X1-X0)**2-A**2-B**2))/(6.*(X1-X0)*EI(N,J) ) 

SK(  IB,MC,J)  =  SK( I,M, J) 
M_  =_M_  -  _1_ 

MC  =  MC  +  1 
LF  (M.OF  *1± 


GO .J.Q  20_ 


0043_ 

0044 

0045 


1  =  1  -_1 
IB  =  IB  +  1 
JF_iI_.GT.IY)   GO  TO  10 


0046 


JJ-__=_2_. 


0047 

0Q48_ 

0049 

00_0„ 

0051 

00  52 


0053 
005  4_ 
0055 
0056 
0057 


I  =  IS(L)  -  1 
55  M  =  IX 

MC  =  1 
57  SKU,MtJ)  =  SK(M,I,J) 

SK(IB,MC,J)  =  SK(MC,IB,J) 
_M_„__M  -_1 

MC  =  MC  +  1 

IF(M.GT.I  )   GO  TO  57 

1  =  1-1 

IB  =  IB  +  1 

IF(  I.GT.IY)   GO  TO  55 


0058 

0059 

Q060_ 

0061 

.0062 

0063 

0064 
006  5 
0066 


0067 
006  8 
0069 
£070 
0071 

0072 
0Q7  3 
0074 


81 

80 

60 

X 


IF  (NSUP.LE.2)  GO  TO  80 

10  =  ML  +  1 

11  =  10  +  IS(L)  -  1 
00  81  1=10,11 

_DD.  fll  ML  =  10.,  II     

SK(I,M,J)  =  SK(I,M,J)+(X(I )-XO)*DELX(L)/( EEI(L,J)*(X1-X0) )* 
1      ( (Xl-X(M)  )*MD(L,I,J)  +  (X(M)-X0)*MD(L+1, I, J) ) 

CONTINUE 

CONTINUE 

CONTINUE 


IF  (NDBUG(3).EQ.O)   GO  TO  87 
WRITE  (10,1010) 

1010  FORMAT  ( "OXSUP  • ) 
WRITE  (10,1011)   XSUP 

1011  FORMAT* 1X,10E13.4) 


WRITE   (10,  1025) 
1025  FORMAT( • OE I • ) 

WRITE  (10,1011)  EI 


0075 
0076 


ISL  = 
WR.ITE 


IS(L) 
1__U1D19) 
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ORTRAN  IV  G  LEVEL   21 


STIFF 


SK«  ) 


DATE 


73170 


00/31/2 


1019  FORMAT( «0    M 

DO. 69  J=1,NB 

DO  70  M=1,ISL 

WR4-7E  (10,1020)  M, J, ( SK( I ,M, J ) , I = 1, I SL ) 

10  20  FORM AT (IX, 2 14, 10E12. 3/ IX ,8X , 10E12.3) 

70  CONTINUE — 


0083 
_00B4 

0085 
J30&6 


69 

CONTINUE 

-87. 

CONTINUE 

CALL  STIFFL 

(L) 

CALL  STIEEA. 

CD 

0087 
0088 


RETURN 
END 
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STIFFL 

The  subroutine  calculates  the  lateral  unit  load  deflections 
6'  as 


2   „  2    ,„     „   ,2 


(Y    -Y  )     Y   (Y     -  Y    -  (Y    -  Y  )     ) 
•      v  SUP   Bj;   Bj v  SUP     Bj    v  SUP    Bj'  ; 


jr  6YSUPEIL 

No  subroutines  are  called  by  this  routine. 
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STIFFL 


4/ 


FORM  SKP  MATRIX 


SKPj,r  + 


(YSUP-YB.)  (YB.)  (YSUP2-YB.2-(YSUP-YB.)2) 
6  (YSUP)EIL 


PRINT  OUT 

PROGRAM 

VARIABLES 


JL 


RETURN 


YSUP,SKP 


Figure   24«   STIFFL  Program  Flow  Chart 
111 


:ORTRAN  IV  G  LEVEL   21 


STIFFL 


DATE 


73170 


00/31/2 


0001 


SUBROUTINE  STIFFL  (L) 


0002 


REAL  MD,  MDT 


000  3 


IN,  10,  ISS,  ITER, 
T,_  WEIGHT,  XLBf  YBCL, 


COMMON  /BLK2/  DT ,       THRSH,   IE,  IEVENT, 
1 HAP^- W6-,- NO-,  -N4-^-NL r  -NSI+ -NSUP ,-.  STEADY* 

2         NXPTS,    NDBUG(IO), 
3-     -   DELW<1-Q),DELX(10),     EEIUO.IO),     EIL(IO),     IS(10),    NCS(IO), 

4         XSUP(IO),    YB(10),    YP(4),    YCL(IO),    ZCL(10,10),     ZCU(10,10), 
6         E(3,10),    EI13,1Q),   J(CS(3,10),    ZBL(3,10)  ,     ZBU(3,10)  , 

7         F(80,10),    MD(  11,20,  10),    MDT(11,10),    WP(50,4),    X(80),XP<50) 


0004 


COMMON    /BLK4/    XK ( 80,80 ) ,XKI ( 80,80) ,    ATEMP( 80, 80 ) , 
._ _        _    ,ML(20),    MB(20)»    SIGB(20,10) 


0TEMP(80) 


0005- 
0006 


0007 

oooa- 

0009 
HQ10 


_D_IM£NSIQR   SK(20,2D,J.Q)  ,     SKP(  10,10) 
EQUIVALENCE    ( SK( 1 , 1 , 1 ) , ATEMP ( 1 , 1 )  ) , 


1010 
___77_ 


IF     (NDBUG(3).EQ.O)       GO    TO    77 

-WRixE  uoaoioj      _ 

FORMAT    («0  YSUP         J 

CONTINUE 


(SKP( 1, 1), ATEMPf 1,51) ) 


B       JR 


SKPM 


0011 


0012 
iKLL3„ 


YSUP 


?*YBCL 


_&0 


J    =    NB 
.COMTJNUE 


0£LL4^ 

0015 

0016 


a_.=  _-YSUE 

JR    =    J 
QO    r.niVITTNIIF 


Y  B  (  J  ) 


JHLL7__ 

0018 

_0_019_ 


.A__=L  YB(JR) 
SKP(J,JR) 
SKP(JR,J) 


=  (B*A*(YSUP**2-A**2-B**2) )/( 6. *YSUP*E IL ( L ) ) 
=  SKP(J,JR) 


-OO20- 


0021 
0022 
0023 

0024 
0025 


0026 
0027- 


1012 

87 

_C 


IF  (-NDB1LG(3J.^Q.0L  -  GO  _I0  _8J 

WRITE  (10,1012)  YSUP,  J,  B,  JR ,  A,  SKP(J,JR) 

FORMAT  (1X,E12.3,  14,  E 12. 3, I 4,2E12. 3) 

CONTINUE 


=  JR  -  1 
--UWUN1  .-OJ- 
J  =  J  -  1 
IF  (J.NE.O) 


JR 


GQ-TO  90 


GO  TO  80 


-0Q2& 
0029 


RETURN 
END 
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STIFFA 


This  subroutine  includes  the  coupled  effect,  as  described  below, 
and  assembles  the  flexibility  matrix.   The  coupling  is  accomplished 
at  the  same  time  as  the  assembly. 

For  the  main  diagonal  deflection  coefficients,  that  is,  at  the 
point  of  coupling,  it  can  be  shown  that  the  total  deflection  at 
the  point  of  coupling,  6  j .  ,  is 

6: .  6. . 

x'       -      1:L   ia- 

hi      Tr-^— 

where  6'.    is  the  deflection  of  the  diaphram  at  ii 

11  # 
due  to  a  1   load  at  ii  and  <s .  .  is  the  deflection  of  the  effective 

#  ia" 
main  beam  at  ii  due  to  a  1   load  at  ii. 

The  proof  is  as  follows: 

given: 

f !  . 
ii 


«u=51Ti-  '  A<  "  c'u   c'  "  fii  sii 
«ii=eM—  A2  "  £ii  c  "  fii  6ii 


and 


then 


f ! .  +  f . .='1# 
ii     ii 


f !  .  =  1  -  f .  . 
ii         ii 


substituting 
and 


«ii  =    (1  -  fii)6ii 


£ii«ii=  (1  -  fii)6£i 

ii    ii        ii  ii    ii 

81. 
f 


ii  6i±+Sii 
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Substituting  for  f.. 


6ii5ii 


11   6ii+6ii 


To  analyze  the  physical  implications  we  change  variables  to  the 
stiffness  coefficients,  as  if  each  beam  were  uncoupled,  that  is, 

n     1/k' . .  1/k. . 

1      '   11  '  11 


Simplifying: 


lc      1/k..+  1/k! 

ii 


k-  •  =  k. .  +  kJ . 
^n    ii    ii 


This  reflects  a  parallel  stiffness  effect,  analogous  to  two 
parallel  springs  of  different  strengths,  which  is  exactly  the 
case  at  the  point  of  coupling. 

Hence,  the  total  unit  load  deflection  at  the  point  of  coupling 
can  be  determined  from  the  two  uncoupled  unit  load  deflections. 

For  a  far  coupled  point  in  the  lateral  direction,  the 
deflection  at  a  point  ij  due  to  a  unit  load  at  ii  can  be 
determined  by  extending  the  above  analysis.   Since  f '  is  that 
fraction  of  the  unit  load  being  applied  to  the  diaphram, 
then 

S.   =  f  !  .  6  !  . 
ID     ii  ID 

where  6 '  .  is  the  uncoupled  deflection  at  ij  due  to  a  unit  load 
at  ii,  and  6.  .  is  the  resulting  deflection  at  ij  due  to  a  coupled 
effect  at  ii  and  a  unit  load  applied  at  ii. 

Since 

6ii 
f »  .  = 


n    6.  .  +  6.". 
un   ii 
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5  .'  .  6. .  . 
(5-  •=  — ^ /  and  again 

3         6. .+6! . 
n   n 

the  coupled  deflection  can  be  determined  from  the  unit  load 

uncoupled  deflections. 

To  analyze  the  physical  implications  we  substitute  the  uncoupled 
stiffness  coefficients, 

,      1/k!  .  1/k.  . 

1    _       IT      11 

~~=      1/k-  +l/k!. 
k .  .       n    ii 
■*"  J 

Simplifying 

k!  . 

k.  .  =  ,43 —  (k!  ,+k.  .)  . 
in    k! .      ii   n 
J     ii 

This  implies  again,  a  parallel  effect  but  modified  by  a  proportion- 
ality factor  determined  from  the  relationship  of  the  uncoupled 
stiffness  coefficients  at  the  point  of  coupling  and  at  the  far 
point. 
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B 


r=r+l 


STIFFA 


k. 


N=NB.IS. 


\L 


ZERO  K 
ARRAY 


k. 


j  -  i 


i 


>/ 


k_ 


p=(j-l)IS1+i 


k- 


q  =  (r-l)IS1+i 


61      6.  .  . 

6pq  "  5r.+   6.. . 
**      ID    ii] 


Figure  25.    STIFFA  Program  Flow  Chart 
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Figure   25.    STIFFA  Program  Flow  Chart  (Continued) 
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PRINT  LATERAL  AND 
LONGITUDINAL  UNIT 
LOAD  DEFLECTION 


PRINT 

FLEXIBILITY 

MATRIX 


± 


m  =  [5] -^5-  +  [i; 

At 


A 
71 


5,  6' 


y 


Figure  25.   STIFFA  Program  Flow  Chart  (Continued) 
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Figure  25.   STIFFA  Program  Flow  Chart  (Continued) 
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:ORTRAN  IV  G  LEVEL   21 


STIFFA 


DATE 


73170 


00/31/2" 


0001 


SUBROUTINE  STIFFA  (L) 
ASSEMBLE  STIFFNESS  MATRIX 


0002 
0003 


INTEGER    Pt    Q,    R 
REAL-NO-,-  MOT 


0004- 


_CQMM©N-/BtKl/    XKBAR(80,8Q) 


&QQ5- 


-CQMMQW  /BL-K2-/    DT^- 


THRSH, .  -  IE,    IEVENT,     IN,    10,     ISS,    ITER, 


1         NAP,    NB,    ND,    NI,    NL,    NSI,    NSUP,     STEADY,     T,    WEIGHT,    XLB,    YBCL , 
-2 NXPTS,    NDBLIG1  10)  , . _. _  .  


-moi>- 


3         DELW(IO),    DELX(IO),    EEI(10,10),     EIL(IO),     IS(10),    NCS(IO), 
_4 JCSUP-UQl,    YBILO),    YP(4),    YCL(IO)*    ZCL(10,10),    ZCU(  10,10), 

6         E(3,10),    EI(3,10),    XCS(3,10),    ZBL(3,10),    ZBU(3,iO), 
.J E4-8Jl*JLQ-)-i  .MD( 11,20,10),    MDT( 11,10),    WP(50,4),    X(80),XP(50) 

COMMON    /-B4JC4/    XK<  80^8X)-)rXIUX80,43ai^-ATEMP4aiV8Q)-,    DTEMP180) 
1  ,MI(20),    MB(20),    SIGB(20,10) 


0007 

aao&_ 


DIMENSION    SK(20,20,10),    SKP(10,10) 

£OU_I  VALENCE    1  SK  (1,1,1  ),ATEMP(  1,1)),    (SKP(  1, 1 )  ,  ATEM?<  1 ,51  )  ) 


0009 


_N_ 


nr*t<;(i  ) . 


_Q_Q_10_. 
0011 
001? 

nnn 


DQ_JLQ1)_JP =1  iN  „ 

DO    100    Q=1,N 
J.00_JUUP,_Q.)    .=    0.0_ 

nn  400  j=i.nb 


0014 
0Q15 
0016 

0017 
-04118- 


_C_ 


0019 
0020 
0021 

0022 
-OA2-3- 


0024 
0025 

0O26 
0027 


ISL  =  IS(L) 

DO  400  1=1, ISL 

P  =  (J-1)*IS(L)  +  I 

DO  200  R=1,NB 

Q=M-R=-l)*LSiL) .-+-1 ■  -  -  ■ 

IF  (P.EO.Q)   GO  TO  200 

XKJ(P,Q)  =   (SKU  ,  I,J)*SKP(  J,R))/(SK<  I, I, J)  +  SKP(J,J)) 
200  CONTINUE 

DO  320  M=1,ISL 

0--  =— ( J=  U-*lS-(  L  i_±--H_ _.  

IF  IP.EQ.Q)   GO  TO  300 


XKI(P,Q)  =   SK(I,M,J) 

(SK(IfltJ)  *  SKP  ( J, J)  )/( SK(  I,  I, J)  +  SKP(J,J)> 


_...  GO  TO   320 
300  XKI (P,Q)  = 


0028 
4H129 

0030 

-O03-1 


320  CONTINUE 
400  CONTINUE 

IF(NDBUG(3).NE.l)   GO  TO  450 

WJIUE  4-10,1400)  ISL* ISL, NB 
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STIFFA 


DATE  =  73170 


00/31/2 


0032 


1400  FORMAT  ( ■ 1SK  ' ,316) 


0033 

00  34 

0035 

-0036 


0037 
00  3A 


DO  420  1=1, NB 
DO  410  J=1,ISL 
410  WRITE  ( 10,1320) 
WRITE-  (  IQ.,.1 4-014 


1401  FORMAT  (  'O'  ) 
420  CONTINUE 


( SK(K, J, I ) ,K=1,ISL) 


003-9 
0040 
0041. 


WRITE  (  10,1405)  _NB,  NB 

1405  FORMAT  ( 'OSKP',316) 
. DH^43IL  ^1=J._,.N.B 


0042 

0043 

0044 
0045 


430  WRITE  (10,1320)   (SKP(K,J),   K=1,NB) 
-C- ._ 

450  CONTINUE 

_C 

MAXIT  =  3 
LERR.._=_JQ .  . 


_0Q46_ 
0047 
004  8  _ 


; WRITE  OUT  THE  STIFFNESS  MATRIX  FOR  EACH  SPAN 

IF  (NDBUG(3) .NE.l)   GO  TO  350 

WRITE( 10,1139) 
1139  FORMAT(  '1'  ) 


0049 
0050 
0051 
0052 
0053 


_   WRITE( 10,1139) 
DO  340  P=1,N 
WRITE  (10,1140)  P 
1140  FORMAT  ( IX,  'XKI ' , 14) 
WRIJE_  (  I0_,1320)  UKIJP.Q)  ,Q=1,N) 


0054 
0055 


0056 


0057 
0058 
0059 
0060 
JD0A1 
0062 
0063 


0064 
0065 
0066 
006  7 
0068 
-0069- 


1320  FORMAT  (1X,10E13.4) 
340  CONTINUE 


350  CONTINUE 


800  CONTINUE 

DO  850  1=1, N 
._  DO  850  J=1,N 

XKBARUtJ)  =  XKI(I,J)  *  ((4.  #  DELW(  L )  )7DT**2) 

IF(I.EQ.J)   XKBAR(I.J)  =  XKBAR(I,J)  +  1.0 

850  CONTINUE 

CALL  MATINV  ( 80  ,N ,XKBAR ,XK , ATEMP ,  DTEMP,  MAXIT,  IERR) 


WRITE  OUT  THE  STIFFNESS  MATRIX  MODIFIED 


IF  (NDBUG(3  ).NE.l)   GO  TO  860 
_ WRITE( 10,1139) 

DO  830  P=1,N 

WRITE (  10,1830)     P 
1830    FORMAT! IX,  •  XK  '  ,  14, •       MODIFIED') 
.WRIXE  JJOaiLO)        (XK(P,Q)  ,Q  =  1,N) 
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0070 
00-71- 


830    CONTINUE 
860- CONT  INUE— . 


^07-2 

0073 


-R-E-TDRN 
END 
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The  basis  of  the  continuous  main  beam  (girder)  solution  util- 
ized by  this  program  is  the  division  of  a  bridge  into  free  bodies 
that  appear  to  be  simple  spans.   In  order  to  reconstitute  the 
bridge  as  one  body  continuity  must  be  re-established  after  gener- 
ating the  simple  span  stresses.   This  is  accomplished  by  develop- 
ing a  set  of  unit  load  discontinuity  moment  vectors  in  the 
initialization  portion  of  the  program  and  operating  in  these 
vectors  during  the  dynamic  calculations  with  the  actual  live 
load  to  determine  the  actual  discontinuity  moments.   The  contrib- 
ution of  the  actual  discontinuity  moments  to  the  individual  grid 
(element)  dynamic  stress  is  then  calculated  and  added  to  the 
simple  span  stress.   By  using  this  approach,  i.e.,  initialization 
of  unit  load  discontinuity  moment  vector,  calculational  time  is 
minimized.   This  is  especially  significant  because  of  the  inversion 
of  a  matrix  in  the  calculations. 

This  subroutine  calculates,  in  the  initialization  portion  of  the 
program,  the  unit  load  discontinuity  vector  for  each  grid  (element) 
as  if  a  unit  load  were  applied  on  that  grid.   It  further  generates 
a  set  of  these  vectors  for  each  effective  beam,  i.e., 

where  A  and  IB  were  defined  in  the  basic  reports. 

Since  it  must  calculate  the  value  of  EI  for  each  effective  beam 
and  span  in  order  to  calculate  the  unit  load  discontinuity  mom- 
ents, it  also 

1.   Interpolates  for  EI  as  a  function  of  position  of 
change  in  the  effective  beam's  cross  section  for  each 
span,  e.g.,  in  the  case  of  a  bottom  lever. 
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2.   Calculates  the  cross  sectional  moduli,  ZCU,  and  ZCU,  , 
for  each  span  and  beam. 

No  subroutines  are  called  by  this  routine. 
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CLEAR  AA,  SA, 
SB,  A,  AINV 
MATRICES 


X  = 
o 

I  = 
o 

i  = 

L  = 


A 


XL  =  XSUPL+1-XSUPL 


H 


-  [(Xi-XSUPL)/XL](XSUPL+l-Xi) 


Al  =  H/2  (Xi-XSUPL) 
A2  =  H/2  (XsupL+1-X.) 

AA.L=  A1+A2 

Ha  =  Hb  =  H/2 

XA1=  (X.+XSUPT )/2-XSUPL 
X       L 

XB1=  (X.+XSUPT^, )/2-XSUPL 
l      L+l  ' 


a   = 


HaXL     1 

XA1-XL   Hb/XBl-Ha/(XAl_xL) 


b   =  XL-a 
bAA 


lL 


Bl1"  Vl 


B 


Figure  26.   DMOMNT  Program  Flow  Chart 
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PRINT  OUT 

PROGRAM 

PARAMETERS 


i=i+l 


VXSUPL+1 


±L 


L=^L+1 


7 


ABAR,    CBAR 
DBAR 


Yes 


-X 


^^T) 


Figure  26.    DMOMNT  Program  Flow  Chart  (Continued) 
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V 


K> 


J  "  1 


N  =  1 

X  =  0 
o 

k  =  0 

1  =  0 
o 

L  =  1 


_±_ 


X  =XSUP  , . 
o      L+l 


Yes 


±C 


EEI=EI 

L    n,j 

ZCU  =ZBU 
L     n,3 

ZCL  =ZBL 
L     n,j 


No 


Figure   26.    DMOMNT  Program  Flow  Chart  (Continued) 

127 


EEI,    = 


ZCU.    = 


ZCLT    = 


XCS       .-XSUPT 

n,j  L         EI 

XSUPTJ,-XSUP_ 

Jj+ J.  Li 


n,D 


XCS       .-XSUPn 
n,j  1 

XSUPL+1-XSUPL      n,j 


ZBU 


XI 


XCS   .-XSUPT 
n'3      L 

XSUP_  ,n-XSUPT 
LtL       L 


XCS   .-XSUPT 


ZBL 


n,D 


Yes 


H 


EEI, 


ZCU, 


XSUPT,,-XCS 
EEIT  +  ^±i *LJ. 


"L    XSUP  +,-XSUP, 
XSUPT.,-XCS 

zcu,  +  5d^ !ia 


EI 


n+l,j 


ZCL    =  ZCL„  + 


XI 


BL, 


=  I 


XSUPT.,-XCS 

L+l    _  n / 3 

L    XSUP   ,-XSUP 

XSUPTJ,-XCS 
BL_    XSUPT^..-XSUPT 

Li  L+L         Li 


ZBU 


n+1,  j 


ZBL 


n+l,j 


"BL 


n,D 


Figure  26.   DMOMNT  Program  Flow  Chart  (Continued) 
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EEIL  =EEI 


XCS  L.   .-XCS 
+     n+1'3     n>l  EI 

XSUPL+1-XSUP        n+l,j 

XSUPr^,-XCS  ,.   . 
L+l     n+1,3 

xsupl+1-xsupl 


ZCU 


.  -XCS 

=zcu  +  n+1'l 

^^UL      XSUPT  , 

L-rl 


XCS  ,_   .-XCS 

n+l,j    .11,3 

XSUPT  , -XSUPT 


EI  ^0  . 
n+2,3 

ZBU  .,  . 
n+1,3 


ZCL 


L 

XSUPTJ,-  XCS  ..   . 
+      L+l      n+1,3   7Rn 

XSUPT^-,-XSUP.        n+2,j 
L+l      ii  J 

XCS  ,,  .-XCS 
=  7n      +  n+l,,j     n,] 

L      XSUP  ,-XSUP 
Jj  +1      L 


ZBL  , .   . 
n+1,3 

ZBL  ^0  . 
n+2,3 


XIBLT 


xsuPL+1-xcsn>. 
xsupl+1-xsupl 

xcs   ,.    .-x^0 

n+1,3   CSn,3   I 
XSUPL+l"XSUPL     BLn+lj 


=I   +  n+1,3 

BLL    XSUPL+1  XSUPL 


-X 
-X 


,    XSUPL+rXCXn,j         z 

XSUPL+l"XSUPL     BLn+2j 


r 


n=n+2 


n=n+l 


No 


L  -L+l 


7 


X  =XSUP_ 
o      u 


Figure  2-6.     DMOMNT  Program  Flow  Chart  (Continued) 
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COMPUTE  [A] .  DIAGONAL  ELEMENTS 
3 


xsupl+1-  xsupl 


'1,1-1 


XIBL, 


,  2(^L±2r 

1,1  *\         XIBL 


+1-XSUPL      +   ^SUPL+2"XSUPL 


XSUPL+2   "    XSUP 


XIBL, 


) 


1,1+1 


L+l 


XIBL, 


Figure   26.   DMOMNT  Program  Flow  Chart  (Continued) 
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[MD]k    =    -6[A]T1[B],      . 
3  J  k/D 


X. 


k=k+l 


yes 


Yes 


M 


PRINT    OUT 

PROGRAM 

PARAMETERS 


XCS,EI,EEI, 
ZBL,ZCL, 
'ZBU,ZCU,A, 
A-MD 


Figure  26.    DMOMNT  Program  Flow  Chart  (Continued) 
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ORTRAN  IV  G  LEVEL   21 


OMOMNT 


DATE 


73170 


00/31/27 


0001 


0002 

00O3- 

0004 


0005 


0006 


0007 


SUBROUTINE  OMOMNT 

INTEGER  P 

REAL-IB4  

REAL  MD,  MOT 


IN,  10,  ISS,  ITER, 
T,  WEIGHT,  XLB,  YBCL, 


COMMON  /BLK2/  OT ,       THRSH,   IE,  IEVENT, 
-1 NAP,  N&,  ND,  Ml,  NL,  NSI,  NSUP,  STEADY, 

2    NXPTS,  NDBUG(IO), 
-3 DELW-ILQ1,  DELXUO),  EEI(10,10),  EIL(IO),  IS(10),  NCS(IO), 

4    XSUP(IO),  YB(10),  YP(4),  YCL(IO),  ZCL(10,10),  ZCU',10,10), 
_6 E(3,10),  Ella,  10  )f-XCS  (3+ -L0)^_ZBU-3^  10),  ZBU(3,  10), 

7    F(80,10),  MD(  11,20,  10),  MDT(11,10),  WP(50,4),  X(80),XP(50) 


X 


COMMON    /BLK6/     IBL(3,10) 

DIMENSION    XIBL(10,10),    AA(20,10),     SA( 20, 10 ) , SB ( 20, 10 ) , 
4 A4r9^9  )  >     B  (-9-)-,- A I-NV -I  9^94 —        -- 


0O08- 


0009 

-ooio- 

0011 
0012 


0013 
0014- 


0015 

0016 
0017 


0018 
0019 


-DO    5-5-L-lrlO 
DO    5    1=1,20 

AAIUL-J 

SA(I,L) 
SB-Cl^ 


DO    54    1=1,10 
.54-  X1BL1 UL-l   .=-  Q.  0 
55    CONTINUE 


DO 
-OO 


J  =  l,9 

A(I,J)    =   0.0 
AINVU  ,J)    =    0.0 


0020- 
0021 


J0_  = 
I    = 


0022 
0023 
0024 

0025 
0026- 


0027 
0O28 
0029 
0030 
0031 
0032 


IF  (NDBUG(4).EQ.0)  GO  TO  15 

WRIJE  (10,1005) 

1005  FORMAT  ( 'ODMOMNT • , /// , •     I 
—   1--   •     Al  A2 

15  CONTINUE 

Hr-  ^-0 

DO  20  L=1,NSI 
=  XSUP(L) 
=  XSUP(L+1) 
=  11  +  1 

II    =    10   +    IS(L)-1 
ZO-I^+GtII 


XO 
XI 
10 

-DO 


XO 


AA 


SA 


SB') 


0033 
0034 
0035 
0036 
0037- 


16  CONTINUE  - 

H  =  (X( I )-X0)/(Xl-X0)*(Xl-X( I ) ) 

A2  =  H*(X1-X( I ) )*.5 

Al  =  (H/2. )*(X(I )-X0) 
AA-<-trL)  -=-Al     +  A2- 
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0038 
0039 
0040 
0044 
0042 
0043 


0044 

0045 

nntA 
uu  t  o 

0047 

004  8 


XA1  =  (X( 
XB1  =  (X( 

HA  =  H/2. 
HB  =  HA 
XL  =  Xl-X 


I  ) 

I  ) 


0 


+  X0)/2.-  X0 
+  XD/2.  -  XO 


--G 


-SA  ( I  ,L  )— ^4*L*X-B14Vf  XL~X-A1+*B1  I 
SB(I,L)  =  XL  -  SA( I,L) 


IF  (NDBUG 
-WIVRE  (  10 


1018  FORMAT  (1 
18  GONTINUi- 


(4).EQ.0)  GO  TO  1 
,1018)  L,  It  X(I ) 
X,  214,  2F10.3,  5 


,  XO,  Al,  A2, 
E13.4) 


AA(  I,L)  ,  SA(  I,L)  ,  SB(  I,L) 


4)049- 


2Q--&QNT  I NUE 


0050- 


0051 
0052 


0053 


0054 
-00545- 


0056 


0057 


-00  500  J=1,NB 
N  =  1 


10  =  0 


DO  150  L= 
-XO  -=    XSUJ? 


1,NSI 
(L+L) 
XI  =  XSUP(L) 


100  IF(XO.LE.XCS(N, J) )   GO  TO  120 


0058 
-00-59- 
0060 
0061 


XS  =  XSUP 
-EE44-L,J4 
ZCL(L,J) 
ZCLHL,..!) 


0062 
0063- 


0064 


C-- 


0065 
0066 


XIBL(L,J) 

-XS_= 

IFtXO.GT. 

EEI  (L,J) 

ZCLIL-tJ4- 


(L) 

-  IXCS(N,J)  -  XS) 
=  (XCS(NtJ)  -  XS) 
=  (XCS(N,J)  _-_XSl 

=  (XCS(NtJ)  -  XS 

XSUPIL+1) 
XCS(N+1, J) )  GO  TO  110 


/(XO  -  XS)*EI(N,J) 
/(XO  -  XS)*ZBL(N,J) 
/  (xa  _-  X  SJ  * Z BU  { M,  J 1 
)/(XO-XS)*IBL(N, J) 


0067 
006  8 
0069 
0070 

-00-74- 


ZCU(L,J) 
XIBL(L,J) 
N    =    N    +     1 
GO    TO    150 


=  EEKL,  J)  +  (XS 
~_ZCL44  »4J  +  IXS 
=  ZCU(L,J)  +  (XS 
=    XIBL(L.J)    +    (X 


-  XCS(N, J) )/( XO    -    X1)*EI(N+1,J) 
=_.XCS(N,4)4/(X0    -    X1J*ZBL(N+1,J) 

-  XCS(N,J))/(XO    -    X1)*ZBU(N+1,J) 
S-XCS(N, J) )/(XO-Xl)*IBL( N+1,J) 


443—EE44  W44-  ^ 


0072 


0033- 


-0O7-4- 


XJ07-5-  ■- 

0076 

0O77 

0078 

0079 


-120 


1  +     (XS 

ZCL(L,J) 

1         +    (XS 

---Z-CUUt-JI- 

1         +     (XS 

-  XIBL(L,J4 

1  +    (X 

..  ..N-=  M   *-2 

GO    TO    150 

EEI  (L,J> 

ZCU(L,J) 

_Z-C4.-44.-rJ4- 


-EE-I  (LrJJ  -+-4XCS 

-  XCS(N+1,J) )/(X0 
=    ZCL(L,J)     +     (XCS 

-  XCSIN+1, J) )/(X0 
=    ZCU(L,J)     +     (XCS 

-  XCS(N+1, J) )/(X0 
_=_X  IB4  XL-  ,-4)  *  4-X 
S-XCS(N+1, J) )/(X0 


=    EI  IN,  J) 
=    ZBU(N,J) 
=_iat4M-MU— 


(N+1,J)-.-    XCSI.N,  J)  )/(X0    -    XI)     *EI(N  +  1,J) 

-  Xl)*EI(N+2fJ) 
(N+1,J)     -    XCS(N, J) )/( XO    -    X1)*ZBL(N+1 ,J ) 

-  Xl)*ZBL(N+2, J) 
(N+1,J)     -    XCStN, J) )/(X0    -    X1)*ZBU(N+1 ,J) 

-  Xl)*ZBU(N+2, J) 
CS  (41+  1 ,  J4-XC  S(N,J))/(X0-X1)*IBL(N+1,J) 
-X1)*IBL(N+2,J) 
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0080 
0081 
0082 


0083 
GQ84- 


0086 
0087 


0088 
-0O89- 
0090 
0091 


XIBL(L»J)     =    IBL(N,J) 


150    CONTINUE 


NSIM1 


NSI    -    1 


0085 


DO    200 

DO  2QQ 

200    A(LtM) 


L=1,NSIM1 
M=4-,NSIM1 
=   0.0 


DO    210    L=ltNSIMl 
XO    =    XSUP(LJ 


XI  =  XSUP(L+1) 
_JC2_a5L_XSLLRiL-t2) 

IF(L.EQ.l)  GO 
4U-UI-J.J--    iXl- 


TO    205 

-    XO)/XIBL(L,  J) 


0092 


009  3 
D094 


009  5 


0096 


0097 
0_Q98_ 


205    A(LiL)    =    2-*qiX4-  ^-XO^X-I-B4.-l-W-J4- 

IF  (L.EQ.NSIM1)       GO    TO    210 
A(l.,L+JJ-=— UCZ,-  -XJJ-ZXTBLt  L+ 1 , J ) 


+-4X2--    XUVXIBUL+UJU 


210    CONTINUE 


CALL    MATINV     (9,NSIM1,    A,    AINV,     ADUM,    AD,    3,     IERR) 


11=0 

nn    301     IN: 


l,NSiMX 


0099 
01O0- 
0101 
010? 


10=11+1 
JL1  =  10    +   XSXLM1  - 
DO    302    K=I0,I1 

on  300  i  =:u.niimj. 


0103 
QlOA 
0105 
0106 

0107 


300 


IF  (XSUP(L+1 )  .LE.X(K) )   GO  TO  300 
IF  (XSUP(L).GT.X(K) )   GO  TO  300 

B(L)=AA(K,L)*SB(K,L)/( ( XSUP ( L+ 1 ) -XSUP( L ) ) *X IBL ( L, J ) ) 
CONTINUE 


0108 

0109 
0-1  ID 
0111 
0112 


0113 
QUA 

0-115- 

0116 

D44-7- 


00-3  54.  -U,   =_1+NSIM-1 

351    MD(LL+1,K, J)=0. 

DO    350    I=1,NSIM1 
350    MD(LL+1,K,J)    =    MD(LL+1,K,J)    -    6.    *    AINV    (LL,I)*B(I) 

MD(1,K,J)=0. 

MD4NSUP-»-K,,-U->  ~=_Q». —  -    -      

302    CONTINUE 
301    CONTINUE 


0118 
OILS 
0120 
0121 
0122 
XU£3- 


IF  (NDBUG(4).EQ.0)  GO  TO  160 

WRITE  (10,1150)  J 
1150  BQBMAX— (-LQB£AM-=.L»IA4 

WRITE  (10,1151) 
10.51  FORMAT  HO  XCS«  ) 

WRITE  (  10,1000)  XCS 

WRITE    (10,1152) 
1152    FORMAT    ( '0    EI     • ) 
WRITE-(ia,1000l   -EI-  134 
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0124  WRITE  (10,1153) 

0125  1153  FORMAT  («0  EEI') 

0126  WRITE  (10,1000)  EEI 

0127 WRITE  (10,1154) 

0128  1154  FORMAT  (»0  ZBL1) 

-0L2-9 WRITE  UO-t-10-QiU  -  XBJ 

0130  WRITE    (10,1155) 

0131-  1153  FORMAT    CO    ZCL  '  ) 

0132  WRITE    (10,1000)    ZCL 

0433 WHITE    (  I0^1L56J 

0134  1156    FORMAT    («0    ZBU«) 

0135 WRTTF     (TO,  1000)     7BU 

0136  WRITE    (10,1157) 

0137 U5JL.EQRMAT    (  •  0    ZCUi  ) 

0138  WRITE    (10,1000)    ZCU 

0139 WRITE-  (  10 ,  2020  ) 

0140  2020    FORMAT( 'OXIBL' ) 

0141 WRIT  F    (IQtlQQll M  ftj. 

0142  WRITE    (10,1158) 

0143 1158    FORMAT    (  '0    A       •) 

0144  WRITE    (10,1001)    A 

0145 WRI1E__(  10,1159) 

0146  1159    FORMAT    ('0    AINV) 

0147 kimE_iJLGLtl001_L  A  imt 

0148  1001     FORMAT     (1X,9E14.4) 

01_49 1000    FORMAT    (1X,10E13.3) 

0150  WRITE     (10,1060) 

0JL51 _    _    WRITE     (10,1001)       B 

0152  1060    FORMAT! • OB • ) 


0153  WRITE  (10,1160) 

0154 ilj>0  FORMAT  («0  MD '  ) 

0155  DO  170  K=l,20 

015j6 WRITE  (10,1170)  (MD(  LL,  K,  J  )  ,  LL=  1 ,  1 1 ) 

0157  1170  FORMAT  (1X,11E11.3) 

0155 17CL XONT  INUE  .  .  - 

C 
0159 .   160  CONTINUE 

C 
0160 500  CONTINUE 

C 

0161 R-EXURN  _._. 

0162  END 
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PNTLOD 

The  subroutine  PNTLOD  calculates  the  vehicle  tracks  across 
the  bridge: 


Yi 


=   Y 


AYr 


CL1 


Lane  1 


Y      +   AYT 
CL1      -ai 


=   Y 


CL2 


-   AY 


T 


Lane  2 


Y      +   AYT 
CL2      -=± 


Axle  loads  are  distributed  laterally  to  the  tracks  to  which 
they  apply. 


Axle  in  Lane  1 


w.    , 
1/1 

= 

w.  . 
11 

2 

w.    0 
1,2 

= 

W.  . 

n 

2 

W.     ., 
1,3 

= 

W.  . 
n 

2 

w.     . 

W. 
li 

Axle  in  Lane  2 


i,4 


where  i   is  the  load  index  (only  axles  actually  on  the  bridge) 
and  ii  is  the  input  axle  load  index  and  where  i  <  ii. 

No  subroutines  are  called  by  this  routine. 
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Y'i=  YCLry 


y-2=  ycl1+ 


T 


Z 

Y,3+  ycVy 


Y'4+  ycVy 


k. 


ZERO   W    ARRAY 
i    =    1 

XMI    =    -1 

ii      =      1 


Figure      27.         PNTLOD   Program  Flow  Chart 
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2slZ- 


Wi-1,3=  Wii/2 


Wi-1,4=   Wii/2 


V 


W! 

1 

,1= 

=  0 

1 

,2  = 

=  0 

i 

,3= 

=  0 

W! 

l 

,4= 

=  0 

-tr 


i=i+l 


ii=ii+l 

XMI=X. 

n-1 


NAP=i-l 


Yes 


f       EXIT      J 


No 


Figure   27.    PNTLOD  Program  Flow  Chart  (Continued) 
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PNTLOD 


DATE  =  73170 


00/31/2 


0001 


C_ 

c 

-C- 

c 


-0002- 


SUBROUTINE  PNTLOD 

THIS  SUBROUTINE  ESTABLISHES  THE  LOCATION  AND  WEIGHT  OF  THE  POINT 
-LOADS  REPRESENTED  BY  THE  INCOMING  DYNAMIC  LOAD  DATA 

REAL  MP,  -MOT _ 


C003 


IEVENT, 
STEADY, 


IN,  10,  ISS,  ITER, 

T,  WEIGHT,  XLB,  YBCL , 


-COMMON  /BLK2/  DT ,-       THRSH,   IE, 
1    NAP,  NB,  ND,  NI,  NL,  NSI,  NSUP, 

2 WLBJS*  -NOBUG-llOJ,   __......_.  . 

3    DELW(IO),  DELX(IO),  EEI(10,10),  EIL(IO),  IS(10),  NCS(IO), 

4 XSURllOU  -JY-RULOJ-,-  \RXJt\  ,_YCLX10_U.  ILL  LIO,  10),    ZC  U  1.1.0  dfl. ) » 

6         E(3,10),    EI(3,10),    XCS(3,10),    ZBL(3,10),    ZBU(3,10), 

2 FiaO,lQ),_MD.<  11,20,10),    MDT(11,10),    WP(50,4),    X(80),XP(50) 


0004^ 


000  5 


0006 
J&0Q7- 


000  8 
0009 


—COMMON  V_BLK3/    SDMHRX,    DTL,     IEVNTX,    NA,  NTRUKX,    LAST, 

1  LTYPE(20),    WGTX(20),SPDX( 20),    LLANE(20),    TIMEX(20), 

-2 XP0SXl5O-)+-LNUMXt50)  ,    ME4-TX.(-504,    DXEOSX(  50)  ,ACCLRXL5Q) 


DYT=3.0 

YP(1 )=YCL(1)-DYT 
YP(2)=YCL(1)+DYT 

YP(3)=YCL(2)-DYT 
YP<4)=YCI  12J+DYT- 


_Q010_ 
0011 
0012 
0013 


0014 
O0J.5- 
0016 


0017 
-OOiS- 


0019 
0020 
0021 
0022 
002  3 


100 


X.MJL_=_-1.D 

DO    100    J  =  l,4 
DD-J.00    I=1„5Q 
WP(I,J)     =    0.0 


0024 
0025  - 
0026 
0027 
0028 

OD29 

00  30 
0031  —  _ 


1  =  1 
DO 
IF 


1    II=L,NA 
(XPOSX( II  ) .LT.0.O.OR.XPOSX( II  ).GT.XLB) 


GO    TO    1 


IF    (XMl.EO.XPOSXl  I  I  ) )    GO    TO 
_i£_4-LNUMX<  I  14-^0^2) -GO  J0--3 
WP(I ,1 )=WEITX( II  1/2. 
WP(I,2)=WEITX(II )/2. 
WP(I ,3)=0.0 
WP(  I,4)  =  0.0 
GO    TO   4 


WPU  ,1)=0.0 

WP(  I,2)  =  0.0 

WPU  ,3)=WEITX(  II  )/2. 

WP(  I  ,4)=WEITX( I  I )/2. 

CONTINUE 

XP  (  I  )--=-  XJ>QSX4  1 14  


1  =  1  +  1 
GO    TO 


0032 
0033 
D034- 


WP(  I-1,3)=WEITX( II )/2. 
WPU-1  ,4)=WEITX(  II  )/2. 
_XMi  =  XPOSXUI  ) 
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ORTRAN  IV  G  LEVEL   21 


PNTLOD 


DATE 


73170 


00/31/27 


0035 
Q036 


1  CONTINUE 
NAP=I-1 


00  37~ 
0038 
0039 


0040 

00-< 

0042 

^043- 

0044 

004  5 


0046 
004J- 


0048 
0049- 


0050 


0051 

005? 

0053 

00  54, 

0055 

0056 


0057 
Q0_5_£ 


J_F  (NDBUG(5) .EQ.O)   GO 
WRITE  (10,10)  NA,NAP 


TO  30 


10  FORMAT— ('-Q--.WL  DO WA=_L_,_14^" 


NAP=«  t-lA) 


WRITE  (10,1010) 
LOIO  FORMAT  ( »0YP   • ) 

WRITE  (10,11)  YP 
L4._FQ.RMAT  <  1X-,4E-L6.4  ) 

WRITE  (10,1011) 
1011  FORMAT-  (  '  0--XP-  -1-X. 

WRITE  (10,12)  XP 
—  12-  FORMAT  tLX  ,  1 OE  11.3.) 

WRITE  (10,1012) 
-10X2.-F0RMAT  (  »Q  XPOSX'  ) 

WRITE  (10,12)  XPOSX 


WRITE  (10,1013) 
D_L3._FJ1I_.MAJ  .('0  WP  •) 

DO  6  1=1 ,NAP 
__.6._WRI_I_E_  110,13).  (WP(I,J),J  = 
13  FORMAT  (1X,4E16.4,I6) 
30  CONXIMIE _ 


1,4), I 


RETURN 
.END ....  _ 
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GRDLOD 

The  GRDLOD  routine  calculates  the  force  vector  F  due  to  the 
point  loads  wf  on  the  bridge  deck.   The  point  loads  are  dis- 
tributed to  the  nearest  grid  points  on  the  basis  of  area. 
The  components  of  the  force  vector  due  to  one  load  are  cal- 
culated as 


2,2 


f  r* 


2,1 


=   ""i  (X  " 

a- 


Ax,  '  Ay. 


1,1 


1/2   = 


2,1   =   -w 


2,2 


AX 

*  AY 

Ax, 

•  Ay2 

AX 

•  AY 

Ax2 

•  Ay1 

AX 

•  AY 

Ax2 

•  Ay2 

AX 


AY 


) 

) 

) 
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All  forces  acting  on  the  same  point  are  summed  for  that  point. 

If  the  span  is  a  simple  span  the  bending  moments  are  calculated  as 


>,j  =  E(> 


i=l 


at  X   for  a  load  at  X. . 
m  1 


and  the  static  bending  stresses  SIGB.  .  are  calculated  by 
dividing  MB.    by  the  lower  cross  sectional  modulus  ZBL 

(below  the  neutral  axis) . 
No  subroutines  are  called  by  this  routine. 
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GRDLOD 


ZERO 

F(I,J) 

BLOCK 


ii  =  1 


jj  =  1 


£ 


NJ 
NI 


33 
ii 


Yes 


ILO  =  1 
IEND  =  IS(1) 
L  =  1 


H 


Figure   28.    GRDLOD  Program  Flow  Chart 
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V 


DX  =  X(L)/2 
DX1=X' 


DX2=X 


NI 
NI 


dx3_xilo"x'ni 
dx4=xilo-x'ni 

IX  =  ILO 


DY  =  1 
DY1=  1 
DY2=  1 
DY3=  1 
DY4=  1 
JY  =  1 


Yes 


Yes 


Figure  28. 


GRDLOD  Program  Flow  Chart  (Continued) 
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DY=1 

DY1=1 

DY2=1 

DY3=1 

DY4=1 

JY=     NB+1 

C3=    1 


■zz 


Yes 


DY=YB(J)-YB(J-1) 
DY1-Y]^J-YB(J-1) 

DY2=YB(J)-Y' 
NJ 

DY3=YB(J)-Y' _ 
NJ 

DY4=Y'    -YB(J-l) 
N  J 

C3=l/3 


J 
JY 


=    2 

=   J 


(r 


J    =    J+l 


Yes 


±. 


PRINT    OUT 

ERROR 

MESSAGE 


si 


STOP 


"  .  Y  VALUE 
(TRACK)  NOT! 
ON  BRIDGE 
DECK" 


Figure   28.    GRDLOD  Program  Flow  Chart  (Continued) 
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A=DX . DY 

Al= 

=DX1. 

DY1 

A2= 

=DX2. 

DY2 

A3= 

=DX3. 

DY3 

A4^ 

=DX4. 

DY4 

IX 


=C3  (1-—  1W' 
,JY-1  ^J  \£    A  J   NI, 


+F 
NJ   IX,JY-1 


JL 


Figure   28.    GRDLOD  Program  Flow  Chart  (Continued) 
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Yes 


No 


A2 

fix/jy=c3 (1"a-)w'ni,nj+fix-i/  JY-1 


No 


k. 


Al 
F  =C3fl-— }W  +F 

riX=l7JY-l    ^     KX    A    '        NI,NJ       IX-1/JY-1 


i. 


Figure  20.    GRDLOD  Program  Flow  Chart  (Continued) 
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No 


JL 


DX  =  AX(L)/2 

DX1=  x'ni"xiend 
DX2=  x'ni"xiend 

DX3=DX-DX1 
DX4=DX-DX1 
IX=IEND+1 


DX=AX(L) 
DX1=X'NI-XI_1 

dx2=x'ni-xi-i 

|DX3=Xi-X'ni 

DX4=XT-X' 


h 


NI 


IX=I 


Figure  28.    GRDLOD  Program  Flow  Chart  (Continued) 
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M 


IL0=IEND+1 
IEND=IEND+IS 

1  =  1+1 


1  +  1 


Figure   28.    GRDLOD  Program  Flow  Chart  (Continued) 
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CHANGE  SIGN 
OF  F 


Yes 


M 


N 


PRINT  OUT 

PROGRAM 

PARAMETERS 


J  =  1 


-> 


i  =  1 


M  =  1 


±- 


O 


1 


X,YB,F 


NO 


H  p 


Figure   28.    GRDLOD  Program  Flow  Chart  (Continued) 
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No 


YPfl 


UB     .=(XSUP9-Xm) 
mj  2.      m 


Uc/i    1)+MB  MB     ,=X    ( 

\XSUP~  /         m,j  m}      m\ 


4L 


fF.  .  (XSUP„-X. ) 
1/3 2 


XSUP. 


•9 


+MB 


m: 


Figure   ^8^    GRDLOD  Program  Flow  Chart  (Continued) 
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CALCULATE  BENDING  STRESSES 
MB 


SIGB.  •  = 


12. 


ij   ZBL 


n,D 


Yes 


PRINT  PROGRAM 
PARAMETERS 


I   RETURN 


*      MB,  SIGB 


Figure  28.     GRDLOD  Program  Flow  Chart  (Continued) 
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ORTRAN  IV  G  LEVEL 24 


ORDLOO 


DATE  =73171 


21/01/13 


0001 


0002- 


SUBROUTINE  GRDLQD 

THIS  SUBROUTINE  CONVERTS  POINT  LOADS  TO  GRID  LOADS 

REAL    MB,    MOT, MB _..___      .,_    _ 


0003^ 


IEVENT, 
STEADY, 


IN,     10,     ISS*     ITER, 

T,    WEIGHT,     XLB,    YBCL , 


00Q4 


-  COMMON    A&IK2/    DT,  THRSH,       IE, 

1         NAP,    NB,    ND,    NI,    NL,    NSI,    NSUP, 
-2-      NXPTS,   NDBUGtlOW 

3         DELH(IO),    DELXUO),    EEI(10,10),     EIL(IO),     IS(10),    NCS(IO), 
_4 XSUP(IO),    YB(10>,    Y4M4),    YIU.41QU  2CL440»10)  +   ZCUUO.lQl, 

6         E(3,10),    EI(3,10),    XCS(3,10),    ZBL(3,10),    ZBU(3,10), 
-J- Ft80*104, -MO  111,20,10),    MDT(11,10),    WP(50,4),    X(80),XP(50) 

C0MM0N-7BUC4/_  XK< 8Q*8Q) ,XKI ( 80,80) ,    ATEMP(  80,  80)  ,    DTEMP(80) 

1  ,MI(20),    MB(20),    SIGB(20,10) 


0005 
0006 


OOQZ 


0008 
0009 


DIMENSION    SK(20,20,10),     SKP(10,10) 

EQUIVALFNCF  _t_SK ( 1 1 1 ,  1  > » ALEMP  (4, 1J  ) ,    i  SKP(  1,1 ) ,  ATEMPt  1,51)) 

ILO=L  _ 

IEND=IS(1) 

DO    1    L=-l,NSI 


0010 

0011  - 

0012 

0043 

0014 

0015 


DO   2    I=ILO,IEND 

J10-^^.=4-,NB    

F(I,J)=0.0 
4L0=1END+1 
IEND=IEND+IS(L+1) 
DO    3    11=1, NAP 


0016 
0017 
0018 
0049- 
0020 
0021 


DO    4    JJ=1,4 
-NJ=JJ 
NA  =  II 

IF  4WP(II,JJ).EQ.0.0) 
IL0=1 
_IEI40=NXPTS 


GO  TO  4 


0022 
0023 
0024 
0025 
0026 
0027 


0028 
0029 
0030 
0031 
0032 
0033 
0034 

0035 

0036 
0037 
3038 


IF  (XP(NA).GT.X(ILO))  GO  TO  10 
—  IL1  =  0 

DO  5  L=1,NSI 

■  III    =    IL1  +  IS(L) 

IF  (IL0.LT.IL1)   GO  TO  55 

5  CONTINUE - 

55  CONTINUE 

DX  =  DELX(L)/2. 

DX1=XP(NA) 

DX2=DX1 

DX3=X(IL0)  -  XP(NA) 
0X4=0*3 

IX=ILO 


11  IF  (YP(NJ).GT.YB(D) 

aval* 

DY1=1. 
0Y4  =  1^ 


0039 


DY2=1.0 


GO  TO  6 
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ORTRAM- IV  .G  LEVEL  -21 


GRDLOO 


DATE  =  73171 


21/01/13 


0040 

0041 

0042 

0043 

0044- 

0045 

0046 

0047 

0048 
0049- 


0050 
0051 
0052 
0053 
0054 
0055- 


0056 
0057 
0058 

0059 
0060 


0061 
0062 
0063 
0064 
0065 
QQ66__ 


0067 

006  8 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0073- - 

0079 

0080 

0081 

0082 

0083- 

0084 

0085 

0086 

0087 

0088 

0089 

0090 


0Y3=1.0 
C3  =  1.0 
JY  =  1 

7  A=DX*DY 

— A1=DX1*DY1 
A2=DX2*DY2 
A3=DX3*DY3 
A4=0X4*DY4 


IF  (IX.EQ.NXPTS+1)   GO  TO  80 

IE  UY.EQ.JMB*!-)   GO- TO  ai  .     .  

F(IX,JY)=(1.-A3/A)*WP(NA,NJ)*C3+F(IX,JY) 

IF  (JY.EQ.l)   GO  TO  80 

F(IX,JY-1)=(1.-A4/A)*WP(NA,NJ)*C3+F(IX,JY-1) 

IF  UX.EQ.l)   GO  TO  4 

IF  (JY.EQ.NB+1)   GO  TO  82 

EXIX-T-lWY-J=ll»-A2yA)^WI>{NA,MJ)*C3+F(IX-l,JY)  .. 

IF  (JY.EQ.l)   GO  TO  4 

F(IX-1,JY-1)=(1.-A1/A)*WP(NA,NJ)*C3+F(IX-1,JY-1) 

GO  TO  4 


ai 


80 


82 


8 


100 


>  IF  (YP(NJ).LE.YB(NB))  GO  TO  8 
_-QY  =  l^ 

DY1=1.0 

DY4=1.0 

0Y2=1.0 

DY3=1.0 

C3  =  1.0 
_JX=MB_+1 . 

GO  TO  7 

J0=2 

00  9  J=JOtNB 

JY=J 

IF  (YP<NJ).GT.YB( J))  GO  TO  9 

DY=YB(J)-YB(J-1) 

DY1=YP(NJ)-YB< J-l) 

DY2=YB(J)-YP(NJ) 

DY3=YB(J)-YP(NJ) 

DY4  =  YP(NJ)-YB(J-H 

C3  =  1.0/3.0 

GO  TO  7 

CONTINUE 

WRITE  (10,100) 

FORMAT  (1X,«  Y  VALUE  (TRACK)  NOT  ON  BRIDGE  DECK  ') 

STOP 


la  IE    1XP  LNAJ..LT  .X  (  IEND  )  > 
DX=DELX(NSI )/2. 
DX1  =  XP(NA)    -    X(  IEND) 
DX2=DX1 
DX3=DX    -   DX1 
DX4=DX3 
IX=IEND+1 
GO    TO    11 


GO    TO    12 
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QRTRAN    IV    G   LEVEL 21- 


__GRDLOD 


DATE    =    73171 


21/01/12 


0091- 

0092 

0093 

0094 

00*5- 


12 


0096 
0097^ 


0098 
0099- 


75 
Z6 


0100 
0101 


IL0=IL0+1 

DO    13    I=ILO,IEND 

IF    IXP(NA).GT.XU))    GO    TO    13 

IL1    =    0 

DO    75-L=l,NSl 

IL1    =    IL1    +    IS(L> 

IF    (I       .LT.IL1)      GO    TO    76 

CONTINUE 

CONTINUE  

DX=DELX(L) 

D_X1=XP(NA)-XU--U 


0102 
0103 
0104 
0105 
0106 
0107 


0108 
0109 
0110 
0111 
0112 
0113 


0114 

0115_„ 

0116 

OUT 

0118 

0119 


0120 


DX2=XP(NA)-X(I-1) 
DX3=X( I ) -XP (NA) 

DX4=X(I)-XP(NA) 
U=J_  _ 

GO   TO    11 
13    CONTINUE 


4  CONTINUE 

_  __3_.  CONT 1NUE      ._ 

IF     (NDBUG(6).EQ.0>    GO    T019 

WRITE    (10,102) 

102    FORMAT    ( IX , »GRDLOD • > 

JMA=0— 

DO    14    L=ltNSI 

NA=NA+1S(LJ 

DO  15  1=1, NA 
DO  15  J=1,NB 
WRITE     (10,103)    X(I  ),YB(J),F(I,J) 

FORMAT    (1X*JlE16^4) 

19   CONTINUE 


14 


15 
103 


__C_ 
C 


SIMPLE    SPAN    BENDING    MOMENTS    AND    BENDING    STRESSES 


0121 
0122- 


IF     (NSUP.GT.2) 
XS1--_1SUJ 


GO    TO    50 


0123 
0124 
0125 
0126 
0127 


0128 
0129 
0130 
0131 
0132 
0133- 


DO    45    J=1,NB 
DO    20    1=1, IS1 
SIGB(UJ)    =    0.0 
MB  (I  )         =0. 

^2-©. --&QNT  IWUF- -  -     - 

DO    30    1=1, IS1 

TANTHA    =     (F(I,J)*(XSUP(2)    -    X ( I ) ) ) /XSUP( 2) 

TANTHB    =    (F(I,J)*X(I))/XSUP(2) 

DO    25    M=1,IS1 

IF     (X(M).LE.X(I ))       GO    TO    23 


0134 


_MB4M4 
GO    TO 


=  -(XSUP-(2)-0UMU*-TANTHB   +    MB(M) 


25 


0135 
0136 
0137 


0138 


23 
25 
30 


MB(M) 

CONTINUE 

CONTINUE 


=    X(M)*TANTHA    +    MB(M) 


DO   45    1=1, IS1 
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QRTRAN  IV  G  LEVEL 21- 


GROLOD 


DATE  =  73171 


21/01/12 


0139 
0140 
0144 
0142 
0143 


41 


NU    =    NCS(J) 
DO    41       N=1,NU 
IF(X4I).LE.XCS(N,J)) 
CONTINUE 
_£I0P_44- 


GO    TO    42 


0144 


-42   CONTINUE 


0145 
0146 
0147 


0148 
01  49 


0150 
0151 


0152 
0153 


0154 
01  £5  _ 
0156 
0157 


0158 
0159 


SIGBUtJ)    =    (MB(I)       /ZBL(N,J))/144.0 

IF     (NDBUG(6).EQ.O)       GO    TO    45 

WRITE    (  Tn,?0A5J-L,..),N,mUJ ...  ,ZBLIN,  d)  ,  SIGBil,  J  ) 


,316, 


MB 


2045    FORMAT    ( •    I    J   N 

WRITE    (10,1040) 

WRITE    (10,1041) 

45_  £0N T  I  NU  E 

IF    (NDBUG(6).EQ.O)      GO    TO    50 

WR1T-E    (  10  ,10421 

WRITE    (10,1043)  SIGB 

_10AQ_  FORMAT    l'OHB         •  ) 

1041    FORMAT    (1X,10E13.5) 
__1042  FHRMAT_J'0  SIGB    •) 

1043    FORMAT    (1X,10E13.5) 
SO   r.flNTINUF 


MB      ZBL         SIGBS3E12.4) 


0160. 
0161 


RETURN 
END 


156 


LODDIS 


This  subroutine  develops  the  actual  discontinuity  moments, 
i.e.,    due  to  the  imposed  live  load.   The  constant,  during  a 
given  run,  unit  load  discontinuity  moment  vectors  are  operated 
on  by  the  applied  force  vector,  i.e., 


M 


DT 


J  3 


IS(L) 


E 


i  =  1 


hi 


M. 


ID 


No  subroutines  are  called  by  this  routine. 
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LL  =2 

I  =IS(2) 
o 

i  =  1 


-* 


MDT  T  T  . =MDT   . +F .  . *MD .   . 
LLD     LLD   ID    !/D 


f 


i=l+2 


Figure   29.    LODDIS  Program  Flow  Chart 
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LL  =LL+1 
I0=Io+IsLL-l 


Yes 


PRINT  PROGRAM 
PARAMETERS 


RETURN 


MDT 


Figure   29.   LODDIS  Program  Flow  Chart 
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-ORTRAN    IV    G    LEVEL       21 


OOOl 
00  0-2 


SUBROUTINE    LODDIS 
REAL    MD,    MDT 


LODDIS 


DATE 


73170 


00/31/2 


00O3- 


-  COMMJDN    /BLK2/    DT  , 
1         NAP,    NB,    ND,    NI,    NL, 
Jg NJU>JS-r-ND8UG4-104  r 


3 

6 
.2- 


THRSH,   IE,  IEVENT,  IN,  10,  ISS,  ITER, 
NSI,  NSUP,  STEADY,  T,  WEIGHT,  XLB,  YBCL, 


DELW(IO),  DELX(IO),  EEI(10,10),  EIL(IO),  IS(10),  NCS(IO), 
XSUPQQ),  YB(10),  YP(4),  YCL(IO),  ZCL(10,10),  ZCU(10,10), 
E(3,10),  EI(3,10),  XCS(3,10),  ZBL(3,10),  ZBU(3,10), 
EXBOflQ),  MD(11#J20,10),  MDTX11,  10),  WP(50,4),  X(80),XP(50) 


0004 


DO  50 0=1^*1*6-- 


0005 

DO  10   LL=1,NSUP 

0006 

C 

10 

MDT(LL,J)  =  0.0 

0007 

0008 

10  =  1 

11  =  0 

oonp 

C 

nn  30  i.i  =?,nst 

0010 
0011 



10  =  10  +  11 

11  =  11  +  IS.lLL-1) 

0012 
noi  3 

DO  20  1=10,11 

MDTILI  ,J)  =  MDTd  L,  J)  + 

0014 

C 

20 

CONTINUE 

0015 
0016 

30 
50 

CONTINUE 
CONTINUE 

0017 

C 

IF  (NnniJG.(7).E0.0)   GO 

0018 

QQ19. 

0020 

J)£L2:L_   __ 
0022 

1001 

L0Q2 
60 

c 

WRITE  (10,1001) 
FORMAT  («OMDT») 
WRITE  (10,1002)  MDT 
FORMAT  (1X,11E11,3) 
CONTINUE 

0023 
0024 

RETURN 
END 

JFJ_U.J.L_*.  MD(JLL^I,  J) 


.id.  6  a 
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ASEMBL 

The  subroutine  ASEMBL  calculates 
The  internal  moments 


MI.   .  T 


The  bending  moments 


X   —  X 
IB.   .  _  =  MI.   .  +  11- |MDTL 

xsup]^1-xsup. 


:) 


x.  -  xn 

1  MDT. 


XSUPL+l"XSUPI      '  L+1/J 


Deflections 


6 


X.  MB.   .  AXT 

i    i>]    L 


i/j    /    '   \         EEI 


Lj 


) 


The  subroutine  STIFF  is  called  to  calculate  the  longitudinal 
and  lateral  deflections,  include  the  discontinuity  con- 
tribution, to  assemble  the  flexibility  matrix  per  span 

(free  body) . 
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The  following  subroutines  are  called 
STIFF 
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(   ASEMBL 


ZERO  BTEMP 
X_=XSUP 


o      L 
X  =XSUPn 


DEL. .=0 


XL 


K  =  II 


-> 


v 


TANAF  =  F 


TANBF= 


(X1-Xi+I1_1) (X(k)) 
kf j      (Xi"V 

(x,l7x(k))x.+i1.1 

k,j     (X--X  ) 

J.   ..  Q 


jk. 


Figure  30.     ASEMBL  Program  Flow  Chart 
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MI.      .=MI.      .+    TANAF 
i/D  I'D     . 


MI.      .=MI.      .+TANBF 
l/D  i/D 


MB 


.  MDT 
-XSUP_  /  L,j 


•  =MI„      .  +  U-Xi+I1"X° \ 

1,3         1/3     \      XSUPL+1-XSUPL  / 

/Xi+Il-XQ  w  \ 

+      \ XSUPL+1-XSUP  J(  MDTL+1 , j/ 


MB 


SIGB 


l,]  ZCL1 


JLiJL 


Figure    30.      ASEMBL   Program  Flow   Chart    (Continued) 
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Figure  30.    ASEMBL  Flow  Chart  (Continued) 
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:£RTI 
0001 
O002 
0003 


_2JL. 


ASEMBL 


DATE  =  73171 


21/01/1 


SUBROUTINE  ASEMBL (L) 
REAL  MI, MB, MOT 


COMMON  /BLK2/  DT, 


NL 


.ItiR£fcU_ -IE.,. 
,    NSI,    NSUP, 


JEVEbHU_  IN,.  10,    ISS,    ITER,-. 
STEADY,    T,    WEIGHT,    XLB,    YBCL, 


EIL(IO),     IS(IO),    NCSUO), 


NAP,    NB,    ND,    NI, 

NXPTS,    NDBUG(IO), 

DELW(IO),    DELX(IO),    EEI(10,10), 

iUPtlOX,    YB(10),    YP(4)^  YCLUO),    ZCL(10,10),    ZCU(10,10), 
E(3,10),    EI(3,10),    XCS(3,10),    ZBL(3,10),    ZBU(3,10), 
Pt»nrmit   Mnint?ntini.    mdtmi.io).   wp(50f4)t  *(  80)  ,XP 150). 


0004 


0005 
000ft 


COMMON./BLK4/   XM80, 80) ,XKI (80,80) ,    ATEMP( 80, 80) ,    DTEMP<80) 
1  ,MI(20),    MB(20),    SIGB(20,10) 


DIMENSION    SK(20,20,10),    SKP(10,10) 
FQI1IVAI..FNCF_1S]<LL1,-Ulij-ATEMP11?1)  ) ,. 


(.SKPJJL,1IIJ.ATEMpa.f5JJJL 


C 
C_ 


0007 

_0008_ 

0009 


II    =    0 

XQ__=    XSUP(L) 

XI    =    XSUP(L+1) 


0010 

ooTi" 

Q012 
0013 
001* 


10 


0015 

C 

0016 

C       _ 

C 

.00X2 

0018 
0019 
0020 

0021 .__     L6 

0022 

0023 . 

0024 

C 
0025 

C  ..    ... 

0026 

0027 


LM1    =    L   -    1 

IF~ (LM1.EQ.0)       GO   TO    11 
DO    10    1  =  1, LM1 
11=11    +    ISC  I  ) 

10  =    II    +_LS(L) 

11  =    11+    1 

ISL    =    IS(L) 


..DO  1Q0   J    =   A,NB.         __ _..     

DO    16    1=1,20 

MB(I)    =    0. 

Mid  )  =    0. 

CONTINUE 

DO    50    1=11,10 
_IANIHA_  =    LF  II  ,JU*lXlrX(I_)  )  )  /(XlrXO) 

TANTHB    =    (F(I,J)*    X( I ) ) /( Xl-XO) 

KK    =    0 

DO    20   K=I1,I0 

KK    =   KK— +_ 1 


0028 
0029. 
0030 

0031 
00^2 
0033 


IF  (X(K).LE.X(I ) > 
MIIKK)  =  MKKK)  + 
GO   TO   20 


GO  TO  15 
(X1-X(K))*TANTHB 


.  £_ 


15 


CONTINUE 
UliKK-L  =. 


20  CONTINUE 


MLKK1  +.X(KJ*TANTHA 
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G  LFVFL 21 


ASEMBJ__ 


DATE  =  73171 


21/01/1: 


0034 


_5.CL_C.0NT_  I NUE 


0035 
0036 
0037 


0038 

0039 
0040- 


II  =  0 

00  70  1=11,10 
_n  =  II  *^4 


MB  (II)  =  MI  .II )  + 
-I + 

SIGBIIrJ-  =  (MB(II) 
-__.__J1X0NI  INU.E 


(1.0  -(X(I)      -X0)/( X1-X0) )*MDT(L,J) 
((X(IJ     .  -X0)/(X1-X0) )  *  (  MDT(L+1,J)) 
/ZCL(L,  J)>/144.0 


-QOAJL 


0042 

_0043_- 

0044 

-004-5— 

0046 

-0047 — 


0048 
-004-9- 


n__XNnfl»f_i3-i---E.Q . 0 ) Gn  Tn  1  oo_ 

WRITE  (10.1040) 

WRITE  (10.1041)   MB  

1040  FORMAT  CO  MB    •  ) 
104.1_-F_0J_.MAT.'  UX.LQE.1_U5.L-. 
WRITE  (10,1044) 

WRITE— U-ULQ41) MI—   — — 

1044  FORMAT  CO  MI    •  ) 
UNJJE-   ._ 


-0050- 

0051 

-0052- 


0053 

-00-54-- 

0055 


0056 


l_L.J.NDBUG(3).EQ.O)   GO  TO  60 

WRITE  (10,1042) 

WRITF  (1Q<1043) SJ-GB 


1042  FORMAT  ( »0  SIGB  *  ) 
1043_  £QRMAT_QX*10E13.5) 
60  CONTINUE 


CALL  STIFF  (L) 


0057 
_QQ5_B__ 


RETURN 
ELND 
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BRGDYN 

This  subroutine  calculates  the  static  and  dynamic  deflections, 
the  dynamic  load  factor,  and  the  dynamic  bending  stress. 

Given  the  nonhomogeneous  equation  of  motion,  for  forced 
undamped  vibration, 

[Fk]  =  [K]  [6t]  +  [M6t] 

where  t  is  the  time  increment  index,  that  is,  time  t  occurs 
at  t  and  time  t-At  occurs  at  t-1,  K  is  the  constant  assembled 
stiffness  matrix  for  a  span,  M  is  the  mass  of  each  element,  6  the 
vertical  position  of  each  element,  6  the  vertical  velocity  of  each 
element  and  5  the  vertical  acceleration  of  each  element. 
The  positive  vertical  direction  is  downward. 

In  order  to  minimize  computer  time  a  presolved  simultaneous 
solution  of  the  numerical  integration  equations  and  the 
equation  of  motion  has  been  utilized.   A  double  trapezoidal 
integration  formula  was  used  as  a  basis,  that  is, 

i«ti  -  [«t-i]  +  t-  lK]  +  ^t-i1 

and 

At 


[6t]  =  [6t_±]  +  £§-   [jt]  *  ta^] 


We  first  multiply  the  flexibility  matrix,  [K]   ,  through  the 
equation  of  motion, 

[K]_1[Ft]  =  [6t]  +  [K]_:L[M6t] 


Solving  for 


[6t]  : 
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[«t]  =  [k]  1    j  [Ft]     [M5t- 


Substituting  in  the  integration  formula 


'Vi1  +  4  j1^1  +  [l-i]|  = [K]_1  ![Ft]  -  [aK] 


Solving    for 


[St]  = 


[K]    -    At    (tK]-1jtFfc)     -    [M6t]   j   -     [6^])-    [6^1 
Substituting    for    [6. ]    in   the    remaining    integration    formula 

lK-l]    +  %   j[S'tI    +    "t-l'j-   At    ([W-1|cFt]    -    [M6t]j 


[6t-l]  -  [6t-li 


Solving    for 


[6t]: 

[6t]    = 

[I] 

4M 
At2 

4 

• 

>  I  AtZ  C  At 

At    ^t-l'    "    ^'t-11! 

Subsequently,  [<S  ]  and  [S._]  are  determined  from  the  integration 
formulas . 


Hence,  an  explicit  solution  of  [6,  ]  exists  which  is  equivalent 
to  an  iterative  double  trapezoidal  numerical  integration  of  the 
nonhomogeneous  differential  equation  of  motion. 
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The  dynamic  bending  stresses  are  then  calculated  as 


a  (t) .  .  =  DLF.  .  asB. 

where  asB.  .  are  the  previously  calculated  static  bending 
stresses. 


This  subroutine  does  not  call  any  subroutines. 
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[6(t~At)] 

[«(t>] 

[5 (t-At)] 

[6(t)] 


Figure   31.    BRGDYN  Program  Flow  Chart 
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Calculate : 


[6(t)]  =  — i-j  ][6(t-At)]  +  At  [6(t-At)]  + 


X  4M 


At' 


(» 


S(t-At)]  + 


|  [K]-1[F(t)])| 


[6(t>]  =  j|  j[6(t)]  -  [6(t-At)]|-  [6(t-At)] 


[6(t)]  =  ~ •  {  [6(t)]       :   t)  \)        [6(t-  t)] 


j& 


Figure  31.   BRGDYN  Program  Flow  Chart  (Continued) 
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Yes 


DLF. .=0 
13L 

a(t)ijL=0 


J  =  1 


\/_ 


i  =  1 


6(t)  .  . 
D   .  .  =  ^ 

lfijl   6  . . 

sijl 
0(t)iJL=  DLFijLaBijL 


vJ*' 


1  -i •+  1 


No 


J  =  j+1 


No 


Figure  31.     BRGDYN  Program  Flow  Chart  (Continued) 
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^QRT-RAN--IV    G  -LEVEL-     21- 


8RGDYN 


DATE    =    73171 


21/01/1 


0GO1 


4Ki©2- 


C 
& 
C 
-G 
C 


SUBROUTINE    BRGOYN(L) 

THIS    SUBROUTINE    CALCULATES 

U    THE    STATIC    DEFLECTIONS 

1    THE  ^>WAMTG   DEFLECTIONS 

3    THE    VELOCITY    OF    THE    DECK 

^  THR    ACCELERATION    OF    THE    DECK 

5    THE    DYNAMIC    LOAD    FACTOR    AT    EACH 

6-T44E    DYNAMIC    BENOING    STRESS 

REAL    MO,    I40-T- 


POSITION    AND    AT    EACH    TIME 


0003 


LOGICAL    STEADY 


0004 


COMMON    /BLK2/    DT , 


THRSH,       IE,     IEVENT,     IN,     10,     ISS,     ITER, 


0005 


0006 


0007- 


000  8 
0009 


0010 
0011 


04U2- 

0013 

0014 


4.         NAP,    NB,    ND,    NI,    NL,    NSI,    NSUP,    STEADY,    T,    WEIGHT,    XLB ,    YBCL, 

2         NXPTS,    NDBUG(IO), 
-3 -&£LW-U-0-)-^-DELX(  1W,    EEtU0,10i,    EIL(IO),     IS(10),    NCS(IO), 

4         XSUP(IO),    YB(10),    YP<4),    YCL(IO),    ZCL(10,10),    ZCU(10,10), 

6  £13, LO),    EI(3,10),    XCS«3,10),    ZBL(3,10),    ZBU(3,10), 

7  F(80,10),  MD< 11,20,  10),  MDT(11,10),  WP(50,4),  X(80),XP(50) 
COMMON  /BLK4/  XK( 80,80) ,XKI ( 80,80) ,  ATEMP( 80, 80) ,  DTEMP(80) 

1  ,MI(20),  MB(20),  SIGB(20,10) 


COMMON  /BLK5/  S  IG{80) ,DFLN ( 80) , DFL( 80, 10) , ITIME,KTIME , 
.1      SIGJM120,3,10>,  SIGC( 10,3,10),  SIGTHS,  DSIG 

2  ,  I0PT1,  I0PT2,  I0PT3,  I0PT4,  NPTS(IO),  PTS(3,10) 

3  ,  IPTS(3,10) 

_  .C0MMHN_^BLK.7/W&T,V£L^DELT,BLEN,NEV,NNL,NBEAM,NP(  10)  ,DMM(3)  ,ITAPE2 
1,SIG0UT(3,10) 

DIMENSION  DLF(80),     DFLD(80,10),  DFLDD( 80, 10 ) , 
1        DFLND(80),  DFLNDD(80),  DFS(80),  FF(80),  DUM(80) 
DATA  ISTDY,JSTDY/0,0/ 

IF(L.EQ.l)  ISTDY  =  0 
ITER  =  3 

INITIALIZE  FORCE  VECTOR  FOR  L  TH  SPAN 

_  IX  =-  0~ 
LL  =  L  -  1 
IF(LL.EO.O)  GO  TO  11 


0015 
0016 


DO  1  1  =  1, LL 
1  IX  =  IX  +  IS(L) 


0017 
0018 
0019 
0020 


11  CONTINUE 

IXU  =  IX  +  IS(L) 
IX  =  IX  +  1 
II  =  1 


O021 
0022 


DO  2  J=-l  ,NB 
DO  2  I=IX,IXU 
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QRIAAiL  UL   G  LFVEJ 21. 


BRGOYM. 


DATE  =  73171 


21/01/1: 


00 13 
0024 

0025 


FF(II)  =  F ( I , J  ) 
2  II  =  II  +  1 

KTER  =  0 


0026 

0027 
0028  - 
0C29 
0030 

0032 


c 


IF  (NDBUG(9).NE.l)   GO  TO  3 

WRITE  (10,1000)  L 
OOQO  -EORMAI I 'OFORCE  VECTOR.  FOR  SPAN  L  =',I4) 

WRITE  (10,1001)  FF 
mm  FDRMAT  MXf  10MO.O) 


J3__NDB  =—lS  LU  *N  B _  . 

ISL  =  IS(L) 

£ 

C    IF  STEADY  STATE  CLEAR  DEFLECTION  STORAGE 


0033 
_QQ34_. 
0035 
_QQ36 ._ 
0037 
J1Q2B 


IF  (ISS.NE.O)   GO  TO  4 
__DQ_5  L=l.t80 

DFS(I)  =  0. 

DFLN(I)  =  0. 

DFLND(I)  =  0, 
ELELNDD(I)  =  0  . 


0039 
Q040 
0041 
0042 
0043 


DO  5  LL=1,NSI 
DFL(I,LL)  =  0. 
DFLD(I,LL)  =  0. 
5_DFLDD(IrLL)  =  0. 
STEADY  =  .FALSE. 


0044 


0045 
0Q46 


4  CALL  MATMULtXKI,  FF,  DFS,  NDB,   1,   NDB,  80,  80) 
_C 

c 

C         SOLVE    NON-HOMOGENEOUS    DIFFERENTIAL    EQUATIONS    SIMULTANEOUSLY    FOR 
C  EACH    GRID    OF    THE    DECK 

_C_ ._.._.     .    _     _   _ 

CI    =    4./(DT*DT) 
C2    =    4./DT 


D047 
0048 
0049- 

0050 
0051 
0052- 
0053 


0054 

0055. 

0056 

0057 

0058 

00  59- 

0060 


_00   6    I=1,NDB 

6  DUM(I)    =    C1*(DFS(I)    -    DFL(I,D)    -    C2*DFLD(I,L)    -    DFLDD(I,L) 
-       CALL  -MATMUL-  1XK_»    DUM,    DFLNDD,  _NDB,       1»-    NDB,_  80*    80) 

DO    7    1=1, NDB 

DFLND(I)    =    DFLD(I,L)    +    DT/2.    *    (DFLNDD(I)    ♦    OFLDD(I,D) 
__DFLN   11)    =   DFL1I,LJ    +   DT/2   *    (DFLND(I)    +    DFLD(I,D) 

7  CONTINUE 


IF     (NDBUG(9).EQ.O)       GO    TO    8 

WRITE    110,2001) 

WRITE    (10,2002)     ( DFLNDD( I ) , I = 1,NDB ) 

WRITE    (10,2003) 

WRITE    (10,2002)    ( DFLND( I ) , 1= 1,NDB ) 

_WR1T£_1 10*20051     _ 

WRITE    (10,1002) 
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:QRTRAN  IV— G-LEVEL 2X- 


BRGDYN 


DATE  =  73171 


Zl/Ql/i: 


0061 
0062 
0063 
0064 
0065 


0066 

0067 
0068 


0069 
0070 


L002  FORMAT  CODYNAMIC  DEFLECTIONS') 

WRITE  (10,2002)  (DFLN ( I  ) , 1  =  1, NDB  ) 
2001  FORMAT! • ODFLNDD • ) 

2003  FORMAT( 'ODFLND') 

2005    FORMAT  (  M3PFLN-U 

2002    F0RMAT(1QE13.4) 
_C 

1003  FORMAT    (1X,10E12.3) 
WJUTE    (1 0 ,L004 ) 

1004  FORMAT    POSTATIC    DEFLECTIONS') 
H R I T£— U-Q-+-LQ03-)  ( 0 F S (U^, J  =  UNDB)  - 


00-71 
0072 

X1Q73- 
0074 

X*OT§- 


0076 

0077- 

007  8 

OOJSl 

0080 

0081 


-CONTINUE 

DO    10    J=1,NB 
-DB   10    LI=1, ISL 

1=11+    ISL*(    J-l) 

DFLDD(  I  ,L4   =-DF4.^0P444- 

DFLD(I,L)     =    DFLND    (I ) 
-DFLU,L>    =    DFLNUJ 

IF    (DFS(I).NE.O)       GO    TO 

DFL(I-»L1    =    0.0 

SIG(I)    =   0.0 

GO    TO    10 


20 


0OB2- 
0083 
XX084 
0085 
0086 


OX)  87 
0088 
0089 


20    CONTINUE 

DLFU  )    =   DFLNUJ/DFSU  ) 

SIG(I)    =    SI GB( 1 1, J)    *    DLF(I) 

IF(ABS(SIG(I  )  ).GT.THRSH)   ISTDY  = 
10  CONTINUE 


COLLECT  OUTPUT  DATA 


29- 


IF     ( ITAPE2.GT.0)       GO    TO    29 
IFU0PT1.EQ.0.AND.I0PT2.EQ.0.AND. 
-GONTtNUF -  


I0PT3.EQ.0)       GO    TO   40 


0090 
0091 
0092 


DO  30  J=1,NB 
NPT  =  NPTS(J) 
IF(NPT.EQ.O) 


GO  TO  30 


O093 
0094 
0095 


oa 


IP 


_25  4-=4_rNPT— 
=  IPTS(I,J) 


IF(IP.LT.IX.OR.IP.GT.IXU)   GO  TO  25 


0096 
0097 
009S- 


0099 
0100 


25 
30 


IPP    =    IP    -    IX    +    1    +    (J-1)*ISL 
SIGP(KT1ME,I,J)    =    SIG(IPP) 
_-U=_  i4TAPE2,£T  . 04    S 1CQUT ( \+  J )  =S I G  (IPP) 
CONTINUE 
CONTINUE 


0101 


0102 


40    CONTINUE 


IF     (NDBUG(9).EQ.O)       GO    TO    100 
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-ORTRAN    IV    G    LEVEL — 2± 


BRGDYN 


DATE 


73 171 


21/01/1: 


0103 
0104 

0106 

0108 


0109 
0110 


0111 

-04i2- 

0113 

0115 


WRITE    (10,1005) 

WRITE    ( 10,1003) (DLF(I ), 1=1, NDB) 
WRITE    < 10,1006) 

WRITE     U0,1003MSIG(I  ), 1  =  1, NDB) 
*AM4G  404D-  FACTOR  «  ) 
1006    FORMAT    ('ODYNAMIC    BENDING    STRESSES') 


C 

c 


CHECK  FOR  STEADY  STATE 


100  CONTINUE 

IF  tL.NE.NSl) — OO-T-O-104 — 

IF  USTDY.EQ.O.AND.JSTDY.EQ.O) 
4STOY  -=_  ISTDY 

101  CONTINUE 
RETURN 

END 


STEADY  =  .TRUE, 
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BROUT 

The  subroutine  BROUT  prints  out  the  optional  output  on  an  event 
basis,  i.e.,  Stress  vs.  Time  Graphs  and  Stress  Minimums ,  Maxi- 
mums and  ranges  for  specified  points.   Maximum  stress  histograms 
are  also  updated. 
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No 


OUTPUT  GRAPHS  OF 
STRESS  VS  TIME  FOR 
SPECIFIED  POINTS 


Yes 


PRINT  OUT  STRESS 
MIN,  MAX,  RANGE  FOR 
SPECIFIED  POINTS 


MAXEXT 


aMIN  aMAX 
RANGE 


Yes 


STORE  MAX 
-»  STRESSES  FOR 
HISTOGRAMS 


Figure   32.    BROUT  Program  Flow  Chart 
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IJV-G-  LEVEL — 24- 


B  ROUT- 


DATE  =  73171 


21/01/1: 


O004~ 
0002- 
0003 


SUBROUTINE  BROUT 


COMMON  /BLK1/  IPAGE ( 120, 46 ) ,  DUM(880) 


COMMON 


1 

^2 

3 

— 4- 


_/SL^^/__OT^ 

NB,  NO,  NI, 


NL,  NSI,  NSUP, 


I-EVENI,- 
STEADY, 


INf  10,-ISS,  ITER, 
T,  WEIGHT,  XLB,  YBCL, 


NAP, 

NXP-TS*    NDBUG(IO), 

DELWUO),    DELX(IO),    EEI(10,10),    EIL(IO),     IS(10),    NCS(IO), 

XSUPUO), -YB4-10X,    YP(4)»    YCL(IO),    ZCL(10,10),    ZCU(10,10), 

E(3,10),    EI(3,10),    XCS(3,10),    ZBL(3,10),     ZBU(3,10), 

F  (  80  ,4Q-^-MD<-UU,2O+-10  U-MOTtl  1,  10)  ,    WP(SO,4),    Xl£0),XP(50) 


000^ 


COMMON    /BLK5/    SIG (80) ,0FLN (80) , DFL ( 80, 10) , I TIME,KT IME , 


1  SIGP(120,3,10), 

_2-  r    I0PT1,     I0PT2, 

3  ,     IPTS(3,10) 


0005 


0007 


C 
C 


DIMENSION 
DIMENSION 


DATA  NTL 


TIME0(13) 
NTL (2  ,46), 


SIGC( 10,3,10),  SIGTHS,  DSIG 
I0PT3,  I0PT4,  NPTS(IO),  PTS(3,10) 


NDASH(IO),  KPT(3) 


/4H   5,,4H000-,4*4H  ,  4H 

_4«-  4.^4^000-^4*4^  -*4H 

4H   3,,4H000-,4#4H  ,4H 

4H   2,,4H000-,4*4H  ,4H 

4H   1,,4H000-,4*4H  ,4H 

4H     ,4H   0-,4*4H  ,4H 

4H   1,,4H000-,4*4H  ,4H 

.  4H-— 2-r+4MO0G-  ,4*4H  , 4H 


4,,4H500-,4*4H 

3,,4H500-r4*4H 

2,,4H500-,4*4H 

1,,4H500-,4*4H 

,4H500-,4*4H 

-,4H500-,4*4H 

-1,,4H500-,4*4H 

-2,,4H500-/ 


0008 
0009 
-0010 

-0044- 


DATA  NDASH  /1H-, 1H-, 1H-, 1H-, 1H-, 1H-, 1H-, 1H-, 1H-,  1H|/ 
DATA  KPT  /1H*,   1H+,   1H./ 
DATA  IBLANK/1H   / 

IF  (  IOP-T4  ,r£O,0  > GO  TO-WO 


OUTPUT  GRAPHS  OF  STRESS  VS  TIME  FOR  ALL  SPECIFIED  POINTS 
CALCULATE  TIME  VALUES 


-O042- 
0013 
0014 
0015 

0016 

0017 
-001-8 
0019 


00-5  J=4,l  3 

TIMEO(J)  =  DT  *  (  ITIME  -  1 

IPTEND  =  (KTIME+9)/10  +  1 


+  (J-l)*10) 


IF  (NDBUG(8).NE.O) 
1WRITE( 10,2005)  ITIME,  KTIME,  IPTEND 
2005  FORMAT  ( » OT IME  KTIME   IPTEND', 316) 


0020 


DO  50  J=1,NB 
NPT  =  NPTS(J) 
IF  (NPT.EQ.O) 

CLEAR  THE  GRAPH 

DO  10  1=1,46 


GO  TO  50 
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ORTRAN    IV    G    LEVEL — 2-4^ 


BROUT 


DATE 


73171 


21/01/1 


0021 
0022 
0O23 
0024 


DO    7    K=l ,120 
1PAGE(K,I)    =    IBLANK 
7    CONTINUE 
10    CONTINUE 


O025- 

0026 

^>0£7- 


6    DO    9    K=l,120, 10 
DO    9    KK    =    1,10 
—9    1PAGE(KK+K-1,31)    =    NDASH(KK) 


PUT    ON    THE    POINTS 


0028 
JHL29- 


0030 
0031 


0032 


0033 
0034- 


0035 
0036 


0037 

0038 
003.9- 


00    20    1=1, NPT 
00-20  K=1»KTIME 
XL    =    31.0    -    ( (3.0/ 
L    =    XL 


500. )*<SIGP(K,I,J)  -83.0)) 


IF  (NDBUG(8).NE.O) 

-JLWRITE  UO,2Q2Q)  K,  I ,  J,  L,  SIGP  (  K,  I  ,  J) 

2020  FORMATC  K  I  J  L»,4I6,'      SIGP(K,I,J)' 
LF__LL.LTJ.:L)   L=l 

IF(L.GT.46>   L=46 
IPAGF(K,I.  )  =  KPT1IJ ._ 


E14.6) 


20  CONTINUE 


0040 
0041 
0042 
^0043 
0044 
0045 


0046 
0047 
0048 


WRITE  (10,1020)  IE,IEVENT,J,  ( KPT{ I ) , PTS( I , J ) , 1=1 ,NPT  ) 
J.O20  FORMAT  (1HL,  •  STRESS  EVENT'  ,  1 4,  •     LOAD  EVENT', 


14, 


BEAM' ,13, 


POINTS', 3(4X,A1,  F8.3)) 


WRITE  (10,1021) 

WRITE  (10,1021) 

1021  FORMAT(IHO) 
DO  3Q  1=1,46 

WRITE  (10,1030) 
1030  FORMAT  LDU2  A4, 


NTH  1,1  ),NTL(2,I  ),(IPAGE(K,I)  ,K=  1,120) 
i_LL*-l£OAl) .   .. 


30 


CONTINUE 

WRITE  (10,1035)  (TIMEQd  )  ,  I  =  1 ,  IPTEND  ) 
1035  FORMAT  (   • OT IME= • ,F5 . 2, 12( 4X,F6 .2  )  ) 


0049 


0050 
0051 


50  CONTINUE 


100  CONTINUE 

IF  (I0PT2.EQ.1.0R.I0PT3.EQ.1)   CALL  MAXEXT 


0052- 
0053 


RETURN 
END 
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MAXEXT 

This  subroutine  extracts  points  of  maxima  from  a  given  stress 
trace.  It  assumes  that  the  stress  values  are  time  wise  ordered 
at  a  constant  interval,  At.   It  further  assumes  that  no  local 
noise  is  induced  in  the  trace.   This  is  predicated  on  the  smooth 
behavior  of  the  stresses,  at  a  given  point,  generated  by  the 
structural  analysis  program,  as  opposed  to  the  type  of  data  ac- 
quired from  a  strain  gage.   Lastly  it  assumes  that  two  conse- 
cutive zero  values  of  stress  imply  that  the  end  of  the  trace  has 
been  reached,  which  is  consistent  with  the  steady  state  criteria 
used  in  the  structural  analysis  program.   The  slope  between  two 
points  is  calculated  as 

c   Q(D  -  a(M) 

&      (I-M) At 

if 

S  >  0 

then 

M  =  I, 
1=1  +  1 

and  the  test  repeated  until  S<0. 

If  S<0  then  a  test  is  made  to  determine  if 

la  (I)  -  a  (M)  I  >a_,    ,  ,  , 
1  '   Threshold 

Where  aml_    ,  ,  ,  is  established  by  a  user  as  the  minimum  decay 
Threshold  J 

in  stress  which  identifies  a  local  maximum.   If  this  threshold 

is  not  exceeded  the  next  stress  value  is  compared  with  the  most 

forward  value  just  tested  to  determine  if  decay  continues,  that 

is,  if 

a(I+l) <a(I) 
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thus  I  is  incremented  and  the  process  is  repeated.   If 

a(I+l)  >  a (I) 

then  both  M  and  I  are  updated  as  shown  above . 

If  the  threshold  is  exceeded  the  stress  class  interval  index,  N, 
for  the  i,j  histogram  is  determined  from 

N=  °M.f  1 

Aa 

where  N  is  an  integer  rounding  to  the  class  interval  containing 
a (M) .   The  class  interval  is  then  incremented  for  the  extracted 
stress  maxima  by  the  previously  determined  incidence  FREQ1,  that 
is 

a  (N,i,j)  =  a  (N,i,j)  +  FREQ1. 

This  subroutine  then  returns  to  the  calling  program. 

No  subroutines  are  called  by  this  subroutine. 
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f    MAXEXT           ] 

w 

m  =   1 

1 

IM  =    2 

1 

a      (IM)    -   o    (M) 

(Im-m)    At 

m  =  IM 
IM  =  Im+1 


NC  = 


_   P 


<5  (m) 


A6 


+  1 


<5  (NC,i,  .)=6  (NC,i,  .)+FREQl 


No 


No 


Figure  33.   MAXEXT  Program  Flow  Chart 
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Figure  33.   MAXEXT  Program  Flow  Chart  (Continued) 
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ORTRAN  IV  G  LEVEL   21 


MAXEXT 


DATE  =  73170 


00/31/27 


SUBROUTINE  MAXEXT 

INTEGER  SIGC 
-COMMON  /BLK2/  DT, 
1    NAP,  NB,  ND,  NI,  I 
-2 NXPTSt-N&SUG^  1-04-*- 


THRSH, 
NSI,  NS 


3  DELW(IO),  DELX(IO),  EEI(10,1 

4  XSUP(IO),  YB(10),  YP(4),  YCL 
6  E<3,10),  EI(3,10I,  XCS(3,10) 
7 -— -F(-80r10)-t  MO  (-1-1,20,  10),  MDT( 


IE,  IEVENT,  IN,  10,  ISS,  ITER,      0 
UP,  STEADY,  T,  WEIGHT,  XLB,  YBCL , 

0),     EIL(IO),     IS(10),    NCS(IO)* 
(10),    ZCL(10,10),     ZCU(10,10), 
,    ZBL(3,10),    ZBU(3,10), 
11,10),     WP(50,4),     X(80),XP(50) 


-0^0^ 


0005 


COMMON    AB44C5/-  S IG  ( &a)-+.DRUl l&O  U 

1  SIGP(120,3,10),    SIGC(10,3, 

JL -,  I0PT1,  10PT2, 

3      ,  IPTS(3,10) 


DFL(80,10),ITIM£,ICTIME, 
10),  SIGTHS,  DSIG 
I0PT3,  I0PT4,  NPTS(IO),  PTS(3,10) 


DIMENSION  ITST(3,10),  S IGMIN( 3, 10 ) ,  TIMMIN(3,10) 


THIS  SUBROUTINE  EXTRACTS  THE 
SXR£S.S-  TRACE  GIVEN  A  MINIMUM 


MAXIMA  STRESSES  FROM  A  DIGITAL 

DROP  IN  STRESS  FROM  THE  LOCAL  MAXIMA 


0007 
0008 


0009 
CLO10- 


JLF_ 
DO 

__D£L 


_(I.TIREJ,N£.l) 

5  J=l,10 

5_K  =  1  ,3 


ITST(KtJ)  = 
CONTINUE 


0 


GO  TO  6 


oaii 

0012 
0013 


0014 


I F„  i  1 0 P T 2  . EG. .0 )   GO 

WRITE  (10,1000)  IE, 
1000  FORMAT  (1HUJSJRESS . 


TO  7 

IEVENT 


LOAD  EVENTi , 14,//) 


7  CONTINUE 


0015 
OQlfj. 
0017 
0018 


0019 


DO  100  J=1,NB 

NPT  =  NPTS(J) 

IF  (NPT.EQ.O) 

J3£)_95  J1=1,NPT 


GO  TO  100 


IF  (  IOPT2.EQ.0)   GO  TO  10 


0020 
00-2L 


0022 
0023 


OG24- 


WRITE  (10,1010)  J,  PTS(K,J) 
1010  FORMAT  (1H0,'BEAM  NO. » ,  14,  • 

1      'STRESSES:  MINIMUM    TIME' 


POINT' ,F10.3,//1X, 
,  10X, 'MAXIMUM    TIME'  ,12X,« RANGE « ) 


1012 
10 


002  5 

0026 

0027 

0028- 

0029 

0030 


FORMAT  (1X,2(7X,F10.2,F7.3),7X,F10.2) 
CONTINUE 

START  THE  SEARCH 


lM-s-2 -    --    -— 

M    =    1 

CONTINUE 

IF     ( IM.GT.KTIME )       GO    TO    90 
-S- =     (SIGP(IM,K,  J)     -    SIGP(M,K,J) 

IF     (S.LT.O.)       GO    TO    1 
-4 F     (I  TS-T-IK^)  .HE  .  0 )-   -  GO    TO    4 

186 


ORTRAN  IV  G  LEVEL   21 


MAXEXT 


DATE  =  73170 


00/31/27 


0031 
0032 
0033 
0034 


IF(S.LT.SIGTHS)   GO  TO  15 
SIGMINIK,.)  =  SIGPlM,KtJ) 
TIMMIN(K»J)  =  DT*( ITIME-1+M) 
ITST(K,J»-  =  1  - 


-0G3S- 


0036 
00-37 
0038 


0039 
0040 


0041 
-00-42- 

jQ0A3— - 

0044 


4    GONTINUE- 


M    =     IM 
^M_=     IM  +    1 
GO    TO    2 


CONTINUE 

IF(ITST(  fcU  J.L_JE__U_CU  - 


-GO-TO    4 


0045 
0046 


0047 


0048 
0049 


SABS    =    ABS(SIGP( IM,K, J)    -    SIGP(M,K,J)) 
LE_  JSABS.GJ.SJGTHS)       GO    TO    3 

I£_1SIGP(1M  +  L»K^J).GE.SIGP(IM,K,  J))       GO    TO   4 

15    CONTINUE _ 


IM    =    IM    +    1 
GIL.TJD_.2_      ._ 

__.. U_PDATE _THE    APPROPRIATE    HISTOGRAM    STRESS    CLASS    INTERVAL 
_3________=_____1G__I_____K__J_)_ZDS,IG_ .+  __._._ _._ 


IF     (NC.GT.10)       NC    =    10 
SIGC.NC.K. J)     =_SIGC(NC,K, J)    +    1 


0050 
0051 
0Q52 


0053 
QQ_L4. 
0055 


80 


_LF_  (I0PT2.EQ..0)   GQ  TO  80 
SIGMAX  =  SIGP(M,K,J) 

SIGRNG  =  SIGMAX  -  SJGMJN.(_K_-JJ  ... 

TIMMAX  =  DT  *  (ITIME  -  1  +  M) 

WRITE  (10,1012)  SIGMIN(K, J) ,TIMMIN(K,J),  SI GMAX, TIMMAX, S IGRNG 

CONTINUE 


0056 
0057 


ITST(K,J)  =  0 
_GiX_XQ__4_ 


0058_ 

0059 

0060- 


90 

95 

1Q0 


CONTINUE 
CONTINUE 
CONTINUE 


0061 

RETURN 

0062 

END 
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INCRMT 

The  INCRMT  routine  updates  the  positions  of  the  axle  loads 

XPOSX.   =   XPOSX.   +   DXP0Sxi 

1  1      — 

NI 


and  updates  the  time 

T   =   T   +   DT 

If  the  debug  output  is  specified ,    the  axle  positions  and 
distance  traveled  over  the  load  interval  are  printed  out 

No  subroutines  are  called  by  this  routine. 
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INCRMT 


i  =  1 


& 


±- 


X. 

X.  (t+  t)=X.  (t)+  -~ 
1        1     N-J- 

i  =  i+1 


No 


Figure  34.     INCRMT  Program  Flow  Chart 
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:ORTRAN  IV  G  LEVEL   21 


INCRMT 


DATE 


73170 


00/31/2" 


0001 


0002 


000  3 


0004 


C 

c 


SUBROUTINE  INCRMT 


THIS  SUBROUTINE  INCREMENTS  THE  AXLE  LOADS 


REAL  MD,  MOT 


COMMON  /BLK2/  DT , 
-1 

2 
-3- 

4 
_6_ 


i  ,di_in^/  u.,       THRSH,   IE,  IEVENT,  IN,  10,  ISS,  ITER, 
NAT?.,  NB,  ND,  NI,  NL,  NSI,  NSUP,  STEAOY,  T,  WEIGHT,  XLB ,  YBCL, 
NXPTS,  NDBUG(IO), 

OELW(IQ),  0£LX( 10),  EEI(10,10),  EIL(IO),  IS( 10) ,  NCS( 10 > , 
XSUP(IO),  YB(10),  YP(4),  YCL(IO),  ZCL(10,10),  ZCU(10,10), 
_E(3,10i-»-Ea.(3,  10U  XC,S(3«-X0U.  ZBL13,10),  ZBU(3,1Q), 
F(80,10),  M0(ll,20, 10),  M0T(11,10),  WP(50,4),  X(80),XP{50) 

COMMON  /BLK3/  SUMHRX,  DTL,  IEVNTX,  NA ,  NTRUKX,  LAST, 
L  _  LTYPE(20),  WGTX«20),SPDX(20),  LLANE(20),  TIMEX(20), 
2      XP0SX(50),  LNUMXI50),  WEITX(50),  0XP0SX( 50 ) , ACCLRX ( 50 ) 


0005 
0006 
0007 
0008 


DO  1  11=1, NA 
_  XP.QLSXC  I  I1  =  XPQSX  (  1 1  )  +DXP0SX  (  1 1  )  /NI 
1  CONTINUE 
___I=L+DT 


0009. 


IF     INDBllGllOl.EO.O) G_CL  TQ    3D 


0010 
HQ.ll- 

HQ12. 


0013 
0014 
0015 


WRITE  (10,1001)  NA,NI 
1001  FORMAT UQ INCRMT    NO  AXLES  =',14,' 

1     •      X  DX  •  ) 

1020  FORMAT  (1X.,2F  12.4) 

DO  20  11=1, NA 
WRITE  (10.1020)  XPOSX( H)»DXP0SX( II ) 


NO  INCREMENTS  =• ,14,//, 


0016 

0017 
QQ18 


20    CONTINUE 

30    CONTINUE 

RETURN 
-EAia 
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MATINV 

The  MATINV  subroutine  computes  the  inverse  of  a  general  real 
Matrix  A  of  order  N. 

The  method  used  is  a  Gaussian  elimination  with  partial  pivoting 
used  to  obtain  a  first  estimate ,    B  ,  to  the  inverse  of  A.   Cor- 
rections may  be  performed,  as  an  option,  using  the  iterative 
scheme 

B  ,_  =  B  (2I-AB  ) 
r+1     r      r 

where  I  is  the  identity  matrix  and  C  =2I-AB   is  the  correction 

-1  r      r 

matrix  needed  to  calculate  the  B  .-ith  estimate.   At  each  step, 
the  one-norm  of  C  -I  is  computed,  i.e., 


[C] 1   =   maxi 


The  user  actually  specifies  only  a  key  indicating  whether  he 
does  or  does  not  want  the  corrections  performed.   A  maximum  of 
three  corrections  will  be  performed  since  no  further  improve- 
ment can  be  expected.   However,  if  the  condition 

|[2I-ABr+1J  r|2I  -  AB^j  <2"45  [21-mJj. 

is  satisfied,  the  iteration  stops  and  B   ,  is  returned  as  the 

final  inverse.   This  condition  indicates  that  the  machine  has 
done  as  well  as  it  can  in  obtaining  an  inverse  and  further 
iterations  will  not  improve  the  result. 


n 

\ 

E 

M 

-i 
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The  calling  sequence  is: 

CALL  MATINV  (NDIM, N, A, B,C,D, MAXIT, IERR) 

where 

NDIM  —  the  formal  row  dimension  of  A,  B,  and  C  in  the  call- 
ing program  (integer) . 

N  —  the  order  of  the  matrix  to  be  inverted,  N  <  NDIM 
(integer) . 

A  —  matrix  to  be  inverted,  whose  row  dimension  must  be 
NDIM,  and  whose  column  dimension  must  be  at  least 
N.   This  matrix  is  not  destroyed  (real) . 

B  —  inverse  of  A,  whose  row  dimension  must  be  NDIM  and 
whose  column  dimension  must  be  at  least  N  (real) . 

C  —  work  matrix  whose  row  dimension  must  be  NDIM  and 

whose  column  dimension  must  be  at  least  N  (real) . 

Note  that  if  corrections  are  performed,  C  contains 
2I-AB  where  B  is  the  returned  inverse. 

D  --  vector  of  working  storage  of  at  least  N  locations 
(real) . 

MAXIT  --  correction  indicator  (integer) . 

=0  —  no  corrections 

^  0  —  maximum  of  three  corrections 

On  return,  MAXIT  contains  the  number  of  corrections 
performed. 

IERR  —  error  return  (integer) . 

0  —  normal  return 

1  —  zero  pivot  found,  matrix  singular 

2  —  norm  of  correction  matrix  greater  than  1,  ma- 

trix may  not  be  convergent,  return  last  inverse 
as  estimate.   The  user  should  carefully  check 
the  inverse  for  this  error  return. 
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No  subroutines  are  used  by  this  routine 
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MATINV 


SET  B  TO  IDEN- 
TITY MATRIX 
SAVE  [A]  IN 
[C] 


1L 


BEGIN  TRIANGULARIZA- 
TION  OF  [C]  AND 
TRANSFORMATION  OF 
[B] 


&- 


PICK  UP  MAX  COL. 
ELEMENT  OF  [C]  ON 
OR  BELOW  DIAGONAL 


INTERCHANGE 
TWO  ROWS  OF 
[C] 


Jsl 


INTERCHANGE 
TWO  ROWS  OF 
[B] 


Figure  35.    MATINV  Program  Flow  Chart 
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A 


IERR=2 


^L 


BACK  SUBSTITUTE  TO 
GET  FIRST 
APPROXIMATION  OF 
THE  INVERSE 


B 


L 


No  J 


1  =  1 
MAXIT  =  3 


4^ 


PERFORM 
CORRECTION 


COMPUTE  NEW 
INVERSE 


1=1+1 


Figure  35.    MATINV  Program  Flow  Chart  (Continued) 
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ORTRAN  IV  G  LEVEL   21 


MATINV 


DATE 


73170 


00/31/2"; 


0001 


SUBROUTINE  MATINV  (NDIM,N , A, B, C, D, MAX  I T, IERR ) 

G-   —..-_— 

C    SUBROUTINE  TO  COMPUTE  THE  ITERATIVELY  REFINED  INVERSE  OF  A  GENERAL 
£ REA   MATRIX 

C       REAL  MATRIX.   TGE  FIRST  APPROXIMATION  IS  OBTAINED  BY  STANDARD 
-G GAUSS  IAN -ELIMINATION  W4-TH-P-ARTI AL  PIVOTING __...__  .     


€ i*jpy-T-PAftARET€RS 

C       NDIM  -  THE  FORMAL  ROW  DIMENSION  OF  A  IN  THE  CALLING  PROGRAM 

C N ^_-S-W€-43F  MATRIX  XO  BE- INVERTED 

C       A     -  MATRIX  TO  BE  INVERTED,  DIMENSION  AT  LEAST  A(NDIM,N). 
C A  IS -NOT  DESTROYED . 


C       B       INVERSE  OF  A,  DIMENSIONED  AT  LEAST  B(NOIM,N)  IN  THE 
X CALLING  PROGRAM  .. 

C       C       WORK  MATRIX,  DIMENSIONED  AT  LEAST  C(NDIM,N)  IN  THE 
X 1  CALLING  PROGRAM 

C       D       VECTOR  OF  WORKING  STORAGE  OF  LENGTH  AT  LEAST  N 
_C MAXIT-  ITERATION  -KEY .   


0  -  NO  ITERATION 

1  -  THREE  ITERATIONS 


_C OUTPUT  ...PARAMETERS 

C       B     -  INVERSE  OF  A 
X MAX  IT-  NUMBFR  QE_  I  TFRATinNS_PER£QRMED_«,-IF_.L.ESS  THAN  THREE. 


C  IMPLIES  NORM  OF  RESIDUAL  MATRIX  CONVERGED 

X LERR  -__.ERR OR  RETURN     .._._...   ..... 

C  0  -  NORMAL  RETURN 

X 1  -  MATRIX  SINGULAR  (ZERO  PIVOT  ELEMENT  FOUND) 

C                   2  -  NORM  NOT  LESS  THAN  ONE.   ERROR  MATRIX  MAY  NOT  BE 
X CONVERGENT.  DO_NOI  PERFORM  ANY  MORE  REFINEMENTS 


0QQ2 


000  3 

0004 
0005- 


oaoi>_„ 

0007 
0008- 


0009 


.DIMENSION  A  (NDIM,  1),  B(NDIM,1),  C(NDIM,1),  D(l) 
REAL*8  DP1,DP2 

NM  =  N  -  1 

IERR-J =— Q 


C   .SET  IDENTITY  MATRIX,  SAVE  A 
C 

DO  5  1=1  ,NM 
IJ  =  I  +  1 
CLLO-l _=_ A( J  ,Xi 


B(It  I  )  =  1. 


0010  DO  5  J=IJ,N 

XKLLJ C(L,JJ  =  A(1,J) 

0012               C(J,I )  =  A(J,I ) 
0013 B(I,J.t  -=_  JL. 


0014 

0015 
00-L6- 


-C- 


5  B(J,I  )  =0. 

C(N,N)  =  A(N,N  ) 
-BIN^N)  =.1. 


X BEGIN  TPIANGHLART7ATTnN  OF  C.    AND  TRANSFORMATION  OF  B 
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OUTRAN  IV  G  LEVEL   21  MATINV  DATE  =  73170 

C 

0Q17   ...  ... 00  100  J=1,NM 

0018  K  =  J 

00J.9 _  -_.T  =  ABS(C(J,J)) 

0020  LL  =  J  +  1 
C- 

C    PICK  UP  MAX  COL  ELEMENT  ON  OR  BELOW  DIAGONAL 
i C~  

0021  DO  20  I=LL,N 

0022- -P--=-  ABS4C-U  ,  J  )  ) 

0023  IF  (T.GE.P)   GO  TO  20 

HQ24 I__=l_P ._ ... 

0025  K  =  I 

on? 6  .20.  CONTINUE 

0027  IF  (K.EQoJ)   GO  TO  35 
C .    ._   _.. 

C    INTERCHANGE  TWO  ROWS  OF  C 

0028  DO  25  I=J,N 

jQQ22 T  .=.  C  <  vb.I  >. 

0030  C(J,I  )  =  C(K,I ) 

QQ3J. 25  C(K,  I)  =  T 

C 
C___lNI£fiCJdAWGE._'LMD_RQHS_QE_.B.  ...  .  _  . 

C 

QQ3.2 DP  30  1=1  ,N_ 

0033  T  =  B(J,I  ) 

0034. .  ._B  ( J  ,  I  )  =  B  ( K  ,  I  )   _ 

0035  30  B(KfI  )  =  T 

0036 3.5__DJ?1  .=  _C-UjlJ1  _ 

C 
.  C.  CHECK  FOR  SINGULARITY 

C 
003.1  ..   ...         IF  (.DP1.EQ.0.D0)   GO  TO  340 
0038                DO  100  I=J,NM 
003Q K  =  I  ,±JL.    .- 

C 
_   .  C    DP2/DP1  IS  THE  MULTIPLIER 

C 
004jO  DP2  =  CtK,J) 

C 
-0041 DO  -95_1J_=LL*JVI .  . 

0042  95  C(K,II)  =  C(K,II)  -  ( DP2*C( J, I  I  ) ) /DPI 
C 

0043  DO  98  11=1, N 

0Q44  98  B(K,II)  =  B(K,II)  -  ( DP2*B ( J , I  I  ) ) /DP  1 

C 
00-45- 100  C(JC,  J-L_=  -DJ>2/ DPI 

C 

C    END  TRIANGULARIZATION  OF  C  AND  TRANSFORMATION  OF  B 

C    CHECK  FOR  SINGULARITY 

C   - 

0046  IF  (C(N,N).EQ.O. )   GO  TO  340 

C -.. 

197 


00/31/2-! 


ORTRAN  IV  G  LEVEL   21 


MATINV 


DATE 


73X70 


00/31/2"/ 


0047 
00A8 


BEGIN  BACK  SUBSTITUTION 

DO  150  11=1, N 
-D(NJ  =  B(N,  II)/C(N,N) 


C049 


00  50 
0-051 
0052 

0053 
0054 


0055 


0059 


DO  140  I=4,NM  - 

0P1  =  O.DO 

_  K  =  N  -I 
L  =  K  +  1 

DO  130  J=LtN 
DP2  =  D(.)> 


130  DPI  =  DPI  +  DP2*C(K,J) 

0056  DP2  =  C(K,K) 

0057 L4Q  .QJUU    =-  .IB  ( K  ,  1 1  )  -  DP  1 )  /DP2 

C 
0058 DO  150  .l  =  L»N 


150  B(J, II  )  =  D<  J) 

END  BACK  SUBSTITUTION 
_.  ME_NDJfl!  HAVE  FIRST  APPROXIMATION  OF  THE  INVERSE 


D060 


JJ_1MAXIIJ.E_QJ0J._  GO  TQ.  350. 


0.0_61_ 
0062 


006  3 


OLDEL  =  0.0 
MAXIT  =  3 

PERFORM  THREE  CORRECTIONS 


DO  320  K0UNT=1, MAXIT 


0064 
Qfl65_ 
0066 

006  7 
0068 
0069 


-_ C- 


DO  240  1  =  1, N 
DO  240  J=1,N 
DPI  =  0. 

DO  230  K=1,N 
DP2  =  A( I,K) 
230  DPI  =  DPI  +  DP2*B(K,J) 


0070 
0071- 

0072 


0073 

0074- 

0075 


0076 
0077- 


C 
C 
C 


IF  ( I  .NE.J )   GO  TO  240 

D  Pi_  =  DP  1  -  2 .  - 

240  C(I, J)  =  -DPI 

NOW  C  =  (2*1  -  A*B) 

DEL  =  0. 

DO-  260  I=1,N  

DPI  =  0. 

COMPUTE  ONE-NORM  (MAX  ROW  SUM) 

DO  250  J=1,N 
15.0-.  DPI-  =  -DPI  +  A£S(C(  I  ,J) ) 
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007  8 
0079 
0080 
0081 


TEMP  =  DPI  -1. 
IF  (ABS(TEMP).LE.DEL) 
DEL  =  ABS(TEMP ) 
260  CONTINUE 


GO    TO    260 


CHECK    SUFFICIENT    CONDITION    FO    R    CONVERGENCE 


0082  IF     (DEL.LT.l.)       GO    TO    270 

0O83- 4-ERR   =    2 

0084  WRITE     (10,1003)     IERR 

-OOS-5 1003    F0RMA-H-'QE&&QPv    • .,  1 4^'-MA  TR  I X     INVERT       -    NORM    NOT    FOUND.    ERROR     •, 


OQ86 


0087 


1     'MATRIX   MAY  NOT  BE  CONVERGENT') 
GO  TO  350 


C 
-C- 

c 


-CHECK. FOR  NORM  CONVERGENCE  EPS 


2*#(-46) 


270  IF  (ABS-^D4L--0L-DE4-4->G-Ty(.14  2E-13)  *OLDEL)   GO  TO  280 


0088 
0089- 


0090- 


c 
c~ 

c 


LL  =  KOUNT  -  1 
G-Q-  TO  345  -  - 

COMPUTE  NEW  INVERSE 


280  00  310-I--WU- 


0091  DO  300  J  =  1,N 

QQ9.Z-^ DPL  =  0. 

C 

009-3 --DO  290  K=l  ,N 

0094  DP2  =  B( I,K) 

0095 290  DPI  =  _DPJL-±-JlP2*ClK»,J.L__ 


QQ96 


300  D(J)  =  DPI 


QQ9J 

0098 


0099 
0100 


DO  310   J=1,N 
310  B<  I  ,J)  =  D( J) 


320  OLDEL  =  DEL 
GO  TO  350 


0101 
0102 
-040-3- 
0104 


340  IERR  =  1 

WRITE  (10,1002)  IERR 
-4002  -FORMAT 4 'OERROR  UI4,  'MARRIX  INVERT   -  MATRIX  SINGULAR') 
GO  TO  350 


0105 

0106 
0107 


345  MAXIT  =  LL 

350  RETURN 
E-ND 
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MATMUL 

The  MATMUL  subroutine  computes  the  product  of  the  two  conformal 
matrices  A  (L,N)  and  B  (N,M)  and  stores  the  product  in  the  ma- 
trix C  (L,M) .   Calling  sequence  is:   CALL  MATMUL  (A,3,C,L,M,N, 
LL,NN)  where  LL  is  the  formal  dimension  of  arrays  A  and  C  and 
NN  is  the  formal  dimension  of  array  B. 

No  subroutines  are  called  by  this  routine. 
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f    MATMUL 


&L 


1  =  1 


2& 


J  =  1 


31 


C(I,J)=0.0 


u. 


K 


-M 


COMPUTE  C(I,J)  = 
C(I,J)+A(I,K)*BCK,J) 


Figure  36.    MATMUL  Program  Flow  Chart 
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00/31/2T 


>0001  SUBROUTINE  MATMUL  (A,  B,  C,  L,  M,  N,  LL,  NN ) 

! C 

C         COMPUTES    THE    ELEMENTS    OF    C(L,M)     AS    A    REAULT    OF    THE    PRODUCT    OF 
i _- C A1L»N)_AMLJJN^M  )__1N  -JLHAT-ORDER,. 

C         LL     IS    THE    FORMAL    DIMENSION    OF    ARRAYS    A    AND    C 
X NN    IS    THE    FORMAL    DIMENSION    OF    ARRAY    B  


0002 


DIM  FN  ST  ON-  A(LL»_L)^  &INN.»1 1*  _C.LLL, 11 


14X003- 


]  0004 

loxios- 


DO   5    I=UL- 


DO    5    J=1,M 

r.cit.n  =  n.o 


0006 


0007 

nfioa_ 


10009 

loom 


on   5   K=lUM 


CU,J>    =   CU,J>    +  AU,K)*B(K,J> 

i  INtJF . 

RETURN 

FND 
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